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Introduction 



The aim of Essays on Microeconomics and Industrial Organisation (2" d edition) 
is to serve as a source and work of reference and consultation for the field of 
Microeconomics in general and of Industrial Organisation in particular. 

Traditionally, Microeconomics is essentially taught as theory and, although 
handbooks illustrate the various microeconomic theories with examples and 
practical cases, they hardly ever offer an estimation of a demand, production and 
cost fimction. In fact, Microeconomics is explained with self-contained theories 
without empirical tests. The editor of these Essays has taught Microeconomics for 
twelve years in the traditional way; for the last ten years -in advanced courses and 
doctorates- he has offered a selection of empirical applications, which have 
complemented the traditional theoretical teaching. These applications have 
emerged from various research projects managed by the editor during the last ten 
years with the financial support of several institutions (DGICYT, DGESIT, 
CICYT, and R&D National Plan). The success in this type of teaching and the 
availability of recent original applications from authors usually collaborating with 
the editor has led the later to compose this text. 

This text combines microeconomic theories with appropriate empirical tests. 
The standardised microeconomic analysis of demand, production and costs 
(supply) is set forth along with appropriate econometric techniques. Moreover, it 
should be pointed out that over the last two decades Microeconomics has greatly 
broadened its field of application. On the one hand, this has been due to the fact 
that the conditions required for existence, unicity and stability of the general 
competitive equilibrium have been met. This was the prevailing focus of 
Microeconomics in the Sixties and part of the Seventies. On the other hand, as big 
samples of inter- and intra-industrial data were increasingly available, a neo- 
classical Microeconomics branch emerged in the mid Seventies traditionally 
called in Britain New Microeconomics or Industrial Organisation, which on the 
structure-behaviour-results paradigm, binds together the earlier and new works on 
structural change and technical progress and applies new techniques -mainly panel 
data- which enable us to observe how the behaviour of the new agents affects 
industrial structure. 

Demand, production and costs are parts of Microeconomics which are greatly 
active at present. Industrial structure and regulation, markets and failures of 
market constitute central nuclei of Industrial Organisation. Therefore, the second 
part of the title of these essays records this expression. Although I dare not give it 
a Roman numeral, a new volume will foreseeably emerge in the fiiture recording 
new advances both in the expression of ideas and in the econometric cointegration 
and panel data techniques used here. 
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Introduction 



The text consists of four parts: Demand, Production and Costs (Supply), Market 
and Industrial Structure and Failures of Market and Industrial Regulation. Each 
part has three chapters. 

Section I deals with demand and Starts with a paper that studies industrial 
demand. Chapter 2 offers a study on air transport demand with respect to the 
remaining modes of transport. Chapter 3 presents the behaviour of the consumer’s 
with respect to the introduction of a new product (High Speed Train AVE). To 
conclude, this section presents an approach to the hub-and-spoke Systems from 
SVARs models, with a practical application to Container traffic between the port 
of Bahia de Algeciras and other ports of the Spanish port System. 

Section II deals with supply. In Chapter 5 a production function is estimated 
with panel data for the hospital industry and the results corresponding to 
asignative efficiency are presented. In Chapter 6 economic efficiency is analysed 
for the case of road transport firms. In Chapter 7 technical efficiency is analysed 
for the case of international air transport (1992-2000). Finally, Chapter 8 offers the 
productor’s behaviour with technological innovation and structural change. 

Section III studies the market and industrial structure. Chapter 9 covers the 
Problems of measurement of inter- and intra-industry trade. Chapter 10 presents a 
categorical aggregation exercise with Spanish trade data. The following chapter 
presents a typical model of Industrial Organisation in which intra-industry activity 
is explained by industrial structure and the behaviour of the agents. To conclude, 
this section presents a paper on economic integration in which the link between 
structural adjustment and horizontal and vertical intra-industry trade is analysed. 

Section IV Starts with a chapter on Organisation and regulation of the port 
industry with an application to Europe and Spain. Chapter 14 presents an 
application of positive theory of regulation to Spanish foreign trade. Chapter 15 
studies the structure and regulation in electrical industry. In Chapter 16 the effects 
of a regulation of working time on labour market are studied. Chapter 17 studies 
transitional dynamics and endogenous growth for the case of public Capital. 
Chapter 18 analyses the relationships among unions, wages and productivity. 
Chapter 19 includes a study on economic impact in Spanish ports. Chapter 20 
studies economic impact in an airport. Chapter 21 includes a theoretical model 
which explains how the lack of co-ordination in the input and output of agents in 
industry may generate complex situations and presents an application. Chapter 22 
describes a theoretical model that allows to calculate congestion effects on 
transport including an application. Finally, in Chapter 23 a cost-benefit analysis 
(CBA) of High Speed Train (AVE) Madrid-Barcelona is offered. 




PART I. DEMAND 




1 Modeling Seasonal Integrated Time Series: the 
Spanish Industrial Production Index 



J. L. Gallego-Gömez 
University of Cantabria (Spain) 



In this paper the Box-Jenkins approach to the building of seasonal time series 
model is extended so that it is adequate to model seasonally integrated time series. 
To this end, the dass of multiplicative ARIMA models is broadened in such a 
way that it allows to describe time series integrated at a few of the seasonal 
ffequencies. Thus, tests for seasonal unit roots are not considered as a rival 
modeling approach, but can be used in the identification stage to decide the 
transformation inducing stationarity. The fit model is used to generate forecasts 
and to estimate unobservable components. The enhanced Box-Jenkins approach is 
illustrated modeling the Spanish Industrial Production Index. 

1.1 

Seasonal Time Series Models 

The Box and Jenkins approach to the seasonal time series analysis can be sketched 
as follows. Firstly, the non-seasonal 1 - B and seasonal 1 - B s differencing 
operators are used to convert a non-stationary series z t into a stationary series w t . It 
is usually necessary to use d-order non-seasonal and D-order seasonal 
differencing, that is, 

w t =(1-B) d (l-B S ) D 

Then, the stationary series w t is expressed, according the Wold decomposition 
theorem, as a weighted sum of current and past values of a white noise process 

QO 

w t =a t + 'Lys j a t _j = y/(B)a t 
j = i 



where 

P. Coto-Millan (ed.), Essays on Micro economics and Industrial Organisation 
© Physica- Verlag Heidelberg 2004 




J. L. Gallego-Gömez 



ys(B) = i + TvjB J 
7=1 

Finally, to achieve parsimonious models the polynomial y/(B) is approximated by 
the rational polynomial 

<f>(B)0(B s ) 

where <f(B) = 1 - faß - ... - (j)pB p and 6{B) = 1 - 0 X B - ... - 6 q B q are the non- 
seasonal autoregressive and moving-average polynomials which describe the 
dependence between consecutive data, and @(B) - 1 - @ x B - ... - & P B P and 
0(B) = 1 - 0] 5 - ... - 0qB q are the seasonal autoregressive and moving-average 
polynomials describing the dependence between data which are s periods apart. 
Therefore, the nonstatinonary seasonal time series z t is described by the general 
multiplicative model 

^(ß)<Z>(5 s )V rf Vfz r = ß(B)0(B s )a t (L1) 

The choice of the seasonal differencing to induce stationarity is based on the fact 
that seasonal time series show a cyclical behaviour with period s (s = 4 for 
quarterly data and s= 12 for monthly data). In the extreme case z t = z t . s , we have a 
homogeneous difference equation whose general solution is a linear combination 

of the 5 roots, e l2nj/s (j = 0, 1, . . ., s-1 and i=4~ 1 ), of the characteristic equation 
1 -^ = 0 , 

5—1 . [ 5 / 2 ] y-k 2 irk 

z t = Y,a i e llni ’ s =« 0 + Z («k cos (—t) + ß k sin(— 0) 



2 7zk 2 Tck 

where a k cos( 1) + ß k sin( 1 ) is a harmonic oscillation with period p k = 

s!2n k and constant amplitude A k = (a 2 + /7 2 ) 12 , generated by the pair of conjugate 



complex roots e l 



-iln ( k-s)/s . 



the harmonic with frequency s/2, a 5/2 cos {n t), 



only arises when s is even and is generated by the negative unit real root; [s/2] = 
s/2 if 5 is even and [s/2] = (s-l)/2 if s is odd. Hence, it is seen the key role played 
by the seasonal difference in representing the seasonal pattem. 

Although the 1 - Bs filter has been commonly used in the Box-Jenkins 
approach, Abraham and Box (1978) pointed out that sometimes some of its factors 
could be sufficient to handle the seasonality. Such factors are found by expresing 
1 - Bs in terms of its roots 

\-B s =Y\(\-e i2Kjls B) 

7=0 

and joining each pair of complex roots into a second-order factor 



(l-e i2>ck/s B)(l-e 



il7c{s-k)l s _ 



ß) = l-2cos(2®t/ s)B + B l 
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Thus, the seasonal difference can be factored as follows 

[s/2] 

1 -B s = ns k (B) 

*=° (1.2) 

where 



S k (B)- 



1 -B 

1-2 cos(2;z£ / s)B + B 2 ) 
1 + B 



k = 0 

* = l,...,[(s-l)/2 

Jfc = [s/2] = s/2 



(1.3) 



Hylleberg, Engle, Granger and Yoo [HEGY] (1990) developed unit root tests to 
identify the simplifying operators S k (B) in the underlying time series model. These 
authors extend the notion of integration to cover the seasonal frequencies. Thus, a 
time series z t is integrated of order d at the frequency k, denoted by z t ~ I k (d), if its 
autoregressive representation contains the factor S k (B) d . Empirical evidence based 
on the HEGY test reveals that most of the monthly and quarterly economic time 
series are integrated of order one at a few seasonal frequencies (see, e. g., Osbom 
1990, Beaulieu and Miron 1993, Hylleberg, Jorgensen and Sorensen 1993). 

In contrast, the empirical evidence found using the Box-Jenkins methodology 
has revealed that most of the monthly and quarterly seasonal time series can be 
adequately described by the IMA(0,1,1)(0,1,1) 5 model, 

(1 - B)( 1 - B s )z t = (1 - 0B)( 1 - ®B S )a t 

implying that z t is I 0 ( 2) and 4(1) for k = 1, . . ., [s/2\. 

Gallego and Treadway (1995) have shown why the empirical results obtained 
with the HEGY test cannot be found with the conventional multiplicative ARIMA 
dass, and how to broaden such a family to allow for unit roots at a few seasonal 
frequencies. Consider for example the seasonal IMA(1,1) 4 model for quarterly 
time series 

(l-B 4 )z l =(l-©B 4 )a, (1.4) 

It is interesting to write (1.4) as 

(1 - ö)(l + B 2 )(1 + B)z t = {\-© xl4 B)(\ + © xl2 B 2 )(\ + © xl4 B)a t (1.5) 

The three nonstationary factors on the left-hand side of (1.5) contribute to 
"spectral peaks" at the frequencies 0, 1 and 2, respectively. So, under the condition 
of invertibility 0 < 1, z t is integrated of order one, z t ~ 4(1), at k = 0, 1 and 2. 
However, if the MA parameter © is positive and strictly noninvertible, 0 = 1, then 
the model contains three common factors whose cancellation implies that z t is 
integrated of order zero, z t ~ 4(0), at k = 0, 1 and 2. In contrants, let the seasonal 
MA(1) 4 polynomial in model (1.4) be replaced with a nonseasonal MA(4) 
polynomial. Then, the resulting IMA(1,4)(1,0) 4 model 

(l-B)(\ + B 2 )(\ + B)z, =(1 -G x B- 9 2 B 2 -0 3 B 3 -0 4 B 4 )a t 



( 1 . 6 ) 
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enables the description of a richer dass of integrated series. Here, the case z t ~ 
4(0) arises when the model contains a common factor with frequency k, that is, 
when 6{B) = 0*(B)S k (B). But, if the factor S k (B) is not included in the MA 
structure, then z t ~ 4(1)- In (1.6) the Order of integration of the process is not 
necessarily the same at all frequencies. In fact, z t ~ I k (d k ), being d k = 0 or 1 for k = 
0, 1, and 2. 

When investigating the Order of integration of quarterly time series, it is 
convenient to express the IMA(1,4)(0,1) 4 in such a way that it shows more clearly 
the presence of common factors. Thus, in the constrained IMA(1,4)(1,0) 4 model 

( l-B)(l + B 2 )(l + B)z t = (l-0 Q B)(l + 0 x B 2 )(l+0 l B)a t (1.7) 

where 0 < 6 k < 1, one can readily to see that z t ~ 4(1) if 0 < 0 k < 1 and z t ~ 4(0) if 
9* = 1. The MA(4) polynomial has been constrained so that it has roots at the 
frequencies k = 0, 1 and 2. This constraint is especially relevant when modeling 
monthly time series because the replacement of a seasonal MA(1 )j 2 polynomial by 
a nonseasonal MA(12) polynomial leads to a nonparsimonious model, whose 
estimation could be problematic. In such a case, the constrained MA(12) 
polynomial can be expressed as the product of the following seven factors 

(1 - 0oB), (1 - B + 0 X B\ (1 - 4^ B+ 02B\ (1 + 0iB 2 ) , 

(1 + -faß + 0 A B 2 ), (1 + 73ö7 B + 0 S B 2 ), (1 + 0(ß) 

where 0 < 0k < 1 . In general, a MA(s) polynomial can be constrained so that it 
has s roots at the frequencies k = 0, 1, . . . , [s/2], which correspond to the 
following factors with real coefficients 

0(B) = [ lfs k (B) (1.8) 

k = 0 



where 



$ k (B) = 



\-S 0 B 



1-2 cos(2;zfc / s)^9^B + 3 h B 2 ) 



1 + ^/ 2 ^ 



k - 0 

Jfc = l,-»,[(j-l)/2] 
k = s / 2 and s even 



(1.9) 



The subscript k of the coefficient 3 k indicates the number of cycles in a period of 5 
time instants. Comparing (1.8)-(1.9) with (1.2)-(1.3) it can be clearly seen the one- 
to-one correspondence between the factors 6 k (B) of the constrained MA(s) 
polynomial and the simplifying factors S k (B) of the seasonal difference. Thus, 
each MA factor acts as an indicator for overdifferencing. 

In summary, seasonal integrated time series can be described by the generalized 
multiplicative ARIMA model 

^(S)0(ß s )V' il n\(5yH = 0{B)®{B s )S{B)a t 

k= 1 



( 1 . 10 ) 
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which is compactly written as (p(B)z t = rj(B)a t . 



1.2 

Forecasting 

In this section the implications that the generalisation of the conventional 
multiplicative ARIMA model has on the forecasting fimction are examined; in 
particular, those derived from the relaxation of the assumption that time series are 
integrated at all seasonal frequencies. To this end, two specifications for the model 
(1.10) are used: (i) the difference equation representation 

Z t+l ~ ( P\ Z t+l-\ ( Pp Z tM-p ~ a t+l ~ r l\ a t+l-\ r lq a t+l-q 

which is useful to recursively compute point forecasts, and (ii) the specification in 
terms of innovations 



z t+i ~ a t+i +v\ a t+i-\ +y / 2 a t+i-i + ••• 



( 1 . 12 ) 



or 



z t+i ~ y / (B)a t+l - (1 + y/ x B + y/ 2 B 2 +---)a t+ i 

which is more convenient to calculate interval forecasts and to update forecasts as 
new data become available. 

The minimum mean square error forecast of z t+i at origin t, for lead time /, is the 
conditional expectation 2 t (/) = E[z t+l / z t , z t _ x ,•••]. Thus, taking conditional 
expectations at time t in (1.1 1) we obtain for / > q the difference equation 

2 t (l)-(p l 2t(l-l) (p p 2 t (l-p) = 0, l > q ^ 

whose general solution is 

2 t (l) = bl l) Ml) + b^f 2 (l) + - + bff p (l), l > max(< 7 - p,0) 

To compute the p coefflcients b\ ,) we need p initial conditions or pivotal 
values 2 t (r + 1) , 2 t (r + 2), . . ., 2 t (r + p) , where r = max(^-p,0). Let the column 

vector b (0 be the p coeficients b f ) , / =' ' 1, . . . , p, then 



b (0 =(ß T ¥y l F T z t 

where F = [f(j) ] is a p x p matrix of fitting functions, i t = [ 2 t (j) ] is an p x 1 
column vector, and j = r+1, . . . , r+ p. 

From ( 1 . 1 2), it follows that 

2t(i) - ¥i a t +v'i+i<*t-i + y / i + 2 a t-2 +••• 

Therefore, the forecast error / Steps ahead 
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^ (0 = a,+i + Wv a t + i-\ + — + Vi - 1 a M 

has zero mean and variance 

V[e t (/)] = cr 2 a {\ + ii/t +■■■ + y/h ) 



The coefficients y/j of the polynomial y/ (B) =rj{B)l(p{B) satisfy the difference 
equation 



Vj -<P\Vj-\ 9 P w j- P =o, j > q 

whose general solution is 

Wj = c x f x {j) + c 2 f 2 {j) + --- + c p f p (j), j>max(0,q-p) 



It is convenient to write 

\|/ = Fe (1.15) 

where V| / = (y/^ u . . . , y/ r+p ) T and c = (c u c 2 , c p ) T . The y/j weights are also used 
in updating the forecasts. As the forecasts of the future observation z t+i made at 
origins t and t - 1 are 

2 t (l) = y / ,a,+ij/ M a t_\+y l+ 2 a t-2 +••• 



2f-l (l + — Vl+\ a t-\ +Vl + 2a,-2 +^/+3«f-3 +- 



the updating formula is 
2f0) = %t - i d + l) + ^ /r / a ? 

Hence, it follows that 
Fb w =Gb (M) +v) ta t 

where G = [fj(j + 1 ) ],j = r+ 1 , . . . ,r+ p. Solving for b (,) we have that 
b w = (F r F)” 1 F 7 'Gb (?_1) +c a, 

Now we can see the role played by the constrained MA operators. While the roots 
of the autoregressive operator determine the form of the eventual forecast fimction 
and the pattem of the y/j weights, the moving average coefficients are relevant to 
provide initical conditions. Since that the replacement of the seasonal MA(1) S 
operator by the nonseasonal MA(s) operator does not change the Orders of the 
extended ARMA model, neither the form of the forecast function nor the number 
of pivotal values are modified. Obviously to the extent the y/j weights changes, the 
pivotal values will change and also the position of the the eventual forecast 
function. However, the main implications arise when some of the MA roots are 
close to unit and cancel with some common simplifying operators, that is, when 
the series is /*( 0) at a few k . In such a case, some coefficients c t will be zero and, 
therefore, the forecasting fimction will contain non-updated or deterministic 
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components and the squares sum of y/j weights, v[e t (l)]= cr^c r F r Fc , will be 
smaller. 

1.3 

Unobservable Components 

Several procedures have been proposed to estimate the trend T h seasonal S t and 
irregulär I t components of a time series described by an ARIMA model (see, e. g. , 
Gallego 2000). I addopt an variant of the approach based on the eventual 
forecasting function. Equation (1.14) can be written as 

^(0 = S^°/t(0+ ±b?f k {l) = t t+ S t (1.16) 

k = 1 k=p\+l 

where it has been assumed that the first pj fitting fiinctions, f k (l) (k= 1, ..., p\) 
are nonseasonal, and f k {l) ( k=p x +\ , ...,jp) are seasonal fimctions; f t and S t are 
the estimates for the trend and seasonal components. 

Assuming that the unobservable components are described by ARIMA models, 
then it follows that 

z t = 2t-\0) + a t> T t = f t _ x (\) + b n S t = ^_i(l) + G ’ It = It-\0) + dt 

where b h c h and d t are white noise innovations. Hence, it is clear that the 
estimation of each unobservable component involves the calculation of its one- 
step-ahead predictor and error. While the quantities T t _ x (1) , § t _ x (1) , and l t _ x (1) 
can be obtained easily by breaking down 2 t _ x (1) , the remaining problem is how to 

estimate the shocks b h c h and d t . 

For the additive seasonal model 

z t = Tt + S t + I t 

it follows that 

2 t _ x (l) = T t _ x (l) + S t _ x (\) + I t _ x (l) 

and 



a t = b t + c t + d t 

It is now apparent that the difficulty inherent to the decomposition of seasonal 
time series is how to isolate the three component innovations of the residual series 
a t . As in other approaches, the identification problem also arises here. 
Nonetheless, the dass of observationally equivalent decompositions is readily 
obtainable from the different allocations of white noise series among the one-step- 
ahead predictors 

T t =t t _ x (\) + e Xn S t = S t _ x (\) + e 2t , I t =l t _ x (\) + a t -e xt -e 2t 
where e Xt and e 2t are white noise process. Hence, the canonical decomposition is 
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T, c =T t -d 1), Sf = $ M (1), If =/ M (l) + fl r (1.17) 

The I t _ x (1) component will be zero for models with p < q, since none of the AR 
roots generates an irregulär movement. In such a case, the irregulär component is 
equal to the residuals of the model fit to the data. In contrast, for top-heavy models 
p < q the eventual forecast function for includes moving average terms, 

which are allocated to the irregulär component. This result is consistent with the 
Observation by Burman (1980) that top heavy models lead to irregulär components 
with moving average structure. 

1.4 

Empirical Analysis 

In this section, I embed the seasonal Integration notion into the Box-Jenkins 
iterative modeling strategy. Although this can be made in several ways I follow a 
"general to specific” approach by assuming that times series, when showing 
seasonality, are integrated of Order one in all of the seasonal frequencies. This 
general assumption does not imply a big change in the Box-Jenkins practice, in 
which the airlines models arises often as a tentative representation for many 
seasonal time series. Then, a more specific model is arrived at by looking for 
overdifferencing at single frequencies. The specific model, which can be thought 
of as a mixed regression-ARIMA model or transfer function model, is used to 
generate forecast and to estimate unobservable components. 

Fig. 1.1 displays several Identification tools for the Spanish Industrial 
Production Index from January 1986 to January 2000, 169 monthly observations. 

Folio wing the Box-Jenkins approach it is concluded that VV 12 z f seems to be a 
stationary series. On the other hand, its periodogram, column 3 in Fig. 1.1, peaks 
at the calendar frequency 0. 348 (close to 7t/3, indicating a trading day Variation 
which can be modelled (see, e. g. , Bell and Hillmer 1983) as 

TD,=j:ß i td it 

/=o 

where td it (i = 0, ..., 6) are, respectively, the number of Sundays, Mondays, ... in 
month t. Since lom t =Y^ =0 td it , where lom t is the length of month /, an 
alternative speciflcation for the trading day Variation is 

TD, = ß Q lom, + X ß*td*u 
1=1 

where ß 0 - (l/7)^ =0 td* it - td it -td Qt , and ß* = ß t - ß 0 measures the 

difference between the i-th day effect and the lom effect. Also it is expected that 
the production decreases due to holidays such as Easter. The Easter effect here is 
modeled as 
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E t =at, Easter 

where g E “' ,er is the proportion of Easter days falling in the month t. Here, it is 
assumed that the length of the Easter is four days. 




Fig. 1.1. Plot, Periodogram, Correlogram and Partial Autocorrelation Function for the 
series: z u Vz f y VV 12 Zf 



Since the sample autocorrelation function (SACF) and the sample partial 
autocorrelation function of VV 12 z, can be distorted by the calendar and Easter 
effects, instead of identifying a tentative multiplicative ARMA structure, the 
airline model with deterministic variables 

z t = ß’olom, + Iß-td* + a4f aster + N, 

i=l 

VV 12 N t = (1 - 0B)(1 - 0B 12 )a t 

is estimated as starting point. The exact maximum likelihood estimates obtained 
with the AMB program (Gailego 2000b) are 





14 



J. L. Gallego-Gömez 



z t = -6.93^f aster + 2.741om t + 0. 15td* t + 0.23td* 2t + 0.76td;, 
(0.68) (1.10) (0.33) (0.34) (0.32) 

+ 0.73td4, + 0.57tdj, - 1 ,35tdg t + N t 
(0.32) (0.33) (0.33) 



VV 12 N t =(l-0.47B)(l-0.83B 12 )at 
(0.07) (0.08) 

a = .054(.175),a a = 2.19,max(p ij ) = -.58, 

Q(39) = 58,LLF = -343.4, AIC = 1.69 (1.18) 

Some of the usual diagnostic checks reveal model inadequacies: (i) the g-statistic 
for the first n/4 residual autocorrelations is greater than the 5% point for % 2 with 
37 degrees of freedom, 52. 192320; (ii) the SACF and SPACF take large values 

compared with the ±2! 4n bands at some lags. While it is not clear how to 
reformulate this model with conventional operators, some improvement is obtain 
by replacing the MA(1 )i 2 i by a constrained MA(12) polynomial and including an 
AR(2) 12 polynomial to describe the annual dependence structure indicated by lags 
12, 24 and 36. The estimated model is 

z t = -7.23£,f aster + 2.211om t + 0.25td„ - 0.10td* 2t + 1.26td;, 

(0.53) (0.91) (0.26) (0.27) (0.26) 

+ 0.28td; t + 1.35td;,-1.95td; t +N t 
(0.26) (0.26) (0.26) 

(1 - .40B 12 + .21B 24 ) VV 12 N t = (1 - .34B) (1 - .99B) (1 - 1.71B + .97B 2 ) 

(.14) (.09) (.09) (.06) (.04) 

(1 - 1 .OB + 1 .OB 2 ) (1 - .97B 2 ) (1 + .97B + .93B 2 ) (1 - 1 .73B + ,99B 2 ) (1 + .96B)at 
(-.-) (.05) (.06) (.15) (.03) 

a = .033(.174),a a = 2.18,max( Pij ) = -.58, 

Q(39) = 50, LLF = -342.8, AIC = 1 .79 (1.19) 

where the presence of a common factor 1 - B + B 2 on both sides of the model 
equation is discovered. This is indicative of the presence of a single deterministic 
seasonal component with a period of six months or a frequency of two cycles per 





1 Modeling Seasonal Integrated Time Series 



15 



year. Thus, the common factor can be removed, but a pair of trigonometric 
variables c 2t and s 2t are added to the model (see,e. g. , Box, Jenkins and Reinsei 

1994). Furthermore, the simplifying factors 1 - B and 1 - a/3 B + B 2 of the 
seasonal difference are nearly cancelated by the associated MA operators, 
indicating that a constant term and a pair of trigonometric variables c u and c 2t 
must be included. Proceeding in this way we reach the estimated model 

Zf =-7. 06^ m,er +3.26 lom, +.06td* u +.29td* 2t +.66tdl, + l%td\, +.60 td\ t 
(0.54) (0.91) (.26) (.27) (.26) (.26) (.26) 

-\A0td* 6t +\21c Xt +5.90.S], +1.23c 2 , -7.73s 2 , -10.02c 3 , -2.03 s 3t 
(.26) (-31) (.40) (.29) (.22) (.46) (.16) 

+ 2.53c 4r +5.15s 4 , +1.72c 5; -7.66s 2t +N t 
(.16) (.40) (.21) (.74) 

(1 - .28 B u -¥ 23B 24 ) [V(l + B)N t - .40] = (1 - .365) (1 + .9 6B)at 

(.08) (.09) (.19) (.09) (.03) 

a = .033(.174),cr fl = 2.18,max(p iy ) = -.58, 

ß(39) - 50, LLF = -342.8, AIC = 1 .79 



( 1 . 20 ) 

Now the main diagnosis checks do not reveal serious model inadequacies: the Q- 
statistics is less than the 50% point of the % 2 with 40 degrees of freedom, and the 
residual plot, cumulative periodogram, correlogram and SPACF are consistent 
with the white noise hypothesis. 

Several implications are derived ffom model (1.20). Firstly, the time series z t is 
4(1) at k = 0 and 6, but 4(0) at k = 1, ..., 5. Similar results are found with the 
Beaulieu and Miron test for seasonal unit roots showed in Table 1. Note that these 
results can not be indicated by the airlines model. Secondly, the trend component 
is described by a linear trend with stochastic Ordinate, but deterministic slope. 
Thirdly, the seasonal component is the sum of five harmonic oscillations with 
deterministic amplitude and one oscillation with stochastic deterministic. 
Therefore, both the trend and seasonal component are a mixture of deterministic 
and stochastic terms. Finally, a damped cyclical component with period 4. 96 
years is described by the annual AR(2) 12 polynomial. The estimates of three 
unobservable components, following the method described in the above section, 
are shown in Fig. 1.2. 



Table 1.1. Test for seasonal unit roots 





ao fli 


b\ a 2 b 2 


Ö3 


b 3 a 4 




b 4 


a 5 


b 5 


a 6 


t 


-1.49 -2.69 


-3.80 -4.32 -248 


-4. 06 


-1.12 -3 


.29 


2.21 


-2. 93 


.20 


-2. 68 


F 


10.4 


14.4 . 


9. 10 


8. 


48 




4. 34 







Lags: 5, 6, 7, 12, 15, 32. Deterministic inputs: Linear trend, seasonal dummies and calendar 
effects. Diagnostic checks: & a = l. 8301, g(31) = 33. 01, LLF = -252. 91, AIC = 1. 8327. 
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Fig. 1.2. Estimated trend, seasonal and irregulär components for the Spanish IPI 
provided by model (1.20). 
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In this study, a theoretical model for the passenger transport demands in Spain is 
presented. Quarterly data have been used for the 1980.I-1992.IV period. 

Cointegration techniques, which are subject to a wide ränge of tests, are used to 
obtain short and long run equations. Moreover, the product, price and cross 
elasticities of each mode of transport are obtained. These estimations may be used 
to analyse the effects of transport tariffs on income changes, as well as to predict 
short- and long run traffics. 

2.1 

Introduction 

The initial models of passenger transport demand were the aggregate "modal split 
models". In these models, there has been an attempt to determine the number of 
joumeys by a given set of modes of transport for two towns, taking into account 
the characteristics of the passengers. Studies on modal split, such as Quandt and 
Baumol (1966), Boyer (1977), and Levin (1978), have been criticised by Oum 
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(1979) and Winston (1985), among others, for the few variables used to account 
for the motivation in the user’s behaviour, and for using very simple linear pattems 
in their estimations. 

Several models of aggregate passenger transport demand based on the user's 
behaviour have been carried out in Order to improve the previous ones. The user's 
utility is optimised in these models in line with the classic theory of consumer's 
behaviour and demand. The work by Oum and Guillen (1979) based on the user’s 
behaviour is a typical example in which the passenger demand in Canada is 
analysed. 

Some disaggregate research based on the user's behaviour has also been done 
on passenger transport demand. The most significant work on these models is 
McFadden (1973, 1974). In these works, the user takes a discrete choice of some 
of the different modes of transport (railway, air, road transport, etc.) and it is 
assumed that the mode chosen optimises the utility for the user. 

Spanish interurban passenger transport was first studied in the "Elasticidad de 
la Demanda del Transporte Püblico de Viajeros" (Elasticity of the Passenger 
Public Transport Demand) by the Instituto de Estudios de Transportes y 
Comunicaciones (Institute of Transport and Communication Studies) (1978). This 
was analysed by Väzquez (1985) in a work carried out by the Secretaria General 
del Ministerio de Transportes (General Secretariat of the Ministry of Transport). 
In addition, other studies such as that by Inglada (1992), Coto-Millän and Sarabia 
(1994) and Coto-Millän, Banos and Inglada (1997) have been carried out on this 
issue. The elasticities of the modes of transport in the main regions were studied in 
IETC (1978) and Väzquez (1985). Price elasticities have been studied in Inglada 
(1992) for monthly data between 1980.01 and 1988.12, with time series in which 
the residues have been modelled with the Box-Jenkins techniques. Uniequational 
models have been carried out in Coto-Millän and Sarabia (1994) in Order to 
estimate income elasticities, using the Industrial Production Index and the Electric 
Power Consumption, for the 1980.01-1988.12 period, and monthly data have been 
used in the estimations. In these works, the series is also modelled by the Box- 
Jenkins methodology. 

An original model is offered in this paper in order to estimate price income and 
cross elasticities for the 1 980.1- 1992.IV period, applying cointegration techniques 
and using monthly data. Such techniques allow the estimation of short run 
elasticities, which add immediate responses to price and income changing, and the 
estimation of long run elasticities which allow to see the effects of the price and/or 
income changing produced later on. 

This research offers a model according to the second proposal above, being 
based on a microeconomic analysis, which can be considered as classic. Its 
structure is very simple. Point number 2.2 presents the theoretical model for 
Spanish passenger transport demand. In the next point, the data used are 
described. Point number 2.4 presents the estimations of the different demands. 
Finally, the main conclusions are offered in point number 2.5. 
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2.2 

The Model 

Suppose a typical user whose preferences are weakly separable, then modelling 
passenger transport demand constitutes the second stage of a two-stage budget 
process. Therefore, the user's spending firstly falls into two large categories: 
passenger transport Services and all other goods and Services; secondly, the user's 
income is assigned to the goods and Services contained in each of these two 
categories. That is to say, the utility function of the representative user is as 
follows: 

U = U (X 1? X 2 ,..., X k ; X k+1 ,..., X n ) 

where the vector = (X 1? X 2 ,..., X k ) with i = 1, 2,..., k represents passenger 
transport Services; the vector Xj = (X k+1 ,..., X n ) with j = k+l,...n represents goods 
and Services except for those corresponding to passenger transport, and U 
represents a utility function which is continuous and differentiable, monotonic, 
increasing, and strictly quasi-concave. 

The consumer balance is reduced to: 

max U (Xi; Xj) 

subject to: Pi • Xi + P, • Xj = Y 

where Pj = (P b P 2 ,..., Pk) and Pj = (P k+ i,..., P n ) are the prices, and where Y 
represents the user's level of income. 

First order conditions allow obtaining the following typical user's Marshallian 
demands: 

Xj = Xi (P i? Pj, Y) (2.1) 

Xj = Xj (P i? Pj, Y) (2.2) 

Of these individual demand functions, function (2.1) is interesting for us since it 
corresponds with passenger transport Services. 

Equation (2.1) still presents some problems. First, functions such as (2.1) 
should be valid for any income distribution among the different economic agents. 
If this were not the case, this function would provide as many values as income Y 
distributions among the possible users and, therefore, such a function would not 
exist. Another assumption would be that income is distributed under a specific 
rule. Once this rule has been established, integrability conditions are verified and 
the existence of the aggregate Marshallian demand functions is guaranteed 
(Varian, 1992). However, there are no data to go along these lines. In order to 
solve this problem in this study we can assume that all the users have the same 
level of income. 

Function (2.1) is general enough to analyse the passenger transport Service 
demands -long-distance railway, air and road transport- identifying the different 
subindexes for the amounts demanded in each Service. 

From 1980.1 to 1992. IV, passenger transport Services in Spain have been 
provided under different regulation conditions. The govemment Company RENFE 
and Iberia have the monopoly of railway and air national transport in Spain 
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respectively, and road transport is provided by private Companies which have 
exclusive routes under a System called "right of testing". It can be said that trump 
road passenger transport, which has a low incidence in quantitative terms, is the 
only mode of transport, which has not yet been regulated. However, given the 
impossibility to obtain quarterly Statistical data on passenger road transport, and 
with the aim of adding interregional transport on the user’s own vehicles, the 
premium petrol consumption variable has been used. The premium petrol 
consumption has also been regulated by the govemment during the period of this 
study. Under such regulatory conditions and with the aim of avoiding any 
Problems arising from the supply-demand simultaneity, we assume that supply is 
exogenous in relation with prices and income and is determined by the decisions 
of the govemment. 

2.3 

The Data 

Data on the series of passenger departures and arrivals of national fligths in 
Spanish airports (AERV) have been obtained from the series provided by the 
Infomies de Coyuntura del Ministerio de Transportes, Turismo y Comunicaciones 
(Reports on the State of the Transport, Tourism and Communication Ministry) No 
data on road transport passengers are available and a "proxy" such as premium 
petrol consumption has been used in order to approximate the transport on the 
user’s own vehicles. This variable, QGAS, has been obtained from the Direcciön 
General de Previsiön y Coyuntura del Ministerio de Economia y Hacienda 
(General Management of the Finance Ministry Forecast and State). The gas-oil 
consumption variable (QGLEO) has also been used with the aim of approximating 
the behaviour of regulär and tramp passenger transport in public Services. 
However, the results obtained are significantly anomalous and the reason for this 
may be that this variable shows the behaviour of road transport of goods (which is 
much more important in terms of consumption), rather than of passengers. 

The data on the series of long distance railway prices (PF) and air transport 
fares (PA) have been obtained from the monthly series worked out from the fares 
of the Boletines Oficiales del Estado (Official State Reports), evaluated within the 
period in which each tariff is in force. 

The data on the prices of premium petrol (PGAS) have been obtained from the 
Direcciön General de Previsiön y Coyuntura del Ministerio de Economia y 
Hacienda, as monthly data, also evaluated within the period in which each fare is 
in force. 

Data on the prices of gas oil (PGLEO) have been obtained from CAMPSA until 
1992. From then onwards, the data from the Compama Logistica de Hidrocarburos 
(Hydrocarbon Logistic Company) have been recorded for further studies. 

The data on the income variable have been obtained considering the Spanish 
quarterly GDP as "proxy”. The series used for the 1980.1- 1989. IV period come 
from Contabilidad Nacional Espanola (Spanish National Accounting). 
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2.4 

Marshallian or Non Compensated Demands of Interurban 
Passenger Transport: Air and Road Transport 

In this paper we use cointegration analysis and error-correction modelling 
techniques to estimate transport demand fimctions. The analysis of cointegration 
offers a generic route to test the existence of long-term relationships, suggested by 
the economic theory, among various integrated variables, which has provided the 
most satisfactory results of the various approaches previously attempted (Inglada 
(1992), Coto-Millän and Sarabia (1994)). Engle and Granger (1987) formalize the 
equivalence between cointegration and error-correction. 

Before an error-correction model can be constructed according to Engle and 
Granger’ s representation theorem, we need to establish that variables in demand 
equations are integrated, at most integrated of order 1, 1(1), and are cointegrated. 
First of all, we examine the properties of individual time series by means of 
Standard tests. The results of these tests suggest the presence of a single unit root 
in all series and this implies that the levels of these variables are nonstationary but 
the first difference is stationary. 

The existence of a long run relationship among variables has been verified in 
this paper by applying the Engle and Granger (1987) technique and the Johansen 
(1988) and the Johansen and Juselius (1990) procedure. The Engle and Granger 
approach to test the existence of cointegration is characterized by a Dickey-Fuller 
statistic, which is used to test the existence of a unit root in the residual of static 
cointegration regressions. Johansen (1988) and Johansen and Juselius (1990) 
develop a maximum likelihood estimation procedure that has several advantages 
over the procedure suggested by Engle and Granger, because it eliminates the 
assumption that the cointegration vector is unique and it takes into account the 
error structure of the underlying process. 

We have estimated some equations from the specifications in model (2.1) by 
adjusting the variables to each mode of transport. All variables headed by letter L 
are in natural logs and those headed by letter D are in differencess. The Statistical 
“t” is presented within brackets under each coefficient. 

2.4.1 Air Transport Demand 

Long run 

The estimated equation of long run balance cointegration has provided the 
following results: 

LAERV t = -1.88 -1.38 LPA t + 1.48 LGDP t 
(-2.13) (-6.43) (21.14) 

R 2 adjusted = 0.91 ; S.E. = 0.04 ; D.W. - 1.25 ; D.F. a = -4.14 ; D.W. b = 1.79 
a Indicates Statistical significance at the 5% level. 

b Is the Durbin- Watson from equation used to computed the DF (Dickey-Fuller) 
statistic. 
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Table 2.1 presents the results of cointegration tests applying the Johansen 
methodology. To test the null hypothesis of at most r cointegrating vectors, 
Johansen (1988) proposed the trace test . The critical value for these statistics has 
been generated by Simulation by Osterwald-Lenum (1992) and the Pantula 
principle has been employed to determine the cointegration rank. Then, if 
Johansen technique is applied to a VAR along with three lags and a restricted 
constant, it is concluded that there is only one cointegration vector. 



Table 2.1 



Number of cointegration vectors 


Under H 0 


Under Hi 


Trace test 


95% critical value (a) 


r = 0 


r > 1 


41.14 


34.91 


r < 1 


r > 2 


19.60 


19.96 


r < 2 


r =3 


6.83 


9.24 



(a) Critical values are from Osterwald-Lenum (1992) 



After normalisation, the following long run solution is obtained: 

LAERV t = -1.76 -1.41 LPA t + 1.17 LGDP t 

In both estimated equations, the long run elasticity of air transport demand with 
respect to the GDP is close, somewhat higher than the unit ad takes 1.17 and 1.48 
values as it would correspond to normal goods and particularly to "luxury" goods. 
The estimated long run own-price elasticity of goods is negative with values 
ranging from 1.38 to 1.41, which reflects a significant response of the demand to 
price changing. 

Short run 

All the demands of passenger transport have been specified, according to the 
Granger representation theorem, in the form of a model with error correction 
mechanism (ECM). This model incorporates the long run relationships, contained 
in the ECM, as well as the dynamics implied by the deviations from this 
equilibrium path and the adjustment process to recover it. The coefficients of the 
variables in differences represent short run elasticities. The joint non-linear 
estimation presents the following results: 

DLAERV t = -0.51 (LAERVt.! + 2.24 +1.12 LPA M - 1.47 LGDP M ) 

(-3.92) (1.96) (2.06) (-9.8) 

- 0.43 DLPAERV t _ 4 + 0.45 DLPGAS t - 0.78 DLPA t 
(-3.07) (2.14) (-2.36) 

R 2 adjusted = 0.95 ; S.E. = 0.048; F = 147.13; D.W. = 2.10 
Serial Correlation: Ljung-Box: Q(l) = 0.30; 
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Q(2) = 0.19; 

Q(3) = 1.36; 

Q(4) = 2.44 

Residual Normality: Bera-Jarque: N(2) = 1.08 

Heteroscedasticity: ARCH (1-4) = 1.27 

The long run elasticities obtained for this and the previous model do not differ 
from each other significantly. Then, long run income elasticity is now 1.47 in 
comparison with the former values 1.17 and 1.48, as it corresponds to luxury 
goods or Services. The negative value of the own-price elasticity of goods is 1.12 
in comparison with the former 1.38 and 1.41 values. 

Short and long run elasticities are once more slightly different. Short run 
elasticities clearly present the inelastic feature of the demand, and a Substitution 
effect of road transport, which has never been revealed before, is detected. Gross 
and net Substitution relationships between air and road transport result once more 
from these estimations. 

2.4.2 Road Transport Demand 

Long run 

In the inter-city road passenger transport demand equation, the dependent variable 
is the amount of premium petrol, in logs, LQGAS, and the independent variables 
are log real rpices of premium petrol, LPGAS and LGDP already defined. The 
equation of the long-term balance cointegration estimated, yielded the following 
results: 

LQGAS t = -3.80 - 0.13 LPGAS t + 1.11 LGDP t 
(-3.21) (-1.94) (8.29) 

R 2 adjusted = 0.94 ; S.E. = 0.03 ; D.W. = 1.51; D.F. a = -5.52 ;D.W. b = 2.01 
a Indicates Statistical significance at the 5% level. 

b Is the Durbin- Watson from equation used to computed the DF (Dickey-Fuller) 
statistic. 

Table 2.2 presents the results of Johansen's cointegration test. This test strongly 
reject the null hypothesis of no cointegration (r = 0), but not the null hypothesis of 
at most one cointegrating vector, so there appears to be a single cointegrating 
vector, wich implies, after normalization, the following long run solution: 

LQGASt - 2.85 - 0.47 LPGAS t + 0.3611 LGDP t 
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Table 2.2 



Number of cointegration vectors 


Under Ho 


Under Hi 


Trace test 


95 % critical value (a) 


r = 0 


r > 1 


41.87 


53.12 


r < 1 


r > 2 


19.19 


34.91 


r < 2 


r =3 


7.90 


19.96 



(a) Critical values are from Osterwald-Lenum (1992) 



The results obtained from the long run estimations provide elasticities of 0.361 
and 1.11 with respect to the GDP, relationships that characterise these Services as 
basic goods rather than as luxury goods, always within the context of normal 
goods. The own-price elasticities of the goods take the negative values 0.13 and 
0.47, once more referring to essential goods with inelastic demand and slight 
demand variations as a response to tariff changes (if we consider such changes as 
proportional to premium petrol price changing).The gas-oil demand equation 
QGLEO presents very similar values with respect to its price and to the GDP 
variable. 

Short run 

Finally, we estimate an error-correction model to integrate short run dynamics 
with long run equilibrium, which presents the following results: 

DLQGAS t = - 0.69 (LQGAS m +3.88 + 0.15LPGAS t _! - 1.1 1 DLGDP M ) - 



(-4.85) 


(2.20) (1.68) (-6.19) 


-0.36 DLPGASm 


+ 0.34 DLPA, 


(-2.73) 


(2.23) 


R 2 adjusted = 0.95 ; S.E. = 0.036 


' ;F = 212.45 ; D.W. = 2.13 


Serial Correlation: Ljung-Box: 


0(1) = 0.28; 




Q(2) = 1.91; 




Q(3) = 4.81; 




Q(4) = 4.82 



Residual Normality: Bera-Jarque: N(2) = 4.16 

Heteroscedasticity: ARCH (1-4) = 1.17 

The long run elasticities obtained for the ECM equation do not differ from the 
previous model. The value of the GDP long run demand elasticity is now 1.11, 
equal to the Engle and Granger approach, and the negative value of the long run 
own-price elasticity of goods is 0.15, while the former values were 0.13 and 0.47. 

The estimated short run own-price elasticities of goods have the negative value 
0.36 and a cross elasticity of 0.34 with respect to air transport price. In the short 
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run, it is possible to speak about gross Substitution relationships between road and 
air transport. However, it is not possible to meet any conclusion with respect to the 
net Substitution or complementary relationships of these transport Services without 
any further assumption. 

2.5 

Conclusions 

A theoretical model of air passenger transport demand has been presented in this 
paper. With quarterly aggregated Spanish data, equations of inter-city passenger 
air and road transport demand have been specified for 1980.1 and 1992.IV. 

Moreover, different demand function estimations have been carried out using 
cointegration techniques, and have been subject to a wide evaluation which allows 
us to check the adequacy of this method with respect to others used in earlier 
works by Inglada (1992) and Coto-Millän and Sarabia (1994). 

Each specific demand may require more detailed studies, especially road 
transport. However, having carried out the estimations, it is possible to draw 
conclusions as regards income, the own-price elasticity of goods and cross price 
elasticities: 

- Long-term income elasticities are all positive and all the Services are normal 
goods. Income elasticities are very close to the unit for air transport, and 
slightly below the unit for road transport. 

- The own-price elasticities of goods increase parallel to the quality of the 
Service, since they increase with fares, and present values close to the unit for 
air transport. They are clearly inelastic for road transport. 

- All cross elasticities present positive values and they are below the unit. Gross 
and net long run Substitution relationships between air and road transport and 
gross Substitution relationships between road and air transport can be 
guaranteed, but net Substitution relationships between these cannot. 
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3.1 

AVE: Characterisation 

The high-speed train (AVE) is mainly characterised by a new infrastructure, which 
uses the adequate mobile material to be able to obtain highly operating speeds. 
Moreover, its way of management is generally different from that typical of the 
conventional train. In line with this, Plassard (1992) proposes the following defining 
characteristics of this mode of transport: 

- A high speed of 250-300 km/h. 

- A high ffequency as the traffic requires it. 

- A weak capacity (lower than 400 people per train). 

- Slightly higher fares than those of the conventional train. 

With these Operation conditions, he States that the AVE must be the link between 
two big cities, excluding every intermediate Service, since this is the only way to 
combine the two requirements of high speed and satisfactory rates of utilisation. 

Bonnafous (1987) agrees with this assumption and points out that the French TGV 
resembles more a plane than a conventional train if we consider some factors such as 
the distance covered, the capacity and the speed. He finally concludes that its 
structural effects essentially affect the urban poles with the highest population, as the 
air transport does. 

In relation with the Madrid-Seville corridor, which is the object of our research, 
we carry out a more detailed analysis, which enables us, to distinguish between 
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the following three clearly differentiated sub-products within the so-called high 
speed trains: the Shuttle, the Long Distance and the Talgo . Chronologically, the 
first train used was the Long Distance. The Madrid-Seville route started to be 
operated on the 21 st of April 1992 with intermediate stops in Ciudad Real, 
Puertollano and Cordoba. The demand for this sub-product had a great boom, 
especially the Madrid-Seville route due to the Universal Exhibition (EXPO) held 
at Seville in 1992. In October 1993, immediately after this event concluded, there 
was a fare reduction , which allowed this market to consolidate itselr , thus 
compensating the decrease in the demand produced after the closure of EXPO. 
This also led to a rising tendency in the dynamics of the product to an aspect, 
which is always inherent to the stage of maturing of any new product. 

Later, on the 18* of October 1992, a new train called “Shuttle” 4 was 
introduced. This train was characterised by the dramatic decrease in fares at rush 
hours, by the introduction of discounts and the addition of units with a greater 
number of tourist-class seats. The decrease in the average revenues of 
approximately 50% brought about an increase in demand of 275% and 115% in 
October and November of 1993 respectively. Therefore, this new offer, which had 
the initial aim of efficiently using the already existing units, was firmly 
Consolidated. 

Finally, the “Talgo” train started to operate in August 1992 from the 
introduction of the Cordoba gauge interchange Service on the Madrid-Malaga 
route. This line was extended in the summer of 1993 by the introduction of the 
Majarabique (Seville) gauge interchange Service, which offered the possibility of 
travelling from Madrid to Cadiz and Huelva on the new line without changing 
trains. 

Table 3.1 presents the most significant features of each type of supply, which 
determine their differentiation from the High-Speed train. The primary 
differentiating factors between Shuttle and Long Distance are those of demand 
and fares. However, the Talgo differs essentially from the others for having to 
employ varied mobile material because it needs a track with different rail widths. 
This fact, along with the gauge interchange Operation, leads to a lower average 
speed in this train and therefore, to a lower generalised cost reduction than that of 
the rest of the High Speed transport. 

A more detailed analysis shows that the shuttle-service, which has low fares 
and in which discount tickets prevail, has generated a demand essentially for 
commuter joumeys of day retums. This fact is mainly due to the reduction of the 
generalised cost of transport -especially in time and fares -, which also leads to a 
reduction of the price of a supplementary commodity such as the housing. 



The new High Speed line has been built with UIC rail width, which is different from the 
Spanish rail width. The great advantage of the Talgo, is that, due to its variable axle System, 
both types of rail can be used with the minimum waste of time when one infrastructure is 
2 changed into the other. 

^ The average revenue feil from 19 ptas. to 14 ptas. per passenger-km. 

4 This is shown by the level of the load factor which exceeds 80%. 

5 In the Madrid-Ciudad Real, Ciudad Real-Puertollano and Madrid-Puertollano lines. 

For instance, the 171 -km Madrid-City takes slightly longer than 50 minutes, virtually the 
same time as any city center-outskirts joumey takes in Madrid. Moreover, we must take into 
account that the average price of housing in Ciudad Real is almost half that of Madrid. 
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Table 3.1. Discriminating factors of the different sub-products 





SHUTTLE 


LONG DISTANCE 


TALGO 


Routes 


Madrid-C.Real 

C.Real-Puertollano 

Madrid-Puertollano 


Madrid-Seville 
Madrid-Cordoba 
Other routes 6 


Madrid-Malaga 

Madrid-Cadiz 

Madrid-Huelva 

Others 7 


Material 


Gec-Alsthom 


Gec-Alsthom 


Talgo 200 


Infrastructure 


New High Speed 
Track 


New High Speed 
Track 


New and 

Conventional Track 


Fares 
(Average 
revenue/ 
Passenger-Km) 
(PTAs. 1993) 


10.5 


15.4 


10.2 


Occupancy rate 


0.65 


0.84 


0.68 


Type of demand 


Local train with a high 
percentage or commuters 
joumeys 


Long distance 


Long distance 



Source: Own elaboration 



3.2 

Qualitative Analysis 

Table 2 records the results obtained from the successive surveys to the high-speed 
train passengers, both in the Long Distance and in Shuttle trains. These results 
refer to qualitative characteristics. The following conclusions can be reached from 
these data. 

The Long Distance Train 

- Education and Professional level: The Long Distance passenger has higher 
education and Professional level: almost 63% have University education. Of 
these, 22% hold University diplomas while 41% hold higher University 
degrees. As regards professions, 40% are executives and businessmen and 26% 
are technicians. 

- Reasons for travelling: 58% travel for Professional reasons, this is an especially 
significant motivation among the business dass passengers (75%). As regards 
other reasons for travelling, 22% travel for tourism and 20% for family reasons. 

- Fares assessment: A survey carried out before the fares were increased in 
September, 1993 reveals that 64% passengers think that the AVE fares are 



This includes the Cordoba-Seville, Ciudad Real-Cordoba, Ciudad Real-Seville, Puertollano- 
Cordoba, Puertollano-Seville lines. 

This would include lines such as Barcelona-Malaga. 
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adequate. This fact led the Company to raise the fares for the possibility, which 
such a measure offered of increasing retums. 

- Level of and cause for satisfaction: There is an almost complete level of 
satisfaction (96%) and the main causes for satisfaction are: speed (29%), 
comfort (26%), good-quality Service (11%), punctuality (8%), safety (4%) and 
good price (3%). 



Table 3.2. Qualitative characterisation of the demand for high speed Services 



ATTRIBUTES 


PERCENTAGES 




Long Distance Train 


Shuttle Train 


Level of education 


Primary School 


15 


31 


Secondary School 


22 


27 


University Diplomas 


22 


16 


Higher University Degrees 


41 


26 


Professional occupation 


Executives 


17 


6 


Businessmen 


23 


15 


Technicians 


26 


23 


Miscellaneous 


14 


19 


Unemployed 


19 


37 


Reasons for travelling 


Professional 


58 


37 


Tourism 


22 


22 


Family 


40 


41 


Fares assessment 


Very high 


5 


6 


High 


30 


25 


Good 


64 


66 


Low 


1 


2 


Very low 


0 


1 


Level of satisfaction 


Very satisfied 


47 


55 


Satisfied 


49 


39 


Indifferent 


3 


3 


Unsatisfied 


0 


1 


Very unsatisfied 


1 


2 



Source: Own elaboration based on the RENFE surveys (1993) 



The Shuttle Train 

- Education and Professional level: The passenger of this Service has lower 
education in general, 42% of the passengers have University education. Of 
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these 16% hold University diplomas and 26% hold higher University degrees. 
As regards professions, 21% are executives, 21% are businessmen, and 23% 
are technicians, which constitutes lower percentages than in the Long Distance 
train. 

- Reasons for travelling: For those passengers using this Service, business (37%) 
has stopped being the main reason for travelling. However, it is worth 
remarking that 41% of passengers travel for family reasons, which is much 
higher percentage than that for the Long Distance train. 

- Fares assessment: As in the Long Distance train, the survey carried out before 
fares were raised in September, 1993, shows that an elevated percentage of 
passengers (69%) think that the AVE fares are appropriate or low. This 
percentage is even higher than that obtained in the Long Distance Service, due 
to the multiple discounts made in this Service. Therefore, the rise margin of 
fares also seems higher for the Shuttle train. 

- Satisfaction causes and level: Almost all the passengers (94%) show a high 
level of satisfaction by their answers. The main causes for it are the following: 
speed (38%), comfort (30%), punctuality (9%), safety (4%), good-quality 
Service (3%) and good fares (2%). It is worth remarking in this Service that the 
highest values are for speed and punctuality for the great importance given to 
time, in connection with the clear commuter typology, which characterises this 
demand. 

3.3 

The Concept of Generalised Cost 

As far as transport is concemed, the traditional concept of a product price must be 
replaced by a broader concept called the generalised cost of the respective 
transport mode. This cost includes other components apart ffom the purely 
monetary ones. Among these components, which correspond to the various 
attributes which characterise the transport product, the outstanding features, 
according to Garcia-Alvarez (1998), are the frequency, time, punctuality, safety, 
comfort, reliability, intermodality and price. The importance of each one of these 
components when the user chooses a particular transport depends on the type of 
product. 

If we analyse these components in detail, we can define the frequency in 
railway transport as the number of trains, which travel between two points at a 
particular period of time. The demand elasticity with respect to the frequency is 
positive, a fact due not only to the existence of more alternative schedules when 
frequency increases but also to a lower waiting time. The value of this elasticity 
naturally depends on extent of frequency since the more elevated this magnitude is 
the lower the elasticity will be . In addition, one of the most important attributes 



Garcia-Alvarez et al. (1998) study the modal distribution and consider a coefficient due to the 
frequency of each operator measured in the number of Services per day. This modal 
distribution coefficient increases when f increases, which reflects that the market share 
increases when the frequency increases. However, for high values of f, the additional growth 
of this coefficient is increasingly lower (for frequencies of more than 25 Services per day the 
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of transport Service is the travel time, especially in those joumeys, which are not 
for leisure. Obviously, the time value varies depending on the individual and the 
reason for travelling. The total time includes not only the transport time but also 
the time of access and stay at departure or arrival terminals. 

If we analyse fares or the monetary component within the framework of the 
generalised cost of transport, we observe a high demand elasticity with respect to 
price due to different reasons. Firstly, as in the case we are studying, we cannot 
consider the transport product as a first-necessity Service, secondly, there are 
various substitutive products on the market we are analysing and finally, a large 
percentage of the population are car owners, a fact leading to the existence of a 
reference value of the marginal cost of car use, which enables a comparison of 
fares of the different modes of transport. 

On the whole, the safety variable of a mode of transport does not significantly 
and directly conditions the decision of travelling or the choice of the mode of 
transport. In line with this, we only have to observe that the private car is the mode 
of transport with the lowest safety rate but it holds the highest passenger traffic 
rate all over the world. However, for certain individuals, this factor may be 
extremely important when choosing the mode of transport. Their decision is made 
from on a subjective point of view of the corresponding mode of transport rather 
than from an objective consideration. Thus, for example, in spite of the fact that 
the plane is much safer than the car, there are a great number of individuals who 
prefer any other means of transport to the air transport. 

Along with the income level, comfort is increasingly important when choosing 
the mode of transport. This is a complex concept which includes, of course, a 
comfortable atmosphere, an adequate accommodation (dimensions and comfort of 
the seat) and the possibility of the minimum interruption in the passenger’ s day to 
day life during the joumey. In connection with this, the availability of toilets, 
coffee and restaurant Services, television or video significantly contributes to 
increase comfort. 

The data appearing on table 3.3 refers to the reasons for the choice of the AVE 
in comparison with other Long Distance modes competing with it. We can reach 
several conclusions from these data about the importance of the various 
components of the AVE generalised cost. In short, we can state that the main 
causes for the choice of the AVE are, by order of importance, speed, comfort, 
price, novelty and safety. 

A first result, certainly relevant, is the fact that the effect of the comfort 
component is virtually as important (29%) as that of time (30%) and it is even 
higher than the effect of price (11%). This result is especially significant within 
the group of air passengers, where the component of comfort is far more important 
(31%) than that of price (19%), a value which is even equal to the sum of the two 
“traditional” generalised cost components: price and time. 



coefficient increases in a maximum 5%). For f=l frequency, the coefficient value obtained is 
1.187 while for f= infinite frequency, which may correspond for example to the car, the 
coefficient value is 4.55. 
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Table 3.3. Reasons for the choice of the AVE in comparison with other alternative modes of 
transport (% for the Long Distance Train) 





Plane 


Car 


Train 


Bus 


Total 


Speed/Time 


13 


42 


57 


67 


30 


Punctuality 


4 


0 


2 


0 


3 


Comfort 


31 


35 


19 


13 


29 


Price 


19 


6 


2 


2 


11 


Novelty 


11 


3 


9 


3 


9 


Safety/Fear 


6 


10 


0 


0 


5 


City Centre 


4 


0 


0 


0 


2 


Schedules 


6 


0 


8 


10 


5 


Miscellaneous 


6 


4 


3 


5 


6 



Source: Own elaboration based on the RENFE surveys (1993) 



In a more detailed disaggregated analysis of each mode of transport, we can state that, 
the main reasons for the choice of the AVE with respect to the plane are comfort 
(31%), price (19%), speed (13%), novelty (11%) and safety (6%). With respect to the 
car, time is the most important reason (42%), followed by comfort (35%) and safety 
(10%). In relation to the conventional train, the reasons to choose the AVE are speed 
(57%), comfort (19%) and novelty (9%). Finally, with respect to the bus, speed 
(67%) and comfort (13%) are the essential reasons for the choice of the new product. 

3.4 

Comparison Among Different Competing Products 

When we compare the different components of the AVE generalised cost with 
those of competing models, we can observe the great advantages of the former. In 
line with this, we can emphasise the fact that the AVE travel time is much lower 
than that of the remaining modes of transport except for the air transport. This 
time reduction is due not only to its high speed but also to the less time employed 
in accessing to the Service, as well as to a decrease in the distances travelled 
thanks to the new infrastructure. 

Table 3.4 presents the time the main routes take. We must point out that in the 
routes in which the new infrastructure is completely used , the time employed by 
the AVE is less than half that used by the remaining modes of transport, except for 
the air transport. For example, in the Madrid-Seville route, the reduction of the 
travel time in the high-speed train is 2h 45’ with respect to the car, 3h 55’ with 
respect to the bus and 3h 20’ with respect to the conventional train. 



In other lines such as Madrid-Cadiz, Madrid-Huelva and Madrid-Malaga the Talgo is 
employed with High-Speed and conventional track. 
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Table 3.4. Travel time for route 



ROUTES 






TRAVEL TIME 






Air 

Transport 


Car 10 


Bus 


Conventional 

Train 


AVE 


Madrid-Seville 


50’ 


5h 20’ 


6h 30’ 


5h 55’ 


2h 35’ 


Madrid-Malaga 


55’ 


5h 40’ 


7h 15’ 


6h 50’ 


4h 45’ 


Madrid-Huelva 


- 


6h 20’ 


7h 50’ 


7h 40’ 


4h 40’ 


Madrid-Cadiz 


55’ 


6h 30’ 


8h 10’ 


7h 45’ 


5h 


Madrid-C. Real 


- 


2h 5’ 


2h 30’ 


lh 55’ 


55’ 


Madrid-Puertollano 


- 


2h 35’ 


3h 


2h 45’ 


lh 15’ 


Madrid-Cordoba 


- 


3h 55’ 


5h 


4h 25’ 


lh 45’ 



In addition, observing table 3.5, which shows the access time to the Services of 
different modes of transport we see that the time used in the AVE is significantly 
lower than in that used in the conventional train and plane but it is equal to that 
used in the bus. The time of access used in the car, obviously null, is the only one, 
which is lower than that in the AVE. This is precisely one of the main advantages 
of this mode of transport since the user considers the access time more “lost time” 
than the travel time. Moreover, we must take into account that this concept 
includes delay penalties, safety, boarding time and access time to and from 
Services. The advantageous aspects of the high-speed train vary depending on the 
mode of transport with which it is compared. For example, with respect to air 
travel, the new product has great advantages in relation to access - its stations are 
located in the city centre - and to delay since the AVE has a high punctuality rate. 
The advantage of this high punctuality rate is also apparent in the comparison of 
the AVE with conventional trains. 

The use of a new inffastructure also adds a comparative advantage to the new 
product, which is the fact that the different distances travelled are quite lower as 
observed in table 3.6. 

Finally, when comparing the prices of the different products appearing in table 
3.7, we must remark that, on the whole, the AVE is - after air travel - the most 
expensive mode of transport. However, in some lines such as the short-distance 
ones, included in the segment of the Shuttle Train, differences with respect to the 
car and the conventional train are minimal. In every case the bus is the mode of 
transport with the lowest price. 



In order to obtain the V speed, the V = 48 + 72 (l-i/c)' /2 equation has been used for highways, 
while in the conventional road the V = 90 - 22 (i/c) alternative expression has been used, 
where i and c are the intensity and hour capacity respectively. Both expressions come from 
the “Highway Capacity Manual” from the Transportation Research Board (1985). 
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Table 3.5. Access time to Service 11 


Plane 


Car 


Bus 


Conventional 

Train 


AVE 


lh 55’ 


_ 


lh 


lh 15’ 


lh 



Table 3.6. Distances 




Road 


AVE 


Conventional Train 


Madrid-C. Real 


190 


171 


255 


Madrid-Puertollano 


228 


210 


293 


Madrid-Cordoba 


400 


343 


442 


Madrid-Seville 


538 


471 


565 


Madrid-Malaga 


544 


528 


627 


Madrid-Cadiz 


663 


633 


727 


Madrid-Huelva 


632 


581 


675 



Table 3.7. Prices of different modes of transport (PTS. per passenger, 1993) 



ROUTES 








FARES 










Plane 12 


Car 13 


Bus 


Conventional 

Train 


AVE 










l“cl. 


2 nd cl. 


Pref. 


Tourist 


Madrid-Seville 


12650 


7346 


2210 


7880 


5450 


10800 


7900 


Madrid-Malaga 


13700 


7655 


2680 


8750 


6050 


9200 


6600 


Madrid-Huelva 


- 


8630 


2630 


9450 


6520 


9700 


6900 


Madrid-Cadiz 


13950 14 


9181 


2600 


10100 


7250 


10200 


7300 


Madrid-C. Real 


- 


2539 


1440 


3690 


1425 


2500 


2000 


Madrid-Puertollano 


- 


3080 


1780 


4280 


1635 


2900 


2300 


Madrid-Cordoba 


- 


5424 


1650 


6130 


4480 


7900 


5800 



For the Madrid-Seville route. For other routes with departure or arrival in Madrid or Seville, 
the access time is reduced by 10 minutes for the bus, the conventional train and AVE. 
Moreover, for the remaining lines, the access time is reduced by 20 minutes. In planes, this 
concept includes lh 15‘ for access, 30‘ for safety and boarding and 10‘ for delay penalty. In 
the bus and the AVE, this access times includes 15‘ for safety and boarding and 45 ‘ for the 
transport to and from the Station. Finally, a 15‘ delay penalty has been considered for the 
conventional train respect to the AVE. 

Price for tourist dass without weekend discounts. 

The occupation for vehicle is of 1.8 individuals and half the depretiation is due to the number 
of Kilometers travelled. 

Madrid- Jerez route. 
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3.5 

Induction and Substitution Effects 

The introduction of the AVE produces two clearly differentiated effects on 
transport demand viz. inducing and Substitution effects. These effects correspond 
to those joumeys which would have never been made if this new Service had not 
existed and to those which would have been made but using a different mode of 
transport. 

3.5.1 Induction Effect 

Within the inducing or generation effect of new joumeys, we must include not 
only those passengers who had never travelled on that route but also the 
component, which accounts for the increase in the travel frequency of those 
passengers who already used that route before the existence of the AVE. In table 
3.8, we observe that the average joumeys per year on the Madrid-Seville route 
have increased ffom 11 to 15 due to the introduction of the new product and that 
28% of the users had never travelled on that route before. The inducting effect still 
continued, but at an increasingly lower rate, until 1996, which marked the end of 
the final period of maturing of the Service. The results obtained show the 
importance of this effect, which, according to the studies made, once the jperiod of 
maturing was completed, accounts for almost 45% of the AVE passengers . 



Table 3.8. Travel frequency in the Madrid-Seville route 





Before the AVE 


After the AVE 


Twice of more a week 


3% 


6% 


Once a week 


6% 


8% 


Once a fortnight 


7% 


9% 


Once a month 


16% 


16% 


Once a quarter 


18% 


16% 


Less 


30% 


17% 


Did not travel 


21% 


- 


Travelled on this route for the 


- 


28% 


first time 

Average joumeys per year 


11.1 


15.2 



Source: Own elaboration based on the RENFE surveys (1993) 



15 



A result similar to that indicated by Nash (1991) for the French Paris-Lyon line. 
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3.5.2 Substitution Effect 

The reduction of the generalised cost in the new mode of transport in comparison 
with the substitutive modes generates the change of part of the demand for these 
competitive modes of transport to the High Speed Train. Obviously, the intensity 
of this “Substitution effect” is lower in the Talgo train because the reduction of the 
generalised cost is also lower. 

The comparative analysis made in point 4 about the different components of 
transport generalised cost indicates that the high-speed train offers certain 
advantages in comparison with the remaining modes of transport. In line with this, 
as inferred from the reasons which lead to the choice of the AVE detailed in table 
3, this new product is preferred to the air transport due to its lower price (19%) 
and higher comfort (31%) and safety (6%). Even a large number of passengers 
(13%) prefer this mode of transport for its reduced access time and waiting time, 
since these two aspects are considered a great time loss. 

The large number of advantages presented by the AVE with respect to the 
conventional train, especially as regards time and comfort, causes such a high 
magnitude of the Substitution effect that the demand for this mode of transport is 
virtually absorbed by the AVE. With respect to the car, the main causes for the 
choice of the AVE are the time (42%), comfort (35%) and safety (10%). Finally, 
in comparison with the bus, the main causes for Substitution are time (67/%) and 
comfort (13%). 

3.5.3 The Demand for the New Product 

The demand for the AVE comes from the combination of the induction and 
Substitution effects the above section. When we analyse the time evolution of 
these two effects, a period of maturing is observed during which the level of 
demand for that product adjusts and grows exponentially until it stabilises. This 
period of maturing is common to all new products. As a reference, we must 
remark that when the High Speed was introduced in the Paris-Lyon corridor, its 
period of maturing lasted 6 years. However, that period was lower - 4 years - in 
the Madrid-Seville case for various reasons . A first reason for this time 
reduction was the demand shock generated by the EXPO, especially in the 
Madrid-Seville line. A second reason was the significant decrease in fares 
occurred in October 1992 at the end of the EXPO. This fall of prices, especially 
marked for the new shuttle-service, accelerated the rhythm of demand absorption 
by the ^new mode of transport. Consequently, the period of maturing was 
reduced . 



This reduction is mainly due to the decrease in fares and the increase in comfort with respect 
to the plane as well as by the reduction of time and the increase in comfort with respect to the 
bus, the conventional train and the car. 
lg This period lasted 6 years in the TALGO train. 

However, the TALGO fares were kept almost fixed so we cannot consider them as a 
reduction of the period of maturing. 
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3.6 

Impact on Demand 

Due to the high magnitude of the Substitution effect, the introduction of the high- 
speed train causes very strong effects on the demand for the remaining modes of 
transport competing with it in the Madrid-Seville route. A more detailed analysis 
enables us to observe the highly strong effect particularly on the conventional train 
since this lost almost 78% of its traffic prior to the introduction of the AVE, which 
virtually produces the disappearance of this transport for the route. The 
introduction of the High-Speed train has led to an approximately 50% fall of the 
air transport demand in the Madrid-Seville route. Therefore, we can remark that 
the operating Company has acted efficiently when it adjusted its supply to the 
demand fall by reducing the number of flights and trying to keep the occupation 
percentage almost stable . As regards the car, the losses are lower than in the 
previously mentioned modes of transport, approximately 30% in the Madrid- 
Seville route. Finally, in bus transport it seems that the impact on long-distance 
routes is not great (11% of demand fall) since both products are hardly 
substitutive. However, in the Madrid-Ciudad Real route, the demand fall, which 
reached 34%, is significant due to the low fares of the high-speed train. 

Therefore, we can conclude that the introduction of the high-speed train causes 
a dramatic change in the modal distribution of the demand, and we can speak of a 
pre AVE transport market and a post AVE market. This effect is especially 
important in the long-distance segment, where the railway is the predominant 
mode of transport in the Madrid-Seville route, and where its use even exceeds the 
use of the private car, a most unusual Situation in the Spanish transport market. 
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4.1 

Introduction and Reasons for Analysis 

The hub-and-spoke Systems are frequently connected with maritime, air and road 
transport. The international Container traffic in maritime transport is developed 
under this model. For a port to achieve a “ Hub Status ” in maritime Container 
transport, its traffic should exceed 1,000,000 tons per year and be provided with the 
necessary equipment to use simultaneously three post-panamax cranes for a new- 
generation vessel. Huge Container vessels dock only at Hub ports and international 
shipping Companies transfer their Containers from these ports to their round the world 
and feeder lines. 



It is known that feeder vessels are used for the transport of Containers from or to minor ports 
which are not served by the huge vessels pertaining to the round the world lines, and make 



P. Coto-Millan (ed.), Essays on Micro economics and Industrial Organisation 
© Physica- Verlag Heidelberg 2004 
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National air networks are also established as hub-and-spoke Systems in connection 
with international transport. Small airports have barely any Connections with 
foreign airports or even with other minor national airports, therefore, indirect 
flights through hub airports network are necessary in order to reach those airports. 
We can also find this System either in road transport (with the Logistics Area 
acting as hub) or in railways, whose stations facilitate the transfer between lines. 

Given that we are dealing with an Organisation System, in the hub-and-spoke 
System - which is constantly present in Transport Economy - it is necessary to find 
methodological formulae which enable us to tackle the study of this System in a 
common manner, notwithstanding the means of transport implied. In this paper, 
we deal with one of those possible common approaches ffom the econometric 
analysis of the time series by means of the SVARs (Structural Vector 
Autoregressions) models with Long-run Restrictions. Other alternatives can be 
found in the BVAR (Bayesian Vector Autoregressions) or in the VEC (Vector 
Error Correction Models) with restrictions. All these approaches allow us to 
introduce the information available a priori on the hub-and-spoke System and 
make it compatible with the Statistical information obtained ffom the evolution of 
traffics in the past. 

4.2 

VAR Models: Methodological Approach 

Düring the last 20 years, we have witnessed important changes in the econometric 
modelling of time series, mainly due to the weaknesses of the economic theory 
when dealing with complex models of simultaneous equations. On rare occasions 
does the theory offer us a dynamic specification which clearly explains the 
relationships between variables, even less frequently in fields of more recent 
development, such as transport economics. This fact has led to the application of 
non-structural models which explain the relationships between the different 
economic variables. Perhaps the models most frequently used have been the VAR. 
These models are employed both to predict and analyse the impact over time of 
some variables in comparison with others, and even to demonstrate hypotheses of 
the variables’ behaviour, although, as regards the latter use, their evidential 
valence is, without doubt, lower than that of the traditional simultaneous equations 
Systems. VAR models also avoid having to formulate structural relationships 
between variables since they state that each endogenous variable of the model is a 
function of their lag values together with a set of exogenous variables. Below, we 
show an autoregressive representation of a VAR model: 



only intermediate stops at the large transoceanic stations. Since feeder vessels operate mainly 
in small ports with less gauge, these usually have a lower capacity, between 500 and 700 tons, 
although some may reach 1,500 tons. In contrast with the huge Container vessels and, in view 
of the limited Provision of mobile facilities in the small ports where they operate, the vessels 
are indeed provided with their own loading and unloading means; they usually have two cranes 
2 located in the middle-line or on the port side of the vessel. 

It is easy to have information available a priori on a hub-and-spoke System, even though it is 
only the structure of the lines and their commercial importance measured taking into account 
the volume of traffic. 
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y, = SiAjy,.; + ZjB j x,. j + e t = A,y t .,+ ... + A p y t . p + B 0 x,+ . . . + B q x,. q + e t . Vi = 

l,...,p Vj = 0,...,q. 

Being: 

y t = endogenous variables vector of k x 1 ränge. 

x t = exogenous variables vector of d x 1 ränge. 

Ai = matrix of the endogenous variables ratios with i-lags and whose ränge is k 
xk. 

Bj = matrix of the exogenous variables ratios with j-lags and whose ränge is d 
x d. 

e t = k x 1 dimension innovations vector. These innovations have zero mean and 
constant variance, and they are generally correlated between one another if they 
are Contemporary, but incorrelated if they correspond to different time periods. 

4.2.1 Impulse Response 

As mentioned above, one of the main uses of the VAR is the analysis of the 
dynamic transmission mechanisms of the effects of changes undergone by the 
variables. A shock or change in the i variable will affect all the variables of the 
model through a structure of lags of the VAR model. A priori, the interpretation of 
the Impulse Response fimctions is not so simple since the innovations observed in 
the model (e t ) will generally present Contemporary correlations between one 
another. In line with this, it can be stated that the innovations have a common 
component which cannot be associated to a specific variable. 

In order to make the interpretation of the Impulse Response fimctions easier, 
the innovations must be incorrelated either over time or during the same time 
period. The incorrelation is obtained by applying a transformation to the 
innovations in the following manner: 

Bu t = Ae t being u t = B _1 Ae t - (0,lu t ) 

3 

Being Zu t a diagonal covariance matrix . The Impulse Response fimctions are 
generally estimated to show the response in the model to a shock of a Standard 
deviation in the structural innovations. Due to this, it is necessary that a SVAR 
model be standardised by adjusting the variances of the structural innovations to 1 , 
since, in this way, a shock of a Standard deviation corresponds to a unit innovation 
in u t . Naturally, the Standardisation only affects the scale of the System formulated 
and does not alter any relationship or substantial feature of this System. 

In order to orthogonalise the impulses, we have used in general the Cholesky 
factor of the residual covariance matrix. This System has important shortcomings 
since the results obtained depend on the ordering of the variables so that, the 
responses to the impulses will be affected according to this order. The Cholesky 
decomposition can only be employed when the order of the variables proposed has 
it foundations in strong theoretical assumptions. This problem has been 
traditionally overcome by ordering the variables ffom high to low exogeneity, for 



3 

Technically speaking P = B' 1 A matrix is standardised in such a way that P 2e t P’ = Eu t = I, I 
being the identity matrix. 
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which we should apply Granger’s causality tests. Another solution is to employ 
the SVAR models, as set out in the section below. 

4.3 

Formulating Long-Run Restrictions under a SVAR Model 
Applied to a Hub-and-Spoke System 

Prior to explaining how to formulate restrictions under the SVAR model, we 
should prove its utility. SVAR models have become an increasingly populär tool 
during the past decade. Among their main uses is the analysis of the economic 
cycle’s sources, and they are also employed in the study of the transmission 
mechanisms of monetary policies. Thanks to short- and long-run restrictions, the 
error term can be orthogonalised in a non-recursive way for the analysis of the 
Impulse Response functions. In other words, it is a an alternative to Cholesky 
recursive orthogonalisation for which it is necessary to have sufficient restrictions 
to identify the orthogonal components of the error term. 

The literature of the SVAR models makes a difference between short- and 
long-run restrictions. In this paper we focus on long-run restrictions mainly 
because we deal with annual series. The way in which short-run restrictions are 
formulated usually implies that one of the variables cannot instantaneously 
respond to the changes or shocks produced in other variables. These changes or 
shocks are unpredictable so, in the analysis of the hub-and-spoke models, we 
would be stating that the managers of one of the transport infrastructures of the 
System - let us call it X -, need some time to evaluate the changes in the traffic 
statistics of the infrastructure which has undergone the shock, called Y. In line 
with this, the short-run restriction would mean that, during the following period, 
the response in the traffic of X infrastructure to a shock in the Y infrastructure of 
the System would be 0. 

Given that transport infrastructures (ports, airports, railway stations...) operate 
with the newest information Systems, these hypotheses can only be kept when we 
work with monthly, term or quarterly series rather than with annual series. This 
fact reduces the scope of application of short-run restrictions to port and railway 
traffic; although it is possible to find monthly statistics in port traffic, these are 
usually mere approximations. In spite of having traffic statistics available every 
six months - or even more regularly - it is difficult to keep short-run restrictions 
useful for the analysis of the hub-and-spoke transport Systems. We cannot 
overlook the fact that a positive (or negative) shock in a spoke infrastructure 
immediately implies an increase (or decrease) in the hub infrastructure since the 
hub infrastructure is the natural entry and departure point of traffic in the spoke 
infrastructure. 

On the contrary, there is a priori a wide scope for the employment of SVAR 
models with long-run restrictions. These restrictions can be stated under quite 
realistic assumptions, as can be seen below. Next, we account for the analytic 
formulation of the long-run restrictions included in a SVAR model , the 



5 For the analysis of short-run restrictions, see Amisano and Giannini (1997). 
For a deeper analysis see Blanchard and Quah (1989). 
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components of the C matrix being the long-run accumulated response to the 
shocks in the structural innovations : 

C= 'Poo A'B 



=(I-Ä 1 -...-Ä p ) 1 

being the estimation of the accumulated response to the shocks in the innovations 
observed. 

The Impulse Response functions for a hub-and-spoke System will allow us to 
check how a positive or negative shock in one of the transport infrastructures will 
affect the remaining ones, mainly in the hub System. Within the ffame of this 
general formulation we can distinguish between two main aims: 

1. It is especially interesting to study how a Variation in the traffic of an 
infrastructure which plays a “spoke” role will affect the “hub” infrastructure of 
the System. With this study, we will be able to contribute to a better distribution 
of the “hub” infrastructure’ s capacity thus avoiding the appearance of 
bottlenecks or reducing the appearance of idle time periods. 

2. In our SVAR model, it may happen that the transport infrastructures considered 
share their hinterland or scope of influence, that is to say, the area where the 
traffic is generated because the industrial Companies which depend on the 
existence of the infrastructures studied are located there. Therefore, we will 
have to take into account that a shock in an infrastructure will be composed by 
two effects with different signs, one positive, which implies the creation of 
commerce, and another one negative which means a trade deviation. If the 
hinterland is shared, the positive effect which the favourable shock of one of 
the spoke infrastructures would have on the hub, is reduced by the traffic loss 
produced by trade deviation. We can even find that a positive shock of a spoke 
infrastructure may have a negative accumulated response on its hub 
infrastructure. This will happen whenever we are dealing with a multiple hubs 
System. In the event of this, the market share of the spoke infrastructure which 
undergoes the shock, held by the hub infrastructure analysed, shall be quite 
small. In other words: 

Chs = k(CMhs*Css - DC) (4.1) 

Being: 

Chs = Accumulated response on the hub infrastructure to a shock in the spoke 
infrastructure. 

CMhs = Market share of the hub infrastructure over the set of traffics of the 
spoke infrastructure. This share is calculated as the traffic rate of the spoke 
infrastructure which has its origin or destiny at the hub infrastructure studied. 
Naturally, this rate shall depend on multiple factors, among which we can stand 
out the following: the number of hubs with which the spoke infrastructure 



Therefore, CI 2 = O is interpreted in such a way that the long-run effect on the variable 1 of a 
shock produced in variable 2 is nil. In other words, variable 2 will only have transitory effects 
on variable 1 . 





46 



J.l. Castillo-Manzano, P. Coto-Millän and L. Lopez-Val puesta 



operates - and the geographic proximity between these and the spoke 
infrastructure or their relative importance within the transport category 
implied. 

Css = Accumulated response on the spoke infrastructure to a shock produced in 
it. 

DChs = This item includes that part of Css traffics which does not imply trade 
creation - in other words, new traffics - but they are the result of a trade 
deviation from the hub infrastructure towards the spoke infrastructure. As 
mentioned above, this effect can only be stated under the initial assumption that 
the hub and spoke infrastructures share a geographic area known as “shared 
hinterland”. 

k = Is a positive proportionality constant. 

In order to implement a SVAR model and be able to measure the objectives 
exposed above (eq. (4.1)), it is necessary to establish a criterion under which long- 
run restrictions are to be established, based upon the economic behaviour of the 
hub-and-spoke Systems. A priori, it is logical to assume that in a bivariate SVAR 
model Csh = 0, in other words, that the long-run effect of a shock in the hub 
variable on the spoke infrastructure is nil. Here we are assuming that, a change in 
the traffics of the hub infrastructure (a shock in the hub traffic) is due either to the 
increase in the demand of the hub’s hinterland or to a former shock in another 
spoke infrastructure in connection with it. A different view of this would be 
provided by deriving the following analytic expression from formula 1 : 

Csh = c(CMsh*Chh - DCsh) = 0 (4.2) 

This restriction is based upon the following initial assumptions: 

CMsh = Market share of the spoke infrastructure over the set of traffics of the 
hub infrastructure. Analogously to CMhs, this share is calculated as the traffic rate 
of the hub infrastructure which has its origin or destiny at the spoke infrastructure 
studied. Also in this case, this rate will depend on a multiplicity of factors such as 
the number of spokes with which the hub infrastructure operates, the geographic 
proximity or the cost structure of both infrastructures, among others. Presumably, 
the more important the hub infrastructure is, the lower this coefficient will be and 
therefore, CMsh*Chh will tend to 0 or be slightly over zero. 

Chh = Accumulated response on the hub infrastructure to a shock in it. 

DCsh = Trade deviation from the spoke infrastructure towards the hub 
infrastructure. Again, this effect shall be estimated if the hub and the spoke 
infrastructures share the hinterland. Obviously, this effect shall be less important 
than the DChs since the majority of the variations in the traffics of a hub 
infrastructure shall be due to changes in the demands of the hub or spokes 
infrastructures’ hinterlands in relation with the hub infrastructure. In general, this 



We cannot miss the fact that the demand for transport Services is a derived demand, as all 
demands for productive factors. Therefore, an increase (decrease) in the traffics of an 
infrastructure is a consequence of a former demand for these traffics generated by an increase 
g (decrease) in the economic activity of a hinterland. 

This importance is measured either with respect to the hub’s traffic volume or to the number of 
different connections established by this hub with the spoke infrastructures and with other 
hubs. 
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effect shall be nil, even when a small hinterland area is shared, or, at best, will be 
slightly negative, being compensated with the CMsh*Chh value. Under these 
assumptions, the CMsh*Chh - DC difference will tend to 0 and, therefore, it can 
be maintained, as a long-run restriction, that Csh = 0. 

c = Is a positive proportionality constant. 

Obviously, in order to support this reasoning we must previously assume that 
the hub infrastructure is not congested, in other words, that the spoke 
infrastructure can employ the hub for its transfers as much as necessary. 
Moreover, in order to analytically sustain this assumption, it is highly 
recommended to check that the Csh = 0 value is within the Csh confidence 

9 

interval by applying any recursive orthogonalisation method and elaborating the 
confidence intervals from the estimation of the accumulated responses’ Standard 
errors . 

Under the above reasoning, in the event that the hub infrastructure were 
congested, it would not be of much sense to analyse the Impulse Response 
functions, and all conclusions would be altered. This analysis would only be 
useful as an attempt to simulate what would happen if the hub infrastructure 
extended its facilities, a fact that would explain the positive shock in its traffics. 
Even in this case, if there were congestion problems in the hub infrastructure, the 
above-mentioned restrictions could not be used, in other words, we could not state 
that Csh = 0. 

Under this simple theoretical outline, assuming that no congestion problems 
exist in the hub infrastructure and applying a methodology of tested solvency to 
the dynamic analysis of time series (the SVAR models), a very interesting set of 
possibilities appears in all fields of transport economy. In order to illustrate the 
usefulness of this methodology, a practical application is set forth in the following 
section. 

4.4 

Practical Application to Container Traffic between the 
Bahia de Algeciras Port and other Ports of the Spanish 
Port System 11 . 

As mentioned at the beginning of this work, this methodology could be applied to 
any transport System developed under a hub-and-spoke structure. We should not 
miss the fact that many of the properties of the VAR models - and of the tests 



For this work we have employed Cholesky factorization by ordering the variables from high 
to low exogeneity after applying the Granger’s test and adjusting the number of degrees of 
freedom for the residual covariance matrix. In line with this, the (i,j) element of the residual 
IQ covariance matrix would be I t ei, t ej, t /(t-p), p being the number of parameters per equation. 

This does not imply any difficulty since Standard econometric packages account for the 
possibility of estimating these Standard errors either analytically (asymptotically) or by means 
of the Montecarlo method. For the applications in section 4 we have chosen the analytical 
11 estimation. 

All the estimations of this section have been carried out through the Eviews 4.1 Statistical 
package. 
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12 

applied to them - are for infinite samples . This makes it necessary to seek 
infrastructures with broad time series to which this methodology can be applied, a 
fact which will in general restrict its scope of action to port and airport Systems. 
For this practical application, we have chosen the Container traffic between the 
Bahia de Algeciras port and other ports of the Spanish port System. 

Within the Spanish port System, the most important hub System in the 
international Container traffic is the Bahia de Algeciras port , together with the 
Valencia and Barcelona ports. Although the difference in the volume of Containers 
moved in these three ports is not so big as to justify the selection of the Bahia de 
Algeciras port as the unquestioned hub of the Spanish port System, this 
difference has to be measured due to the extension of their hinterlands. The Bahia 
de Algeciras port has the minor hinterland, which leads to the fact that, 
historically, an average of 90% of its Container traffic does not depart frorn its 
hinterland. This 90% is dedicated to the international or national transfer by means 
of its round the world and feeder lines, and, therefore, this traffic is only due to its 
hub fimction within the Spanish port System. For all these reasons, we have 
chosen it as a reference of our practical application. 

Since time series are not very long (1975-2001) and, in addition, they are not 
stationary, we have decided to work with bivariate VAR models in Order to avoid 
losing all degrees of freedom. This would happen if we worked with a wide 
number of variables. Our objective is to demonstrate the suitability of the two 
scenarios depicted above, in other words, we will study the interdependency of 
these infrastructures in the System whether there be a shared hinterland or not. We 
have chosen a group of applications in which, a priori, we could certainly know 
which are the predicted results in such a way that any deviation frorn these results 
produced by using SVAR models would be difflcult to justify and therefore, the 
utility of these models within the field of the hub-and-spoke Systems economic 
analysis. 

4.4.1 Scenario 1: Existence of a Shared Hinterland between the 

Infrastructures. Practical Applications: Bivariate SVAR between 
the Ports Bahia de Algeciras and Cadiz 

In order to illustrate the first scenario, we have established a bivariate SVAR 
between the Bahia de Algeciras and Bahia de Cadiz ports. The distance between 
these two ports is less than 100 km; it takes one hour to a truck to travel by road 
between these two ports. In line with this approach, the Bahia de Cadiz port would 
be a spoke port with respect to the Bahia de Algeciras port, which would be a hub 



12 

It is widely known that, in time series, the ordinary least squares (OLS) results are the same 
as in linear economy but for infinite series. Therefore, the estimations by OLS of the 
Parameters of a VAR converge in their distribution towards the real values of the parameters, 
provided that the model’s residues are white noise. In order to prove this assumption, it is 
necessary to apply a multivariate normality test to these values. 

The Bahia de Algeciras port, in its function of hub port in Container traffic, has feeder lines 
within the Spanish port System with the Melilla, Alicante, Valencia, Barcelona, Las Palmas 
de Gran Canaria, Santa Cruz de Tenerife, Ceuta and Tarragona ports. 

In 2001, the Bahia de Algeciras port moved 1,414,390 tons in comparison with the 1,144,775 
tons of Valencia and the 947,767 of Barcelona. 
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port. It is certain that the Bahia de Cadiz port is not a spoke port of Bahia de 
Algeciras since it does not have feeder lines with it. 

This fact does not invalidate the application of the above-mentioned 
methodology, indeed, it has clearly favourable effects. On the one hand, it 
reinforces the theoretical justification of the long-run restriction (Csh = 0) for 
various reasons. Firstly, since there is not a feeder line between the Bahia de 
Algeciras and Bahia de Cadiz ports, the shock of the hub infrastructure would not 
have any favourable effect on the spoke port, in other words, CMsh*Chh = 0, 
since CMsh = 0 by definition. Moreover, we must assume that the effects of trade 
deviation derived from sharing a common hinterland (DCsh) will be very low as 
regards the accumulated response of the spoke infrastructure to the hub shock. 
This can be easily supported if we take into account the low magnitude of traffic 
in the spoke ports with respect to the hub ports; for example, the Bahia de Cadiz 
port would account for less than 3.5% of Container traffic in the Bahia de 
Algeciras port, measured in tons. 

Another possibility of analysis for this particular and exceptional case - where 
no hub-spoke relationship exists between the infrastructures covered in the 
analysis - is to consider that DCsh is different from zero. Therefore, we would be 
assuming that a shock in the hub infrastructure has a slightly negative effect on the 
spoke infrastructure. 

Obviously, it is much more difficult to establish a priori by deduction a long- 
run restriction different from zero, either positive or negative. In this work, we 
have used the same analytical method to choose the value a priori of this new 
negative restriction as the one used to test that Csh = 0. To Start with, we have 
established the confidence interval for Csh by - with a view to orthogonalising the 
errors observed - using the inverse of Cholesky factor, then adjusting the number 
of degrees of ffeedom and ordering the variables from high to low exogeneity 
under Granger’s causality model . The highest exogeneity of the variable as we 
expose it, coincides with the concept of causal priority as Granger understands it. 
Once the interval has been established, a slightly negative value within this 
interval must be chosen. 

The first Step to establish a bivariate SVAR model between the Container 
traffics, expressed in tons, of the Bahia de Algeciras and Bahia de Cadiz ports, is 
to survey the series which are to be dealt with. Graphically, we observe a positive 
relationship between the mean and the typical deviation of the variables expressed 
in levels, especially for the case of the Bahia de Algeciras port (see flgure 4.1). 
This evident hint of heterocedasticity suggests the employment of the logarithmic 
transformation in both series. Moreover, without having to resort to a unit-root 
test, we can conclude the nonstationariety in the mean of the level series, as seen 
in flgure 4.2. In Order to avoid these shortcomings, we will establish the SVAR 
model with the following VlogALG and VlogCAD series, in other words, we will 
explain the evolution of the growth rate of each one of the level series depending 
on the last growth rates. Table 4.1 shows the results of the unit-root tests for these 



We consider a variable y 2t to cause, in the Granger’s meaning, another variable yu if the 
future projection of yi t , with all the information available (Yu = yi t + y2t), and the projection 
with only the information of variable yu, are different. 
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series, in which the nil hypothesis is rejected and therefore, the original series 
(logALG and logCAD) are just 1(1). 




MEAN 



Fig. 4.1 . Graph of the mean-typical deviation of the Bahia de Algeciras port series 




Fig. 4.2. Bahia de Cadiz and Bahia de Algeciras port series (standardised data) 



Table 4.1. Unit-root test for VlogALG and VlogCAD 



Test VlogALG 


t-stat 


MacKinnon 

p-values 


Augmented Dickey-Fuller based SIC 


-4.8840 


0.0007 


Phillips-Perron (Newey-West using Bartlett kemel) 


-9.0163 


0.0000 






MacKinnon 


Test VlogCAD 


t-stat 


p-values 


Augmented Dickey-Fuller based SIC 


-5.6475 


0.0001 


Phillips-Perron (Newey-West using Bartlett kemel) 


-5.5837 


0.0001 
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Once the^likelihood ratio tests have been applied to determine the optimal lag of 
the VAR , we obtain the following structure: 

VlogCAD = 0.2 1 76* VlogCAD(- 1 ) + 0.1 193*VlogALG(-l) + 0.0162. 

VlogALG = - 0.3264* VlogCAD(- 1 ) + 0.1522*VlogALG(-l) + 0.1221 



Table 4.2. Optimal lag in the VAR (VlogALG y VlogCAD) 



Lag 


LogL 


LR 


FPE 


AIC 


SC 


HQ 


0 


27.71421 


NA 


0.000438 


-2.057137 


-1.959627* 


-2.030092 


1 


33.32200 


9.869710* 


0.000386* 


-2.185760* 


-1.893230 


-2.104624* 



* indicates lag order selected by the criterion 

LR: sequential modified LR test statistic (each test at 5% level) 

FPE: Final prediction error 

AIC: Akaike Information criterion 

SC: Schwarz Information criterion 

HQ: Hannan-Quinn information criterion 



One data necessary to specify the SVAR model is the confidence intervals of the 
long-run restriction which we wish to establish (Csh = 0). In this case, the ranges 
for an accumulated response of 10 periods are (0.1311,-0.0464), estimated from 
the Cholesky factor by adjusting the degrees of freedom and ordering them from 
high to low exogeneity. Then, it would be possible that Csh = 0 or slightly 
negative taking into account that only trade deviation effects exist. For example, 
an Option would be to state that Csh = -0,010. Table 4.3 accounts for the results of 
the VAR model eventually selected. 

Below appear the graphs (figures 4.3 and 4.4) of the Impulse Response 
fimctions under two choices Csh = 0 and Csh = -0,010. Csh would be Ccad/alg in 
this model. For the results to be coherent with the model exposed, the graphs must 
show that Ccad/cad and Calg/alg will be certainly positive while Calg/cad must be 
negative. The following graphs present the estimations of the accumulated 
fimctions Ccad/cad and Calg/alg under the two hypotheses considered (Csh = 0 
and Csh = -0,010), which, as can be observed, satisfy the assumptions exposed for 
them. 

After analysing the model’ s residues we have applied the Jarque-Bera 
Standardisation test for multivariate series. With the orthogonalisation restrictions 
proposed Csh = 0 and Csh = -0,010, we have obtained p-values of 0,5028 and 
0,4905 respectively. These values are more than sufficient to keep the residues’ 
white noise hypothesis. 
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Table 4.2 shows the results of the tests in order to determine the optimal lag in the model 
studied. The optimal lag for each election criterion is marked with an asterisk. 
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Accumuiated Response to Structurai One S.D. Innovations 

Accumuiated Response of DLOGALG to DLOGALG Accumuiated Response of DLOGALG to DLOGCAD 




Accumuiated Response of DLOGCAD to DLOGALG Accumuiated Response of DLOCAD to DLOGCAD 




Fig. 4.3. Impulse Response functions accumuiated under the Csh = 0 restriction 





Fig. 4.4. Impulse Response functions accumuiated under the Csh = -0,010 restriction 
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Table 4.3. VAR Bahia de Algeciras and Bahia de Cadiz ports 




VLogALG 


VLogCAD 


VLogALG(-l) 


0.152218 


0.119393 




(0.14036) 


(0.14995) 




[1.08449] 


[ 0.79620] 


VLogCAD(-l) 


-0.326435 


0.217643 




(0.16384) 


(0.17504) 




[-1.99238] 


[1.24339] 


C 


0.122171 


0.016222 




(0.03620) 


(0.03868) 




[ 3.37482] 


[0.41944] 


R-squared 


0.305305 


0.067350 


Adj. R-squared 


0.242151 


-0.017436 


Sum sq. Resids 


0.364097 


0.415571 


S.E. equation 


0.128646 


0.137439 


F-statistic 


4.834294 


0.794353 


Log likelihood 


17.39166 


15.73875 


Akaike AIC 


-1.151333 


-1.019100 


Schwarz SC 


-1.005.068 


-0.872835 


Mean dependent 


0.146757 


0.038944 


S.D. dependent 


0.147777 


0.136257 


Determinant Residual Covariance 




0.000308 


Log Likelihood (d.f. adjusted) 




30.12616 


Akaike Information Criteria 




-1.930093 


Schwarz Criteria 




-1.637563 



4.4.2 Scenario 2: Lack of Shared Hinterland between Infrastructures. 

A Practical Application: Bivariate SVAR between the Bahia de 
Algeciras and Las Palmas Ports 

In this scenario we have two transport infrastructures which are sufficiently far 
apart so that no hinterland is shared by them. This makes the analysis much 
simpler since the effects of trade deviation between the infrastructures shall be nil 
and we therefore have that DCsh = DChs = 0. Within this scenario our model’s 
equations would be as follows: 

Chs = kCMhs*Css. 

The long-run restriction being: 

Csh = cCMsh*Chh = 0, since, as explained above, CMsh will tend to 0 in a truly 
hub infrastructure with multiple national and international connections and much 
more global traffic than the one of the spoke infrastructure. 

The current Spanish ports which have a feeder line in connection with the 
Bahia de Algeciras port are Alicante, Barcelona, Bilbao, Cartagena, Las Palmas de 
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Gran Canaria, Melilla, Santa Cruz de Tenerife, Tarragona, Valencia and Vigo. 
However, not all these ports are suitable for our analysis. Firstly, we must discard 
those ports which have hub infrastructures at the same time, since we cannot 
assume that Csh =0 when both ports are hub. Moreover, since these ports would 
compete for the same transoceanic lines within the same geographic area - in this 
case within the westem Mediterranean - the effects of trade deviation between 
them will be strong and would have to be taken into account, even though the 
hinterland shared between the Bahia de Algeciras port and the rest of hub ports 
may be nil. For ^ these reasons, we must exclude the analysis of the Valencia and 
Barcelona ports . Finally, in this first approach, we have chosen the Las Palmas 
port since, being an insular port, it is to be expected that the hinterland shared with 
the port of Bahia de Algeciras is nil. 

Following the same line of the previous model, the first Step to specify this 
model is to survey the series which are to be dealt with. As was the case 
previously, we can observe on the graph a positive relationship between the mean 
and the Standard deviation of the variables expressed in levels, especially in the 
port of Bahia de Algeciras (see figures 4.1 and 4.5). This evidence of 
heterocedasticity suggests us again the employment of the logarithmic 
transformation in both series. In addition, again by merely observing the series, we 
can conclude logALG and logLasP are nonstationary series, as seen in figure 4.6. 
For all this, the SVAR model shall be specified with the following series 
VlogALG and VlogLasP. 




Fig. 4.5. Graph of the mean-typical deviation of Las Palmas port series 



In the case of Valencia, Barcelona ports and, to a lesser extent, Bilbao port, we must take into 
account that they really have a shared hinterland with the port of Bahia de Algeciras although 
indirectly, through the Dry Dock of Madrid. This shared hinterland is the central area of 
18 Spain. 

It is logical to assume that the heterocedasticity will be higher in a hub infrastructure than in a 
spoke one. The traffics of the former will undergone a stronger fluctuation depending on the 
phase of the cycle, while the spoke infrastructure will be more stable since it uniquely 
depends on its hinterland. 
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Fig. 4.6. Bahia de Algeciras and Las Palmas port series (standardised data) 



By merely observing the VlogALG and VlogLasP series, we cannot conclude that 
they are stationary so we must apply the unit-root tests. Table 4.4 sets forth the 
results of these unit-root tests for these series, which show that the nil hypothesis 
is again rejected, and therefore, the original series (logALG y logLasP) are just 
i(i). 



Table 4.4. Unit-root test for VlogALG and VlogLasP 



Test VlogALG 


t-stat 


MacKinnon 

p-values 


Augmented Dickey-Fuller based SIC 


-4.8840 


0.0007 


Phillips-Perron(Newey-West using Bartlett kemel) 


-9.0163 


0.0000 






MacKinnon 


Test VlogLasP 


t-stat 


p-values 


Augmented Dickey-Fuller based SIC 


-5.9147 


0.0001 


Phillips-Perron(Newey-West using Bartlett kemel) 


-6.5017 


0.0000 



Once the likelihood ratio tests have been applied to determine the optimal lag of 
the VAR and the lag exclusion tests, we have chosen the following 
specification: 

VlogALG - - 0.0060* VlogALG(-2) + 0.4472* VlogLasP(-2) + 0.0819. 

VlogLasP = - 0.0702*VlogALG(-2) - 0.3233 *VlogLasP(-2) + 0.1448. 
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Table 4.5 accounts for the results of the tests in order to determine the optimal lag in the 
model studied. The optimal lag selected by the criterion is marked with an asterisk. 
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Table 4.5. VAR optimal lag (VlogALG and VlogLasP) 








Lag LogL LR 


FPE 


AIC 


SC 


HQ 


0 33.18294 NA 


0.000201 


-2.834813 


-2.735627* 


-2.811448 


1 34.76988 2.741085 


0.000252 


-2.615444 


-2.317887 


-2.545349 


2 42.51476 11.96935* 


0.000181* - 


2.955887* 


-2.459958 


-2.839061* 


* indicates lag Order selected by the criterion 

LR: sequential modified LR test statistic (each test at 5% level) 

FPE: Final prediction error 

AIC: Akaike information criterion 

SC: Schwarz information criterion 

HQ: Hannan-Quinn information criterion 







Table 4.6. VAR Bahia de Algeciras and Las Palmas 




VLogALG 


VLogLasP 


VLogALG(-2) 


-0.006029 


-0.070259 




(0.12120) 


(0.11066) 




[-0.04974] 


[-0.63490] 


VLogLasP(-2) 


0.447254 


-0.323318 




(0.24665) 


(0.22520) 




[ 1.81331] 


[-1.43569] 


C 


0.081991 


0.144847 




(0.03998) 


(0.03650) 




[ 2.05073] 


[ 3.96795] 


R-squared 


0.140674 


0.127223 


Adj. R-squared 


0.058833 


0.044101 


Sum sq. Resids 


0.321877 


0.268328 


S.E. equation 


0.123804 


0.113038 


F-statistic 


1.718877 


1.530559 


Log likelihood 


17.68515 


19.86865 


Akaike AIC 


-1.223.762 


-1.405721 


Schwarz SC 


-1.076506 


-1.258464 


Mean dependent 


0.130969 


0.095923 


S.D. dependent 


0.127615 


0.115616 


Determinant Residual Covariance 




0.000193 


Log Likelihood (d.f. adjusted) 




34.54087 


Akaike Information Criteria 




-2.378406 


Schwarz Criteria 




-2.083893 



As can be observed, the first lag has been removed since the p-value associated 
with the Nil Hypothesis of TheWald’s Exclusion Test of the coefficients is 





4 An Approach to the Hub-and-Spoke Systems from SVARs Models 



57 



20 

0.6714 . The following Step is to check that Csh = 0 is within our model’ s 
confidence intervals. In this case, the intervals for an accumulated response of 10 
periods are (0.0135,-0.0261), so we can support this hypothesis. 

Below are set out the graphs of the Impulse Response functions with the long- 
run restriction of Csh = 0, Csh being Clasp/alg in this model. The graphs 
correspond to the assumptions of the model proposed since, as can be observed, 
Clasp/lasp, Calg/lasp and Calg/alg are clearly positive. 




Fig. 4.7. Impulse Response functions accumulated under the Csh = 0 restriction 



Finally, in this case, the analysis of residues again Supports the white noise 
assumption for the residues. We have obtained a p-value of 0,7961 under the 
orthogonalisation restriction proposed. 

4.5 

Conclusions 

The applications detailed in this paper make us be optimistic as regards the 
employment of the SVAR methodology in transport economy. There is clear 
evidence that this category of models may be a useful instrument not only for 
traffic predictions but also for planning inffastructures and reducing the 
congestion costs associated with them. We should not overlook that we are 
dealing with an economic sector - transport - in which , the inffastructure demand 
is growing at a higher annual rate than that of the supply, both public and private. 
This fact makes it necessary to have the widest number of complementary 
mathematical instruments which help to plan the sector. 



20 



Table 4.6 accounts for the results of the VAR model specified. 
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Among the positive aspects of this methodology we can stand out its compatibility 
with other methods to make economic decisions such as the cost-social welfare 
analysis or the economic impact methodologies. As far as its negative aspects are 
concemed, we find that this methodology is initially restricted to hub-and-spoke 
transport structures, since we cannot establish long-run restrictions to carry out the 
SVAR for other ways of transport which do not have a central infrastructure to 
distribute the traffic. For this type of transport infrastructures which do not have a 
hub and operate with a common traffic, other possibilities within cointegration 
analyses can be employed. 

The tests here exposed have shortcomings which force us to be cautious and 
wait for further works which should be orientated not only to the port activity but 
also to other sectors such as the airport System. 
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The OECD economies are undergoing an upward trend in health expenditure. The 
usual explanations for this tendency are based on such elements as the progressive 
ageing of population, the increasing use of technology in health procedures, and 
the fact that public health Services become more extensive (cover more people) 
year after year. 

All this has led govemments to try to curb health spending by reducing the 
degree of inefficiency of public health Services. Such inefficiency has been 
generally admitted by many studies carried out by experts in the area . 

In Spain, health spending is mainly public and about 60% of this expenditure 
goes to hospitals. This draws attention to a sector that has continuously exceeded 
budgets in recent years, requiring extra money and generating a continuous 
pressure on costs. 

There has been an attempt to deal with the Situation from the point of view of 
prospective budgeting, through establishing Program-Contracts. This attempt, 
according to recent papers by Gonzalez and Barber (1996) and Ventura and 
Gonzalez (1999), has led to a significant decrease of inefficiency in hospitals 
directly managed by the national govemment. Ventura and Gonzalez’ s study 
concentrates on technical efficiency by applying Data Envelopment Analysis. On 
the other hand, Gonzalez and Barber follow a double approach: First, they analyse 
global inefficiency on the basis of the estimation of a parametric cost ffontier. 



* I wish to thank professors Guillem Löpez-Casasnovas, Carlos Perez Dommguez and Jose 
Miguel Sänchez Molinero for their comments, which I have found very helpful. However, I am 
1 responsible for any possible mistake. 

Undoubtedly, the "Dunning report" (1991) from Holland is the most renowned. In Spain, the 
Congress entrusted a committee headed by Abril Martorell with a study, and their conclusions 
were recorded in a document known as “Abril’ s Report” (1991). 



P. Coto-Millan (ed.), Essays on Micro economies and Industrial Organisation 
© Physica- Verlag Heidelberg 2004 
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Second, they use Data Envelopment Analysis to quantify the technical component 
of hospital inefficiency. 

However, the allocative component of hospital inefficiency has not been given 
much attention in the literature. This is probably due to information problems 
conceming input prices. In spite of such problems, there have been some relevant 
papers on this matter, such as Eakin and Kniesner (1988) and Eakin (1991). These 
two papers use the same sample of 133 hospitals in the USA and concentrate 
exclusively on allocative inefficiency. Other papers, such as Morey, Fine and 
Loree (1990) and Byrnes and Valdmanis (1994), quantify allocative inefficiency 
for a set of hospitals in California by means of Data Envelopment Analysis. 

For the case of Spain, only Puig-Junoy (2000) has studied allocative efficiency 
in hospitals, using a sample of both public and private hospitals in the region of 
Catalonia. Puig-Junoy’s approach uses the Data Envelopment Analysis technique. 
This author finds an interesting negative association between the proportion of 
public funds in the total budget of each hospital and de degree of allocative 
efficiency. 

The present paper studies the degree of efficiency of the Spanish hospitals 
managed by the INSALUD , by quantifying the extent to which each of the two 
inefficiency components (technical and allocative) increases cost. We have used 
an econometric analysis based on the estimation of a System consisting of a cost 
frontier and a set of input share equations. 

This paper is organised as follows: section 5.1 presents the theoretical model to 
be used; section 5.2 describes the variables and data to be used in the estimation; 
section 5.3 analyses the results, and section 5.4 summarises the final conclusions. 

5.1 

The Model 

Empirical studies on economic efficiency have usually followed the frontier 
methodology, ever since Farrell (1957) set the bases for this approach . Farrell’s 
approach basically Starts ffom the assumption that firms do not follow optimal 
behaviour pattems (they do not maximise profits) due to inefficiencies in the 
management of the production process. 

When calculating a cost frontier, it must be taken into account that firms may 
fail in their attempts to minimise costs for two different reasons: On the one hand, 
firms may not manage to use the minimum quantity possible of each input for the 
obtained level of output; this causes technical inefficiency. On the other hand, the 
factor combination used by the firm may not be the cheapest one, given the input 
prices; this would lead to allocative inefficiency. Both behaviours are shown in 
figure 5.1 for the case of two inputs. 



2 

INSALUD is the public agency which finances and administrates the Spanish hospitals in 1 1 
~ out of the 1 7 autonomous regions existing in the country. 

Forsund, Lovell and Schmidt (1980), Bauer (1990), Greene (1993) and Coelli, Rao and 
Battesse (1998) are good surveys on this technique. 
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A firm operating at point B would present both types of inefficiency: allocative 
inefficiency, since the firm is using an input combination that is not optimal, and 
technical inefficiency, since both inputs are being used in quantities greater than 
those required by output level y 0 , when production process OR' is chosen. The fact 
that OR' is chosen, instead of OR, results in allocative inefficiency, while the fact 
that input combination B is chosen, instead of D, results in technical inefficiency. 

If the firm chose input combination D instead of B, there would be no technical 
inefficiency but there would be allocative inefficiency. This allows us to measure 
allocative inefficiency at point B by the ratio AI B = C D /C*, where C D Stands for 
the total cost at point D and C* Stands for the minimum cost, reachable only when 
both sources of inefficiency have been eliminated (point A). The measure of 
technical inefficiency is given by the ratio TI B = C B /Cr where C B Stands for total 
cost at point B. The product of these two measures ( TI B * AI B = C B /C* ) shows 
the proportion in which the actual cost is higher than the minimum cost. 
Therefore, we can state that C B /C* is a global measure of economic inefficiency. 

Economic efficiency can be measured by the inverse ratio, C*/C B . It shows the 
proportion in which the minimum cost is lower than actual cost. This ratio is the 
product of two measures: that of technical efficiency and that of allocative 
efficiency: 



er _ i 

c s ' ti b ai b 



1 1 
ti b ai b 



= TE b AE b 



According to Aigner, Lovell and Schmidt (1977), the cost incurred by a firm may 
be outside the minimum cost fimetion due to economic inefficiency (either 
technical or allocative) as well as to purely random elements (thus, the term 
"stochastic ffontier"). This is shown by the expression: 

C (y, w) = C* ( y, w, ß ) e v TI AI (5.1) 

where C (y, w) is the effective cost of the firm and C* ( y, w, ß ) is the minimum 
cost fimetion which depends on production, y, input prices, w, and a number of 
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Parameters ß, that must be estimated; the term e v is the random disturbance and 
the product TI AI, is the measure of economic inefficiency as stated earlier . 

In logarithmic form: 

ln C (y, w) = ln C* ( y, w, ß ) + ln AI + ln TI + v (5.2) 

While the random disturbance can make costs go up or down (u may be either 
positive of negative), inefficiency always increases costs, as shown by the fact that 
InAI and lnTI are both non-negative. 

Expression (5.2) could be estimated directly. We would only need to postulate 
a particular fimctional form for the cost frontier and an adequate distribution for 
each of the error terms (v, InAI and lnTI). Nevertheless, this procedure would be 
subject to strong objections. The basic problem here is that allocative inefficiency 
(the InAI term) is clearly related to input prices and this implies that estimators 
cannot be expected to be consistent. 

In order to avoid this problem, we could estimate (5.2) together with the input 
share equations. This will be the procedure followed in this paper. Next, we will 
proceed to specify the input share equations, taking into consideration the fact that 
there may be allocative inefficiency. 

The effective share of the j* input in total cost, Sj (y, w), is given by the ratio: 

Sj (y, w) - x jWj / C(y, w) (5.3) 

where xj represents the amount of the j th input used by the firm, and wj the price 
of that particular input. Now, if we take into account that C(y, w) = Sj XjWj , it 
follows that d C(y, w)/dwj=Xj , which allows us to write: 

31nC(y,w)_ XjWj 

ainwj C(y,w) ( 5 . 4 ) 



Hence, we may define the effective share of the j th input as: 
31nC(y,w) 



Sj(y,w) = - 



d ln w ; 



and taking into account expression (5.2), it follows that: 



(5.5) 



0 , ölnC*(y, w,ß) 5AI 

STy,w) = w K/ + 

Slnwj Ölnwj (56 ) 

since TI and u do not depend on input prices. 

Now, applying Shephard’s lemma, we can identify 51nC*(y,w,ß)/51nWj with 
the optimal share of the f 1 input; that is: 



Sj*(y,w,ß) = 51nC*(y,w,ß)/01nWj (5.7) 

This expression allows us to write: 

Sj(y,w) = Sj*(y,w,ß) + Sj (5.8) 

where the error term, Sj , would include the effects of allocative inefficiency plus a 
purely random element, that is: 



Subscripts have been removed with the aim of simplifying the expression. 
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Sj = dlnAI/dlnwj + (5.9) 

Sj can have both positive and negative values, showing over or under-utilization of 
the j th input with respect to the Optimum level. 

In short, we will estimate the following System: 

ln C (y, w) = ln C* ( y, w, ß ) + ln AI + ln TI + v 

Sj(y,w) = Sj*(y,w,ß) + sj j = l...n-l (5.10) 

It must be observed that the number of share equations is not n, but n-1. This is so 
in order to avoid the estimation problems derived from the fact that the sum of the 
n shares must be equal to one. 

The relationship between the error terms of the input share equations and the 
lnIA term in the cost function is usually approximated by a quadratic form. This 
quadratic expression accounts for the fact that both positive and negative 
deviations from optimality result in cost increases. 

Schmidt (1984) proposed a specification in which the sum of error squares was 
weighed through a positive semi-definite matrix of parameters, F, such that InAI = 
s’Fs. The problem here is how to chose such parameters. On this subject, 
Kumbakhar (1991) showed that the expression InAI = s'Fs cannot be taken 
without restrictions. More specifically, the integrability condition, 

d ln (s'Fs) / d ln Wj = Sj j = l...n (5.11) 

must be fulfilled. As shown by Kumbhakar (1991), this integrability condition 
implies that the parameters in the F matrix are exact functions of the cost frontier 
parameters when a translog frontier is considered. To reach such a conclusion, he 
uses a rather restrictive specification of the share equations. In essence, he regards 
these equations as deterministic; that is to say, the error terms, e/s, are only due to 
allocative inefficiency and do not have any random component. 

In a later paper, Melfi (1984) proposed simplifying F into an identity matrix. 
This would make allocative inefficiency take values systematically close to zero, 
since the absolute value of share deviations is lower than one. Faced with this 
problem, Bauer (1990) suggested weighing the identity matrix through a positive 
Parameter, which, depending on the final value it takes in the estimation, enables 
us to avoid this systematic tendency. This paper will follow Bauer's Suggestion. 
Hence, for us: 

lnAI = c£s, 2 

J=> (5.12) 

In order to carry out the estimation of the System, some assumptions must be made 
about the distribution of the disturbances: lnTI, v and Sj, taking into account that 
they are independent of one another. It will be assumed that v is independently 
and identically normal distributed with zero mean and variance a 2 v ; s is 
distributed independently and identically as multivariate normal with constant 
mean jli and constant covariance L, where s = (si...e n _i)' and ja = (pi,... jan_i)\ 

5 The exact relationship for the restrictive Cobb-Douglas function has been proposed by Schmidt 
and Lovell (1979). It has also been obtained for the translog function -Kumbhakar (1997)- 
although, given its complexity, it has not yet been employed in any empirical study. 
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The fact that the mean of the errors in the share equations may be different from 
zero, as Schmidt and Lovell (1979) suggested, can be interpreted as evidence of 
systematic mistakes in the choice of the input combination. If the error mean is 
shown to be equal to zero, deviations with respect to optimal share levels will be 
merely random. 

As far as technical inefficiency is concemed, we must propose a distribution 
over positive values. In this case, we have chosen an exponential distribution 
with constant mean 1/a, although results show a high correlation with those 
obtained when the chosen distribution is a seminormal . 

Once the errors of the input share equations have been estimated, InAI can be 
obtained from expression (5.12). Therefore, the likelihood function of the System 
(5.10) is given by the product of two density functions: that of the composed error, 
lnTI + v, and that of the errors of share equations, s, that is: 

f [(lnTI + InAI + v), e] = g (lnTI + v) h (e) (5.13) 

The maximum likelihood technique enables us to obtain a consistent and 
asymptotically efficient estimation of the cost frontier parameters and the 
allocative inefficiency of each firm. From the cost frontier residuals obtained, it is 
possible to estimate the mean technical inefficiency of all firms; particular 
estimates of technical inefficiency for individual firms can also be obtained using 
the procedure suggested by Jondrow et al. (1982). It must be observed, however, 
that the estimates of technical efficiency obtained by such a procedure are 
unbiased but non-consistent. 

5.2 

Definition of Variables and Functional Specification 

Most of the data used in this estimation has been obtained from the EESRI, a 
survey of hospitals carried out by the INE (the Spanish Institute of Statistics). All 
data refer to 1994. 

These data refer only to public hospitals managed by the Insalud. The very 
small hospitals (less than 80 beds) and the very large ones (more than 1000 beds) 
have been excluded from the sample. The reason for this is that small and large 
hospitals usually differ substantially in their endowments of technical equipment. 
They also differ widely in the type of cases attended. This fact could be a handicap 
for inefficiency estimation since the highest use of resources could be understood 
as inefficiency when, in fact, it could simply mean a more intensive attendance to 
very complex cases. A total of 67 hospitals have been studied. 

Since we are dealing with a cross-sectional analysis, a short-run cost function 
has been specified, which depends on the output level, the input prices and the 
amount used of the fixed input, Capital. 



The literature has mainly dealt with the following distributions: the seminormal, the truncated 
normal (only for positive values) and the exponential distributions; see, for instance, Greene 
^ (1993) for further details on this issue. 

In the case of a truncated normal distribution for positive values, the mean was not found 
different from zero. 
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5.2.1 Outputs 

The first difficulty encountered when defining the variables involved in the cost 
function is how to measure the final output of hospital activity. Hospitals produce 
health and health improvement is difficult to measure. Nevertheless, since we do 
need a measure, we have to use some proxies. Such proxies could be the number 
of inpatient days or the number of cases attended. We have chosen the number of 
cases rather than the hospital stays. The problem with the number of stays is that 
many of them are presumably unnecessary and cost minimization may often lead 
the hospital to reduce the average length of the stay in each case by administering 
treatment in fewer days, all of which causes a cost increase per stay. To sum up, 
the accommodation Services provided by hospitals may be considered as input in 
the “production of health" rather than an output. 

When discussing hospital activity, it is convenient to distinguish between 
inpatient and outpatient attendance. This provides us with two basic output 
measures: 

CASES : this variable accounts for inpatient hospital activity, and it is 
developed from the number of discharged patients adjusted by the casemix 
complexity. It is a weighed sum of discharges from the different hospital Services. 
The weighs have been carried out using the coefficients defined by the Insalud in 
the UPA . 

AMBU: this variable measures ambulatory hospital care: first and successive 
visits and emergencies with no hospitalization. In this case, weighing is also 
carried out with the coefficients fixed by the UPA. 

5.2.2 Input Prices 

Two inputs, labor and materials, have been considered. The variables that account 
for their respective prices are: 

WAGE: The wage measure used in this paper is the average eamings of all 
hospital employees. Insalud wages are equal for all hospitals within each 
particular labor category (doctor, nurses, etc.), and wage data come from the 
BOE 9 (the Official State Bulletin). 

According to the previous definition, differences in wages among hospitals 
basically show the existence of different mixes of labor inputs (that is, the fact that 
hospitals employ different proportions of doctors, nurses, etc.). In spite of 
everything, our wage variable can be regarded as a global measure of labor costs, 
which allows us to use it as a proxy for the price of labor. If instead of using this 
global measure, we were to use the "real" price of any particular labor category 
(doctors' salaries, for instance) we would observe identical values for all hospitals 
and that would make our estimation impossible. 



These coefficients have been determined by the Insalud in a study that estimated the relative 
consumption of resources of a one day stay in each Service, outpatient visits and emergency 
g attendance. 

21/1993 Act of the Presupuestos Generales del Estado (National Budget Account) for 1994, 
December 29 th . 
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PMAT: This variable accounts for the price of an aggregated input which includes 
the purchasing of varied health equipment, cleaning, food, energy supply..., all of 
which can be called materials. The diversity of components makes it difficult to 
estimate a particular price for this category. Hence, we have resorted to an 
approximation: the total spending in this category per day of stay . 

The variability of PMAT may be related to both differences in prices and 
differences in the casemix of each hospital. This is indeed a drawback of the 
PMAT measure. Nevertheless, it can be argued that most elements included in this 
category (for instance: expenditure in food, clothes, cleaning, heating...) are 
unrelated to the casemix complexity. 

5.2.3 Fixed Input 

BEDS: The amount of fixed input (capital) used by each hospital is approximated 
by the number of beds. The justification for this is that the number of beds in the 
hospital reflects the stock of Capital invested around them. 

In order to use this proxy properly, it is necessary to have a homogeneous 
sample of hospitals. Therefore, we have focused the analysis on general 
hospitals , leaving out big hospitals with expensive modern technologies, and 
small hospitals with less sophisticated equipment. This leads us to consider that 
the hospitals in the sample have invested the same amount of Capital per bed and, 
therefore, that the use of the bed variable seems adequate as an approximation to 
capital. 

The Statistical description of the variables defined above is summarised in table 
5.1. 



Table 5.1. Statistical description of the variables 



VARIABLE 


MEAN 


STD. ERROR 


MAX 


MIN 


COST 


5,784,800 


4,670,787 


20,545,979 


939,649 


CASES 


23,573 


17,342 


72,962 


989 


AMBU 


119,038 


92,075 


533,145 


7,530 


WAGE 


3,347,091 


102,420 


3,560,862 


3,034,687 


PMAT 


19,870 


12,925 


95,672 


12,088 


BEDS 


360 


242 


958 


80 



5.2.4 Functional Form 

As regards the functional form, a translog cost function has been chosen for two 
reasons. On the one hand, for its flexibility, since it hardly places restrictions on 
technological features frorn the beginning, but these restrictions can be tested 



Therefore, it is considered that the stays produced may be an adequate proxy of the amount of 
U factor employed. 

Hospitals specialized in particular health Services have not been included since they may 
show differences with respect to the rest for requiring either too much or too little capital 
12 investment. 

We should not forget that the translog function is the approximation to the real function in 
one point. 
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later. On the other hand, it enables us to introduce more than one output, thus 
accounting for the multiproduct character of hospital activity. 

Then, the following System has been estimated: 

ln(Cost) - oto + a! ln (CASES) + a 2 ln (AMBU) + ßi ln (Wage) + 

!4 8„ ln (CASES) 2 + !4 8 22 ln (AMBU) 2 + Si 2 ln (CASES) ln (AMBU) + 

14 y n ln (Wage) 2 + p„ ln (CASES) ln (Wage) + p 21 ln (AMBU) ln (Wage) + 

r| u ln (BEDS) + InAI + lnTI + v 

S] = ß, + yn ln (Wage) + p,, ln (CASES) + p 2 , ln (AMBU) + s, (5.14) 

13 

in which only the labor share equation has been included; the usual symmetry 
restrictions have been imposed, and the theoretical requirement of homogeneity of 
degree one in input prices has been achieved by dividing costs and input prices by 
the price of the second input, which yields: 

Cost = COST / PMAT, 

Wage = WAGE / PMAT. 

5.3 

Estimation and Results 

Maximum likelihood estimation has been carried out on deviations with respect to 
the mean. Therefore, the parameters of the first-order terms represent the 
respective cost fimction elasticities for the average hospital. All of them are 
positive and, except for the case of the AMBU variable, significant at 1%, as 
observed in table 5.2, which summarises the results obtained. 

On the whole, parameters are highly significant, except those associated with 
the AMBU variable and its cross products. This is due to the strong correlation 
between the AMBU and the CASES variables. However, the existing 
multicolinearity does not affect the precision of the parameters relevant for 
inefficiency; so, the conclusions still stand. 

The significance of the exponential distribution parameter (a) confirms the 
existence of technical inefficiency, which explains 60% of the total variability of 
the composed error in the cost fimction . 

As far as labor share is concemed, the mean error, p 1? is significantly different 
ffom zero, indicating a systematic overutilization of labor. The estimated value of 
Pi means that hospitals, on the average, spend about 16.8 percentage points more 
on labor than required by minimum cost considerations. 



The share equation of the remaining factor has not been included because it does not provide 
14 any additional information (the sum of both equations is one). 

The composed error variability is explained by the variance of each one of its two 
components: lnTI and v. When the exponential distribution parameter, a, approaches infinity, 
the technical inefficiency variance tends to zero. In such a case, the variability of lnTI+v 
would be entirely explained by the random disturbance variance. Then, we would conclude 
that there is no technical inefficiency in that case. 
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Table 5.2. Estimation results 


VARIABLE 


COEFFICIENT 


STD. ERROR 


constant 


-0.214891*** 


0.07496 


In (CASES) 


0.596417*** 


0.11039 


In (AMBU) 


0.052727 


0.06816 


y 2 ln (CASES) 2 


-0.166184 


0.14222 


y 2 ln (AMBU) 2 


-0.189554* 


0.12799 


In (Wage) 


0.549951 *** 


0.07594 


y 2 ln (Wage) 2 


0.090459*** 


0.01210 


In (BEDS) 


0.361394*** 


0.07057 


In (CASES) In (AMBU) 


0.143169 


0.13046 


In (Wage) In (CASES) 


-0.034781 *** 


0.00737 


In (Wage) In (AMBU) 


0.00872 7 


0.00840 


c 


3.002640** 


1.74996 


a 


12.97428*** 


4.59385 


a v 


0.063042*** 


0.01652 




0.168106** 


0.07587 


cr 8 


0.022380*** 


0.00246 


Log likelihood 




222.7361 


Akaike info criterion 




-6.171227 



*** significant at 1%; ** significant at 5%; *significant at 10%. 



The overutilization of labor implies underutilization of the other factor. Hence, we 
can state that the wrong input combination is been used, which is evidence of 
allocative inefficiency, and its Overall impact on costs can be calculated through 
expression (5.12). 

From these results, efficiency indexes, whose Statistical description is 
summarized in table 5.3, are estimated for each hospital. In short, mean technical 
efficiency reached by the hospitals is considerably higher than mean allocative 
efficiency (with a difference of eight and a half points). 



Table 5.3. Estimated economic, technical and allocative efficiency 





Economic ef. 


Technical ef. 


Allocative ef. 


Mean 


78.12 


92.74 


84.22 


Std. Error 


5.77 


5.26 


3.75 


Min. 


55.45 


70.86 


77.37 


Max. 


90.33 


97.99 


93.96 



These indexes can be interpreted as follows: if Insalud hospitals were as efficient 
as the best ones, then their cost would be, on average, only 78.12% of their actual 
cost. If they used the least amount of resources possible, their cost would be 
around 92.7% of their actual cost, and if they chose factors in the right proportion 
(depending on prices), their cost would be 84.2% of their current cost. 
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0.78 0.80 0.82 0.84 0.86 0.88 0.90 0.92 0.94 
Fig. 5.2. Allocative efficiency distribution 




i i i i i i — i — i i i i i i i i i i i i i i i i i i i i 

0.70 0.75 0.80 0.85 0.90 0.95 



Fig. 5.3. Technical efficiency distribution 



Figures 5.2 and 5.3 show the distribution of each of the two indexes among the 
hospitals. Both series are hardly correlative so it cannot be concluded that the most 
technically efficient hospitals also present the highest allocative efficiency. 

Comparing the present paper with earlier works on Insalud hospitals, we may 
conclude the following: our technical efficiency index of 92.7% is slightly lower 
than the indexes estimated by Gonzalez and Barber (1996). These authors use 
Data Envelopment Analysis and obtain estimates of technical efficiency for 1991, 
1992 and 1993 of 95%, 96% and 97% respectively. However, our result is higher 
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than the one obtained by Ventura and Gonzalez (1999), who also used a non- 
parametric analysis. Ventura and Gonzalez studied Insalud hospitals between 1993 
and 1996, and obtained values that go from 80.7% in 1993 up to 84.7% in 1996. 

As regards allocative inefficiency, we agree with Puig-Junoy (1999) in that this 
inefficiency represents a problem quantitatively more important than that of 
technical inefficiency. Puig-Junoy (2000) uses Data Envelopment Analysis on a 
sample including both public and private hospitals in the region of Catalonia and 
obtains a measure of allocative efficiency with a higher dispersion (mean 89.1% 
and Standard error 9.5). 

Table 5.4 shows the cost increase caused by both types of inefficiency. 
According to our previous results, allocative inefficiency has a much higher 
impact on costs than technical inefficiency (18.96% and 8.2% respectively). 
Hence, any study focussing exclusively on technical inefficiency ignores the 
quantitatively more important aspect of hospital inefficiency. 



Table 5.4. Cost increase due to economic, technical and allocative inefficiency 





Economic inef. 


Technical inef. 


Allocative inef. 


Mean 


28.76 


8.2 


18.96 


Std. Error 


10.52 


6.9 


5.25 


Min. 


10.69 


2.05 


6.43 


Max. 


80.31 


41.12 


29.25 



As shown in table 5.4, total inefficiency (both technical and allocative) accounts 
for a 28.76% increase in hospital costs above minimum cost. Previous estimates of 
global inefficiency in Spanish hospitals have shown a wide variety of figures 
ranging from 13% to 58%, -see for instance: Gonzalez and Barber (1996), 
Wagstaff and Lopez (1996), Lopez and Wagstaff (1993)-. 

5.4 

Conclusions 

The aim of this paper has been to quantify separately technical and allocative 
inefficiency in the Spanish public hospitals. Given that a separate estimation of the 
hospitals’ cost frontier is not adequate, we have resorted to a joint estimation of the 
cost frontier together with the input share equations. On the basis of this joint 
estimation, we have obtained two separate measures of technical and allocative 
efficiency. 

These measures allow us to conclude that total inefficiency accounts for a 
28.7% increase in hospital costs above minimum cost. Previous estimates of 
global inefficiency in Spanish hospitals have shown a wide variety of figures 
ranging from 13% to 58%, -see for instance: Gonzalez and Barber (1996), 
Wagstaff and Lopez (1996), Lopez and Wagstaff (1993)-. 

Perhaps the most important aspect of the present paper is that it provides an 
estimate of allocative inefficiency, which is missing in previous papers. The 
present paper also identifies allocative inefficiency as the main component of 
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economic inefficiency. More specifically, allocative inefficiency accounts for 65% 
of the total cost increase due to global inefficiency. 
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This study justifies the Contemporary importance of efficiency analysis. We put 
forward the theoretical concepts of technical, allocative, and economic efficiency. 
We then tackle problems of an empirical nature. In this way, different ways of 
measuring efficiency are presented, with their main disadvantages. To finish, a 
theoretical application is presented for haulage firms operating in Spanish roads in 
six different sub-sectors with panel data. 

6.1 

Introduction 

The production theory of the firm, in most basic Microeconomics textbooks, Starts 
from the premise that the manager behaves in an efficient manner, making the best 
possible use of the resources available. To begin with, there is the hypothesis that 
the firm, given the technology and the inputs available, produces the maximum 
output possible. This relationship is represented in the production function. 
Secondly, the first order conditions for cost minimisation, given the factor prices, 
are introduced through the cost function. Finally, there is the assumption that the 
manager chooses that quantity of output which maximises profits, as represented 
through the profit function. However, in reality it may happen that some or none 
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of these hypotheses are fulfilled. There are various reason for this, ranging from 
bad luck to the pursuit of objectives different to those outlined above. As a result, 
in the empirical estimation of these behavioural functions much research has been 
dedicated to testing firm behaviour by analysing the problem of inefficiency in the 
productive System. 

6.2 

Efficiency 

We now distinguish between the different types of efficiency dealt with in the 
economics literature. We can distinguish between at least three types of efficiency: 
technical, allocative, and economic. We refer to technical efficiency when the 
optimality condition is defined by the production function. This means that the 
exact amount of inputs necessary are used in Order to produce a vector of Outputs, 
and hence there are no redundant (i.e. more than were strictly necessary) inputs. 
To estimate technical efficiency we need data about the quantities of inputs and 
Outputs but it is not necessary to have data on prices. 

Another concept is that of allocative efficiency , which occurs when the 
productive inputs are used in the proportions which minimise costs, that is, when 
the ratio of marginal products is equal to the ratio of their prices. 

Finally, we talk about economic efficiency when both of the aforementioned 
kinds of efficiencies are achieved simultaneously. In Farrell (1957) this concept 
(he refers to it as productive or global efficiency) is defined as the product of 
technical efficiency and allocative efficiency (he refers to the latter as ‘price 
efficiency’). 

We use Farrell’ s original ideas in order to provide a brief introduction to these 
concepts of efficiency. To begin, consider a firm employing two factors of 
production (x l5 x 2 ) to produce a single product. Under conditions of constant 
retums to scale, the technology of the Company could be represented by a single 
isoquant, the unitary isoquant ( y= 1), which can be drawn as in Figure 6. 1 . 

According to the definition of technical efficiency used above, all the 
combinations of inputs lying on the curve y = 1 are technically efficient. The 
combinations below the frontier are unfeasible while the combinations which are 
above are technically inefficient. With reference to Figure 6.1, point A represents 
the Situation of a Company that employs (x^, x 2 °) of inputs to produce the unit y. 
Now, according to the isoquant, y = 1 can be produced with (OB/OA)x!° and 
(OB/OA)x 2 °. Altematively, with (xi°, x 2 °) an efficient Company can obtain 
(OA/OB)y of the product. This means that technical efficiency (TE) can be 
measured using the ratio OB/OA. The maximum value of this index is one, which 
would mean that the firm is operating on the isoquant and is thus technically 
efficient. When a Company is inefficient, then OB ^ OA (as can be seen in Figure 
6.1) and the index takes a value lower than one which informs us about the degree 
of technical efficiency achieved by the firm. 

The definitions of allocative and economic efficiency can also be illustrated 
using Fig. 6.1. Allocative efficiency reflects the ability of a firm to use the inputs 
in the optimal proportions, given their respective prices and the production 
technology. Hence, if the relative prices of the factors of production are available 
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it is possible to know the slope of the isocost line PP'. A Company that allocates 
resources efficiently will thus produce with a input combination given by the slope 
of the isocost where it is equal to that of the isoquant, represented at point C. At 
this point, the firm is producing the output at minimum cost, given the technology 
and the input prices. Hence, the measure of allocative efficiency (AE) is derived 
fforn the ratio OD/OB. 

It is also possible to compare the total cost of producing a unit y using the input 
combination at C with that of producing y using the combination at A. We can use 
the ratio OD/OA to derive a measurement of economic or global efficiency (EE). 
A producer is economically efficient if it minimises the cost of producing the 
output given the input prices and the technology. A producer is economically 
efficient if and only if it is technically and allocatively efficient. In addition, we 
can see that: 

OD/OA = (OB/OA) (OD/OB) 

That is to say: 

EE = (TE) (AE) 




Fig. 6.1 



Technical efficiency as described above could be expressed in different ways 
depending on whether we use the output or the inputs as a reference. For example, 
we could take the production quantity as given and calculate the efficient quantity 
of inputs (input-oriented technical efficiency). Altematively, we could take the 
quantities of inputs as given and calculate the efficient quantity of output (output- 
oriented technical efficiency). Only if the technology exhibits constant retums 
scale will the input-orientated technical efficiency index be the same that the 
output-orientated technical efficiency index (Färe and Lovell, 1978). Although 
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constant retums to scale have been assumed in the explanation above, Forsund and 
Hjalmarsson (1979) generalises this for variable retums to scale. 

It is important to point out that the concept of economic efficiency is very 
wide-ranging. For example, we could also define economic efficiency in terms of 
the objectives of the Company or sector under study. If the objective is to 
maximise profits, the Optimum is associated with a profit function, given the 
prices of Outputs and inputs. If the objective is to minimise costs, the Optimum is 
associated with the cost function, given the prices of inputs and the quantity of 
Outputs. Economic efficiency may also be dealt with from other perspectives. For 
instance, when the price of production is not equal to the marginal cost of 
production, there will be economic inefficiency. Here we are referring to 
deviations with respect to the behaviour of a competitive industry, which 
generates inefficient allocations of resources in the long term. 

6.3 

Classification of Efficiency Frontiers 

Measurements of Companies’ efficiency are obtained by comparing the observed 
values for each productive unit with the Optimum defined by the estimated 
frontier. 

The studies on fiinctions of production frontiers (or those of costs or profits) 
can be classified in the following way: 

i) Deterministic studies: the error term is composed of a single one-sided term, 
which captures inefficiency. 

ii) Stochastic studies: the error term is composed of two parts - a component to 
capture “noise” and a one-sided component which captures inefficiency. 

These give rise to two types of production frontier: deterministic or stochastic. 

Either of these two production frontiers can be estimated using specifications 
that may be: 

i) Parametric specifications: a concrete parametric function is specified (Cobb- 
Douglas, CES, Translog...) 

ii) Non-parametric specifications: no function is specified. 

In addition, to calculate efficiency measurements one can: 

i) Use methods of mathematical programming to constmct the frontier. 

ii) Specify a Statistical relationship to estimate the frontier. 

From here eight possible combinations emerge in order to carry out empirical 
investigations. In the different empirical investigations stochastic frontiers are 
usually used with parametric specifications (for example, translog functions) and 
deterministic frontiers with non-parametric specifications (for example, Data 
Envelopment Analysis). 

If we use deterministic frontiers (see for example Aigner and Chu, 1968), all 
the observed cases below the production frontier represent inefficiency. Then, if 
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we have cross-sectional data with N inputs used to produce a single output, it is 
possible to define the production fimction as: 



yi =/(*«) exp(-w,) 



where y t is the output of producer z, x t is the input vector of N inputs used by 

producer Uf(x) is the production frontier and u t > 0 (0 <e~ U[ < 1). It is assumed 
that the observations on u t are independently and identically distributed (iid) and 
that x is exogenous (independent of ü). The technical efficiency index 
TE t =exp(-w / ) willbe: 



TE. =-¥i— 
fix,) 

In words, the technical efficiency index is the ratio of observed output to 
maximum output. However, the main disadvantage of this definition is that it 
attributes all of the deviation from the frontier to technical inefficiency. 

Stochastic frontiers , proposed by Aigner, Lovell, Schmidt (1977) and Meusen 
and van den Broeck (1977), are also called frontiers of composed errors because to 
every frontier we add an error which is composed of two components: a random 
error, and another error which captures the degree of inefficiency in the Company. 
That is: 



yt =f(Xi) exp (v r ui) 

where v t is a random error with some Symmetrie distribution to capture the random 
effects of measurement error and exogenous shocks, and {/(xj exp (vj } is the 
stochastic production frontier. In this way, technical inefficiency is captured by 
the one-sided error component exp (-Uj), u t > 0. Hence, the technical efficiency 
index is defined as: 



TE: = ^ 

/(x,)exp(v,) 

However, while this model is capable of revealing the average efficiency of the 
sample, it is unable to yield estimates of technical efficiecny for each observation. 

A solution, proposed by Jondrow, Lovell, Materov and Schmidt (1982), was to 
specify the functional form of the distribution of the u t component and to derive 
the conditional distribution (zz/lvj+w*). 

At first, empirical studies using frontier models used cross-sectional data. 
Schmidt and Sickles (1984) noted three difficulties with cross-sectional stochastic 
production frontier models: 

i) It is necessary to make assumptions about the distribution of technical 
inefficiency (e.g., half-normal) and Statistical noise (e.g., normal) in order to be 
able to decompose the error term. It is not clear how robust one’s results are to 
these assumptions. 

ii) It is necessary to assume that the degree of inefficiency is independent of the 
regressors, although it is easy to imagine that if the firm knows the technical 
inefficiency level, this should affect its input choices. 



Given the objectives of this study, the analysis will follow a production function focus. 
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iii) Technical inefficiency of a particular firm (observation) could be estimated, 
but not consistently. 

Schmidt and Sickles (1984) formulated and estimated some methods with panel 
data that would overcome these disadvantages. If each producer is observed over a 
period of time, we can apply panel data techniques. A production ffontier with 
time-invariant technical efficiency can be written as: 

yu =f(*it) exp (v ir Uj) 

where there are T time periods. We can then consider different panel-data models. 
The fixed-effects model treats the technical inefficiency component (u t ) as fixed 
effects, that is to say, as non-random. Thus, it is possible to allow for technical 
inefficiency to be correlated with regressors. However, if we assume that 
inefficiency indices are randomly distributed with a constant mean, we have the 
random-effects model. In this model it is necessary to assume that the technical 
inefficiency is uncorrelated with the regressors. To test this assumption, Hausman 
and Taylor (1981) created a test between individual effects (technical inefficiency) 
and explanatory variables (inputs vector in the case of the production fimction). 

The development of the analysis of econometric models from the 1980s has led 
the way to applications to frontiers with panel data. The economic literature in the 
1990s up to the present has gone in at least two directions: a) The first comes from 
the increasing use of data panels with a dynamic character. In models of efficiency 
there is an attempt to estimate indices of efficiency which vary with time. Studies 
by Comwell, Schmidt and Sickles (1990), Kumbhakar (1990) and and Battesse 
and Coelli (1992) constitute and good example of the literature in this vein, b) The 
second direction has an empirical character. There have been efforts to compare 
research done on the same sector, using the same methods and concepts of 
efficiency, in different countries or groups of countries, and the most relevant 
conclusions are drawn. An excellent survey of stochastic frontier analysis since its 
beginnings was given by Kumbhakar and Lovell (2000). 

6.4 

A Theoretical Application to Goods Haulage Companies 
on Spanish Roads Using Panel Data 

In this section, we calculate technical efficiency for goods haulage firms on 
Spanish Roads, considering six sectors (Full Load Transport, Grouppage-Service, 
Reffigerated Transport, Crane Transport, Special Transport and International 
Transport) over the period 1994-1997. We have chosen a translog fiinctional form 
to estimate the production function. Thus, the econometric specification will be: 

Log GVA= ß 0 +ßi LogCI +ß 2 LogK +ß 3 LogL + ß 4 (LogCI) 2 + ß 5 LogCI LogK 
+ ß 6 LogCI LogL + ß 7 (LogK) 2 + ß 8 LogK LogL + ß 9 (LogL) 2 + ß 10 T +ß„ (T) 2 
+ ß i 2 LogCI T + ß i 3 LogK T + ß i 4 LogL T 

where GVA is the output, measured as the Gross Value Added of the different 
Companies calculated as the sum of the Gross Operating Surplus (GOS) and the 
compensation to employees, where the GOS is calculated as the difference 
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between the income and expenditures, excluding consumption of fixed Capital; CI 
refers to the intermediate consumption, which in stock accounting is called 
provisions; K represents Capital and is approximated by fixed assets dedicated to 
production; and L represents the number of employees. 

All these variables except the variable L appear in millions of pesetas and are 
considered in constant pesetas. To achieve this, the variables VAB and K have 
been deflated by the Consumer Price Index (CPI), and the variable CI has been 
deflated by the Price of Energy Index to take into account that a very significant 
part of such consumption refers to the fuel of transport vehicles. In addition, all 
the variables are expressed in the form of deviations from their mean values, so 
that (at the mean) the first-order coefficients of the production fimction indicate 
how production varies with respect to the inputs. 

In all the estimations a fixed-effects model was chosen, that is, a separate 
dummy is estimated for every firm and the constant term is suppressed. These 
effects will be interpreted as the indices of technical inefficiency for each firm 
(Schmidt and Sickles, 1984). According to this method, the most technically 
efficient firm takes the value one. To correct for endogeneity of inputs we have 
used an instrumental variables approach, using their lagged values as instruments. 
Though the generic expression of the production fimction takes into account a 
temporal trend to measure technical progress, it was finally decided not to include 
it since it either did not prove significant or, as in the case of international 
transport, it showed a coefficient of around 8%, which was not very credible. 

In Tables 6.1 to 6.6, we summarise the results of the estimations. The 
estimations are quite satisfactory for most of the distinct sectors of haulage by 
road. In Table 6.7 the scale parameter and the technical efficiency parameter is 
measured for each sector. As can be observed, in haulage by road there is a sector 
such as Groupage-Service with increasing retums to scale, since the scale 
parameter is significantly higher than 1 . In the other sectors there are decreasing 
retums to scale with an average scale parameter of 0.81, as in the case of Crane 
Transport and International Transport. 

In Figures 6.2 to 6.7 the mean of the level of technical (in)efficiency and other 
statistics are summarised for each sector. The arithmetic mean of the indicators of 
technical efficiency relative to the sectoral mean is 0.80. The economic 
interpretation of this indicator is the following: in the haulage by road sectors, on 
average, the production could have been increased by 20% during the years 1994- 
1997 without increasing the quantity of inputs. 

However, the average levels and the dispersion of the (in)efficiency show a 
great deal of sectoral heterogeneity. The average level varies from 92% in the 
sectors of Full Load Transport and Crane Transport to 44% in the sector of 
Groupage-Service. The sector in which the indicator of technical (in)efficiency 
presents the highest dispersion with respect to the average sectoral value is 
Groupage-Service. This sector also has an average level of efficiency that is quite 
low. The sectors with less dispersion in the grade of inefficiency of Companies are 
International Transport, Full Load Transport and Crane Transport, which are 
mostly sectors with a high level of efficiency. 

In developing this application we wanted to illustrate one of the many methods 
which have been developed and which are available for the empirical study of the 
relative technical efficiency of Companies using the methodology of efficiency 
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frontiers. As has been seen, the topic of relative efficiency 


is a fundamental part of 


the Microeconomics of the Theory of Production and 


it is a subject whose 


contribution is important in terms 


of its applications. 


There remains another 


important economic topic, which is 


the analysis of the causes of the differences 


between efficiencies across Companies of the same sector. 


This is needed in order 


to identify allocation problems and to guide sectoral policies. This is important 


because, as is obvious, inefficiency is not only costly for the productive unit in 
question, but for the economy as a whole. 


Table 6.1. Full load transport 




Coefficient 


t-Student 


ß. 


0.3364 


4.3769 


ß 2 


0.7789 


6.9367 


ß 3 


0.0761 


2.3175 


ß4 


0.5735 


2.4592 


ß 5 


-2.4978 


-2.4874 


ß6 


0.3555 


1.3071 


ßv 


1.2029 


1.9595 


ßs 


0.6404 


1.4513 


ß9 


-0.1562 


-1.4568 


Table 6.2. Grouppage-service 




Coefficient 


t-Student 


ßi 


0.3854 


6.7586 


ß2 


0.8036 


10.744 


ß3 


0.0527 


1.2349 


ß4 


0.3198 


2.9957 


ß5 


0.2315 


0.7376 


ß6 


-0.7004 


-3.6098 


ß 7 


0.2552 


0.5308 


ßs 


0.2529 


0.5032 


ß9 


-0.5415 


-5.7813 


Table 6.3. Refrigerated transport 




Coefficient 


t-Student 


ß. 


0.1771 


4.1177 


ß 2 


0.1993 


3.3153 


ß3 


0.7388 


6.6673 


ß4 


0.2408 


2.7204 


ß5 


-0.0363 


-0.1296 


ßö 


-0.8211 


-2.1997 


ß7 


-0.2916 


-1.9447 


ßs 


0.4617 


0.9723 
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Table 6.3. Continued 








Coefficient 


t-Student 


ß9 


0.0022 


0.0059 



Table 6.4. Crane transport 




Coefficient 


t-Student 


ß. 


0.0635 


1.6573 


ß 2 


-0.0367 


-1.3134 


ßs 


0.8261 


10.678 


ß4 


0.0432 


1.6237 


ßs 


0.0293 


0.4857 


ß6 


-1.0788 


-3.1990 


ß7 


0.0988 


2.4034 


ßs 


-0.4035 


-2.0788 


ß9 


1.8919 


4.9323 



Table 6.5. Special transport 




Coefficient 


t-Student 


ß. 


0.2041 


2.2994 


ß 2 


0.0079 


0.2239 


ß3 


0.5063 


3.7002 


ß- 


0.3246 


2.2737 


ß5 


0.4331 


0.7951 


ß* 


-2.0812 


-3.4231 


ß7 


0.1933 


1.6674 


ß8 


1.2879 


2.1186 


ß9 


-1.0610 


-1.8020 



Table 6.6. International transport 




Coefficient 


t-Student 


ß. 


0.2336 


5.9920 


ß 2 


0.0865 


3.7535 


ß 3 


0.4440 


10.449 


ß4 


0.2874 


1.2973 


ß5 


0.2846 


1.3142 


ß6 


-1.1295 


-2.6446 


ß7 


-0.2230 


-3.6819 


ßs 


-0.0015 


-0.0064 


ß9 


0.9369 


2.5462 
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Table 6.7. Summary of results 


Sector 


Scale Parameter 


zf 


Average 








Technical Efficiency 


Full Load Transport 


1.19 


2.75 


0.92 


Groupage-Service 


1.24 


8.92* 


0.44 


Refrig Transport 


1.11 


0.98 


0.80 


Crane Transport 


0.85 


3.85* 


0.92 


Special Transport 


0.72 


3.21 


0.80 


International Transport 


0.76 


25.44* 


0.89 



*Statistically significant different from one at 5% level. 




Observations 27 




Mean 


0.923007 


Median 


0.938476 


Maximum 


1 .000000 


Minimum 


0.721584 


Std. Dev. 


0.062701 


Skewness 


-1.411545 


Kurtos is 


5.126332 


Jarque-Bera 


14.05251 


Probability 


0.000888 



Fig. 6.2. Technical efficiency. Full load transport 



12 




Observations 22 




Mean 


0.437364 


Median 


0.398774 


Maximum 


1.000000 


Minimum 


0.364753 


Std. Dev. 


0.134262 


Skewness 


3.592300 


Kurtosis 


15.44274 


Jarque-Bera 


189.2370 


Probability 


0.000000 



Fig. 6.3. Technical efficiency. Groupage-service 
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Fig. 6.4. Technical efficiency. Refrigerated transport 



Observations 14 


Mean 


0.803987 


Median 


0.784973 


Maximum 


1.000000 


Minimum 


0.561372 


Std. Dev. 


0.098407 


Skewness 


-0.463339 


Kurtosis 


4.529376 


Jarque-Bera 


1.865338 


Probability 


0.393502 




0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 1.05 



Observations 24 




Mean 


0.921818 


Median 


0.949583 


Maximum 


1.000000 


Minimum 


0.668798 


Std. Dev. 


0.076595 


Skewness 


-1.704514 


Kurtosis 


6.021926 


Jarque-Bera 


20.75351 


Probability 


0.000031 



Fig. 6.5. Crane transport 
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Observations 10 




Mean 


0.801077 


Median 


0.793137 


Maximum 


1.000000 


Minimum 


0.691286 


Std. Dev. 


0.080353 


Skewness 


1.393794 


Kurtos is 


5.151078 


Jarque-Bera 


5.165742 


Probability 


0.075557 



Fig. 6.6. Technical efficiency. Special transport 




Observations 32 




Mean 


0.889072 


Median 


0.893369 


Maximum 


1.000000 


Minimum 


0.778647 


Std. Dev. 


0.041603 


Skewness 


-0.281336 


Kurtos is 


5.002921 


Jarque-Bera 


5.771054 


Probability 


0.055825 



Fig. 6.7. Technical efficiency. International transport 
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7.1 

Introduction 

The main aim of this article is to analyse the technical efficiency of different 
international airline Companies during the period 1992-2000. For a number of 
reasons this analysis is of particular value. Firstly, dealing with the degree of 
market competition, the main studies into the question of efficiency -Forsyth, Hill 
and Trengove (1986), Encaoua (1991)- use a period of analysis during which air 
transport was tightly regulated in many countries and the majority of national 
airline Companies were under public ownership. More recent studies, such as those 
published by Good, Lars-Hendrik and Sickles (1999) or by Inglada, Coto-Millän 
and Rodriguez-Älvarez (1999), analyse efficiency by comparing regulated periods 
with other periods of partial liberaiisation. In both studies, improvements in 
efficiency are documented as the liberaiisation process advances. Nevertheless, 



P. Coto-Millan (ed.), Essays on Micro economics and Industrial Organisation 
© Physica- Verlag Heidelberg 2004 




90 



P. Coto-Millän, V. Inglada-L., B. Rey and A. Rodnguez-Älvarez 



there are no studies available regarding the efficiency of air transport Companies 
that were operational dnring a complete period of liberalisation, through which we 
might discover whether there are trends towards increases or reductions in 
efficiency. 

In second place, when focussing on the ownership of the Companies, most of 
the authors agree that the organisational System of large Companies under public 
ownership diminishes incentives to increase efficiency. Economic models exist 
that consider the lack of inter-company competition to be one of the principal 
determinants for inefficiency and lack of productivity, since there are no 
incentives for them to behave efficiently (see, for example, Leibenstein (1966) or 
De Alessi (1983)). Indeed, at the international level, since the mid- 1980s the trend 
has been towards the privatisation of big Companies. Pioneers in this process, for 
instance, were British Airways, Singapur Airlines (SIA) and Japan Airlines (JAL). 
More recently, in Spain we have witnessed the privatisation of Iberia. In this 
article we intend to make a comparative behaviour study of the efficiency of these 
19 airline Companies, most of which are in private ownership. 

In the light of all these changes that have taken place in the air transport 
industry, it seems logical to pose the question of how this liberalisation process 
has affected Company efficiency. Consequently, this article undertakes an analysis 
of technical efficiency through the estimation of a production function for a data 
panel charting the development of 19 airline Companies in the period 1992-2000. 
The methodology proposed is that of Battese and Coelli (1995), which makes it 
possible to make a model of the technical inefficiency of Companies on the basis 
of a set of explanatory variables that vary over time, as well as to differentiate the 
components of that inefficiency from other factors that are purely random and 
beyond the control of airline Companies. 

The outline of the paper is as follows: in the first instance, we present a resume 
of the most pertinent institutional changes in the air transport industry; we then go 
on to review the methodology we used to make the estimations, based on Battese 
and Coelli (1995); thirdly, a description of the data and variables is provided; in 
fourth place, the economic specification and fimctional form chosen are 
introduced and, lastly, we provide the results of the estimation and the conclusions 
drawn from them. 

Main policy changes in international civil aviation 

The international air transport industry, in the United States and Europe, to be 
specific, has been undergoing profound institutional changes since the beginning 
of the 1980s. The Asian countries, on the other hand, recently initiated moves to 
liberalise their markets and are close to an ”open skies" policy. Although it is 
beyond the reach of this article to provide a detailed background description, it is 
nevertheless worthwhile running through the main transformations this industry 
has undergone, concentrating on the differences between the European countries 
and their American competitors. 

From the early 1950s, the legal structure of civil aviation policy had been based 
on complex Systems of bilateral agreements between countries. These varied from 
deals that introduced great flexibility (in relation to capacity and fares) to others of 
a far more restrictive nature, based on the designation of Companies and as a 
consequence, on pre-established capacity and fares. The USA was at the head of 
the field in its efforts to adopt competitive bilateral agreements and this took on 
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concrete form in the International Air Transport Competition Act in 1979 
(IATCA). In it, the aims of the United States were defined as the striking of 
bilateral agreements allowing for the multiple designation of Companies, free 
access of charter Companies, the elimination of restrictions on capacity and fares 
and single treatment regarding the access of national and foreign Companies to 
airport Services. So from that moment on the American air transport industry went 
into full liberalisation, although, dating from the early 1970s, serious efforts had 
already been taken by the administration to limit restrictions, resulting in 
important gains in efficiency in comparison with previous periods (Sickles, Good 
and Johnson (1986)). 

In contrast, the European countries maintained a more protectionist approach 
until the decade after. Suffice it to say Article 84 of the Treaty of Rome excluded 
both air and maritime transport from the Obligation of having to observe the 
general rules of policy conceming competition. On the other hand, most national 
Companies in Europe were under public ownership; this form of ownership 
reduces the efficiency of airline Companies (De Alessi (1983)) and, furthermore, 
encourages the subsidising of Companies when they suffer losses. 

The first effect in Europe of the application of the IATCA was the unilateral 
agreement signed between the United States and Holland, allowing each Company 
to freely determine the capacity of their aeroplanes without any govemment 
intervention. 

The changes in Europe, however, were slow during the 1980s. Great Britain 
and Holland were the first two countries to attempt to end restrictions on 
competition, negotiating a bilateral agreement in 1984 whereby any change that 
affected capacity or fares would be considered valid if neither of the countries 
expressed disagreement (principle of disapproval). But it was not until December 
31, 1987, that the Official Journal of the European Communities published four 
Council provisions (comprising two Regulations conceming the defense of 
competition, a Directive on fares and a Decision on bilateral air transport 
agreements) that were to constitute the first package of measures in what was 
known as “EEC air transport deregulation”. These provisions were later brought to 
completion with Commission Regulations governing the defense of competition in 
the field of air transport. . 

Subsequent to that date, two more packages of liberalising measures were 
brought into force. To be specific, the period of 1993-1997 coincided with the 
application in Europe of the third package of air transport liberalisation measures, 
which was brought to completion in 1997. From that moment onwards, any 
Company from any country in the European community was at liberty to run 
passenger transport Services between any two destinations within the Community, 
with absolute freedom to establish their own air fares. 

If we compare Asian countries with the USA and Europe, we can see that the 
markets are still very ffagmented and Asian Companies are currently still subject 
to very restrictive bilateral agreements (Oum and Lee (2002)). However, the threat 
posed by the creation of international alliances between American and European 
Companies, is leading these countries towards liberalising their air transport. One 
of the first measures was the introduction of second tier airlines, such as Silk Air, 



* A detailed description of the air transport liberalisation process may be found in Rey (2000). 
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Eva Airways, Japan Asia Airways, All Nippon Airways, Asiana, Sempati and 
Dragon Air which have increased market competitiveness, against the traditional 
national Companies that had been operating as big monopolies. In this fashion, 
there has been a gradual introduction of a “limited or partial open skies” policy, 
begun by Thailand, the Philippines and Indonesia. And, parallel to this, the 
national Companies of these countries have been privatised. 

In 1991 Korea signed a revised Version of its bilateral agreement with the USA, 
described as “open carrier designation, open route and double disapproval 
pricing”, making it one of the most favourable bilateral deals struck in terms of the 
Stimulus it provided to competition. 

The bilateral deal signed between Japan and the United States in March 1998, 
as well as the agreement the USA and China came to in April 1999, bears the 
same liberalising hallmark. These constitute Symptoms that indicate the advance 
of the liberalisation process. Nevertheless, the pressure exercised by the USA and 
Europe in their drive to go beyond the “fifth liberty” in the region is still 
encountering stumbling blocks in many Asian countries. 

7.2 

Methodology 

To use Koopmans’ (1951) definition, a Company is technically efficient if it 
reaches the maximum output possible for the inputs and technology used. Farrelfs 
pioneering work (1957) developed a methodology designed to measure the 
concept of relative efficiency. This methodology makes it possible to “construct” 
an efficient frontier with the best productive units observed, where the 
(in)efficiency of a unit is defined as an index that measures the distance between 
itself and that frontier. This pioneering contribution gave rise to a wealth of 
written literature devoted to the task of discovering the best way to construct an 
efficient frontier (employing both mathematical and Statistical methods). The first 
contribution in the parametric estimation of frontier fimctions came from Aigner 
and Chu (1968), who obtained a deterministic production frontier using 
mathematical programming methods as their starting point. Affiat (1972) and 
Richmond (1974) proposed their estimation using parametric techniques, 
involving the possibility of bringing Statistical inference to bear on the results 
obtained. 

The main drawback of the deterministic frontier is that any deviation from it is 
deemed to be inefficiency and, therefore, there is no consideration of the 
possibility that productivity may be affected by random exogenous shocks that, 
along with possible errors of measurement of variables or problems of poor model 
specification, make it inadvisable to interpret the term of error as an isolated 
measurement of inefficiency. From quite another angle, the parametric and 
stochastic frontier approach (Aigner, Lovell and Schmidt (1977) and Meeusen and 
Van den Broeck (1977)) Starts out from the assumption that the deviation between 
the observed output level and the possible maximum is comprised of two 



Altematively, non-parametric frontiers ean be calculated via techniques for analysing data 
envelopment (DEA) and mathematical programming. 
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components: a Symmetrie term of error to capture the effect of variables that are 
not under the control of the productive unit under analysis, errors of measurement 
in variables and other statistic noise; and a second term that is supposed to capture 
the degree of inefficiency, positioning the production level below the frontier’s 
maximum output, making it necessary, in consequence, to specify an asymmetric 
distribution for this second term of error. Schmidt and Sickles (1984) develop a 
production frontier estimation using panel data. Kumbhakar, Ghosh and McGukin 
(1991) and Reifschneider and Stevenson (1991) establish a stochastic frontier 
model that sets out to explain inefficiency through the use of a specific variables 
vector for different Companies and a random error through a second stage analysis. 
Later, Battese and Coelli (1995) developed^an alternative model to get round the 
inconsistencies in the second stage analysis. 

In this research, the Battese and Coelli model (1995) has been specified, as it 
provides the possibility of analysing the determinants of the development of 
technical inefficiencv of a nroductive unit in^ terms of a set of explanatory 
variables that, in addition, may vary in time. To be specific, delimiting the 
analysis to the estimation of a production frontier using a data panel, the model 
can be expressed in the following form: 

Y it = X it ß + (V it - U it ); i=l,...,N; t=l,...,T. (7.1) 

where: 

Yj is the production of the i-th firm; 

X it is a kx 1 vector of input quantities of the i-th firm in the t-th time period; 

ß is a vector of unknown parameters; 

the V it are random variables which are assumed to be iid. N(0 ,cj v 2 ), and 
independent of the 

Uit which are non-negative random variables which are assumed to account for 
technical inefficiency in production and are assumed to be independently 
distributed as truncations at zero of the N(m it ,a L 2 ) distribution; where: 



m it = z it 5, (7.2) 

where z it is a px 1 vector of variables which may influence the efficiency of a firm; 
and 8 is an 1 xp vector of parameters to be estimated. 

So, the special contribution of this methodological approach lies in the fact that, 
without resorting to a second stage analysis, it enables us to specify the time 
trajectory of the technical inefficiency of a Company in relation to a set of 
explanatory variables that may change over time. 



3 

For a review of this literature, see Kumbhakar and Lovell (2000). 

4 

For possible regressors in the equation of inefficiency, Battese y Coelli (1995) propose both 
the explanatory variables of the production function and any variable able to determine 
changes in inefficiency. 




94 



P. Coto-Millän, V. Inglada-L., B. Rey and A. Rodriguez-Älvarez 



7.3 

Description of the Data and the Variables Used 

The database used in this article is constituted by a data panel formed by 19 
international airline Companies, of which 6 are European (Lufthansa, SAS, 
Finnair, Spannair, Iberia, British Airways), 6 North American (American Airlines, 
United, Delta, Northwest, USAir, Continental) 1 Canadian (Canadian), 2 Mexican 
(Aeromexico and Mexicana) and 4 Asian ( JAL, Korean Air, Cathay Pacific and 
SIA). They are all large-scale Companies operating international flights. The 
observations are annual and cover the period 1992-2000. The basic reference is 
provided by ICAO statistics (Digest of Statistics ffom the International Civil 
Aviation Organisation) and World Air Transport Statistics, published by the 
International Air Transport Association (IATA). Using this source of information, 
the database was constructed with the aim of estimating a production function for 
the air transport industry, utilising the methodology described in the previous 
section. The production function to be considered is as follows: 

Y = F (L, K, N) (7.3) 

where: Y is the production that has been calculated for the number of tons- 
kilometre available; L represents the total number of workers in the airline 
Company; K represents the Capital estimated for the aeroplane capacity available 
(expressed as tons available per plane); and, finally, N indicates the average 
distance covered in a flight stage (total kilometres flown / number of departures) 
and gives us an approximation of a measurement of the routes that make up the 
network of each airline Company. Table 7.1 shows a descriptive analysis of the 
data. 



Table 7.1. Descriptive data analysis 



Variable 


Average 


Typical Deviation 


Minimum 


Maximum 


Production (Y) 


14,297,400 


10,623,000 


757,895 


40,237,300 


Labor (L) 


34,109.21 


25,613.07 


1,837 


87,586.56 


Capital (K) 


6,639.49 


4,699.45 


693.799 


17,290.34 


Average Distance (N) 


1,607.36 


799.55 


702.15 


4,764.49 


Charter (%) 


3.59 


10.53 


0.01 


69.86 



7.4 

Econometric Specification 

The fünctional form of the production ffontier adopted is a transcendental 
logarithmic production function. The choice was based on the flexibility this 
function displays in adapting to all kinds of productive technology without it being 
necessary to impose a priori restrictions on scale capacity (see, for instance, Oum 



The reason for using tons-kilometre available as a measure of output is that this combines 
passengers-kilometre and cargo/tons-kilometre. 
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and Yu (1996)). So the function representing the production of international 
airline Companies is represented as: 



ln Y it = ß 0 + £ ß j ln(X jit )+Üß jh ln(X jit ) ln(X hit ) + 

j=l j=lh=l 

T 

+ Sy T D t +v it -u it 

t=l 



(7.4) 



i = 1, 19 countries; t = 1, 9 years 

where Y it is output and X it is a vector referring to the inputs under consideration. 
Technical progress is incorporated via time dummies (D t ) that register the effect of 
variables and, varying over time, have an equal influence on all the Companies. 
Lastly, v it is the random error and u it represents the term of inefficiency and they 
are distributed following the assumptions defined in the previous section. 
Component u it is defined by means of the equation: 

u it = 8 0 + 8 j (Trend it ) + 8 2 (Trend 2 it ) + W it (1 



i = 1, ..., 19 countries; t = 1, ..., 9 years 
where W it is a random error. 

In equation (7.5), factors have been included that are assumed to have the 
capacity to influence technical inefficiency. In the first place, we introduced a 
trend which registers the effects of the time period used and which, as we 
mentioned above, corresponds in time with the total liberalisation of air transport 
in the countries selected. In addition, with the aim of lending flexibility to the 
effect of time on inefficiency, the variable trend was also included in the table. 



7.5 

Results of the Estimation 



Equations (7.4)-(7.5) were simultaneously estimated for Maximum Likelihood, 
using the Frontier 4.1. program (Coelli (1996)). The variables were taken on 
deviations ffom their averages (approximate fimctional form). Thus, the estimated 
translog function is a second order approximation to the real function in the data 
average and the first order coefficients are the production elasticities of each input, 
evaluated on the sample average. 

The results of the estimation are recorded in table 7.2. The first order 
coefficients of the inputs bear the expected positive sign and are significant. The 
time dummy variables introduced into the model, as already explained, aim to 
register the effect of the technical change that has taken place in the air transport 
sector. 
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Table 7.2. Estimation of the production function (Battese and Coelli, 1995) 


Variable 


Coefficient 


t-Statistic 


Production Function (equation 4) 


Constant 


16.2507 


392.7113 


L(L) 


0.0811 


2.5937 


L(K) 


0.9172 


32.0844 


L(N) 


0.4259 


19.0189 


L(L) L(L) 


-0.1335 


-1.3777 


L(K) L(K) 


0.4362 


3.7914 


L(N) L(N) 


-0.2383 


-3.1769 


L(L) L(K) 


-0.0919 


-0.8810 


L(L) L(N) 


-0.1725 


-2.3922 


L(K) L(N) 


-0.0445 


-0.7103 


D93 


-0.0267 


-0.7170 


D94 


-0.0688 


-1.7424 


D95 


-0.1092 


-2.6007 


D96 


-0.1242 


-2.9545 


D97 


-0.1106 


-2.4538 


D98 


-0.0492 


-1.2264 


D99 


-0.0814 


-1.8071 


D00 


-0.0398 


-0.8571 


Model of Inefficiency (equation 7.5) 


Constant 


0.3966 


5.1966 


T 


-0.1520 


-3.9147 


T 2 


0.0128 


3.3112 


a 2 


0.0173 


7.1838 


Y 


0.9931 


50.0521 


Log. F. Likelihood 


168.75 





Using expression (7.6) it is possible to understand how these time effects have 
affected the production frontier from one period to another: 

TC t+ u = Yt+pYt (7 .6) 

So, a positive TC value would imply technical progress, that is to say, an upward 
shift of the production frontier, making it possible to obtain a higher output 
quantity without incurring a greater consumption of resources, and vice versa. The 
values obtained from eq. (7.6) are presented in table 7.3 and figure 7.1. 
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Table 7.3. Technological change 


Period 


TC 


92-93 


-0.067 


93-94 


-0.042 


94-95 


-0.040 


95-96 


-0.015 


96-97 


0.013 


97-98 


0.061 


98-99 


-0.032 


99-00 


0.041 



Technical progress evolution 




period 



Fig. 7.1. Development of technological change for the period of 1992-2000 



The coefficients obtained reveal that, during a first period (1992-1996), the 
passing of time had a negative influence on air Company technology, even though 
development was favourable. This improvement was particularly evident after the 
period of 1997-1996 where the coefficients indicated discrete technical progress 
(with the exception of the period 1999-1998). These results bear out our initial 
hypothesis that a direct relationship exists between air market liberalisation and 
the productive improvement registered in the Companies that make up the sector. 

Moreover, the interpretation of the coefficients obtained for the inefficiency 
model is particularly pertinent to the objectives of this study. The coefficient 
linked with the time variable is negative and statistically significant, suggesting 
that air Company efficiency improved during the time period under consideration 
and that this coincided, as previously suggested, with the process of deregulation 
in the sector. The squared coefficient of this variable is also significant and bears 
an opposite sign, indicating the suitability of the model proposed, since it makes it 
possible to introduce flexibility into the effect of time on inefficiency. 

The estimated value of parameter y indicates that the proportion of variance of 
Ui t on total compound error is 99.3% and serves to warn us that it is a mistake to 
use average production fimctions in which differences in inefficiency are not taken 
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into account. Furthermore, in order to justify the methodology adopted here, a 
series of specification contrasts were carried out, based on the ratio of likelihood. 
Firstly, we tested out the suitability of the functional translog form chosen. In this 
first contrast (see table 7.4), there is a rejection of the null hypothesis that the 
Cobb-Douglas functional form is preferable to the translog function. We then went 
on to contrast and reject the hypothesis of non-existence of technical inefficiency 
in the term of error. The third null hypothesis, considering inefficiency not to be a 
function of the regressors used, was also rejected. 



Table 7.4. Specification contrasts 



Null Hypothesis 


Likelihood Ratio 


Number of 
Restrictions 


Decision (95%) 


Ho:ß jh =0 


45.54 


6 


REJECTION 


H 0 : y = 8 0 = 5 j = 8 2 =0 


27.49 


4 


REJECTION 


H 0 :8j =S 2 =0 


9.28 


2 


REJECTION 



Finally, the estimation made allows us to calculate the indices of technical 
efficiency for each Company and year via the following expression: 

ET it = exp(-u jt ) = exp[-(s o +8, (trendy ) + 8 2 (trend 2 it /Y it ))-W it ] (7.7) 

In this way, technical efficiency is calculated as the ratio of the production level 
obtained in relation to the achievable maximum (that is, when U it =0) given the 
amount of inputs and the technology. Its value fluctuates between 0 and 1, the 
latter being the most favourable case. A more detailed analysis of these indicators 
is found in figure 7.2 and tables 7.5 and 7.6. In figure 7.2 we see the favourable 
development of the indices of technical efficiency within the time period analysed 
and observe that they stabilise at the end of the period in values close to the level 
of efficiency. Tables 7.5 and 7.6 show the average indices of technical efficiency 
for each Company, classified, respectively, by area. It can be seen in table 7.5 that 
the German Company (LUFTHANSA) reaches the highest level of efficiency 
(0.952), closely followed by the American Company DELTA, both with an 
efficiency index of 0.839. IBERIA (0.767) and USAIR are among the most 
inefficient Companies (0.792). By area, the results displayed in table 7.6 show that 
the Asian and North American Companies are the most efficient, followed by their 
European counterparts. One possible explanation of these results, in the case of the 
Asian Companies in particular, would be found in the growing competition 
introduced into these markets ffom the end of the 1980s onwards. That was when 
the monopoly of the national Companies was broken and the market was opened 
up to tier Companies, whose interests have had a striking impact on air transport 
policy in Asia (Knibb (1993)). Also, Hooper (1996) describes how a number of 
competing Companies sprang up in various Asian countries and the effect this had 
on international traffic in Korea, the Philippines, China, India y Japan. 
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Evolution of technical efficiency in air 
Companies. Period 1992-2000 



1 T 




1992 1993 1994 1995 1996 1997 1998 1999 2000 

year 



Fig. 7.2. Development of technical efficiency 



For these reasons, the Asian air market is po werful and experiencing very rapid 
growth. But, beyond factors due to the introduction of competition, some Asian 
Companies display features that make them more efficient than American and 
European Companies. If we look at SIA or Cathay Pacific, for instance, both 
Companies have an established reputation for the passenger Service they provide, 
for the efficient collection and delivery of luggage, for keeping their planes clean, 
and so on. 

From a cost perspective, their labour costs are lower, and their aircraft are of 
the latest generation. 

Last, but not at all least, we must mention their Computer reservation System. In 
Asia, Abacus, the joint System established by six of the regional airlines, including 
Cathay Pacific and SIA, is already up and running in both Hong Kong and 
Singapore. Of equal signiflcance, Abacus has come to an agreement to establish a 
joint Company with All Nippon in Japan, which will give the partner access to the 
dominant market in Asia. 

The European CRS Systems, although established prior to Abacus, are still 
struggling to create a viable product, and appear to be a long way from having 
something concrete to offer the market. Again, there is a clear danger that the lead 
established by Asian airlines in this fleld might give them a competitive advantage 
that could ultimately enable them to control the market. 
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Table 7.5. Indices of technical efficiency 


COMPANIES 


Average 


Minimum 


Maximum 


LUFTHANSA (GERMANY) 


0.9522 


0.8728 


0.9916 


SAS (SACNDINAVIA) 


0.8631 


0.7305 


0.9498 


FINNAIR (FINLAND) 


0.8022 


0.5831 


0.9903 


SPANNAIR (SPAIN) 


0.9183 


0.8634 


0.9795 


IBERIA (SPAIN) 


0.7679 


0.6863 


0.9952 


BRITISH AIRWAYS (UK) 


0.8873 


0.7515 


0.9737 


AMERICAN (UNITED STATES) 


0.8997 


0.8177 


0.9607 


UNITED (UNITED STATES) 


0.9309 


0.8099 


0.9900 


DELTA (UNITED STATES) 


0.9435 


0.8932 


0.9856 


NORTHWEST (UNITED STATES) 


0.8564 


0.8009 


0.8948 


USAIR (UNITED STATES) 


0.7928 


0.7466 


0.8141 


CONTINENTAL (UNITED STATES) 


0.9309 


0.8531 


0.9828 


JAL (JAPAN) 


0.8951 


0.7655 


0.9908 


CANADIAN (CANADA) 


0.8623 


0.7239 


0.9786 


AEROMEXICO (MEXICO) 


0.9533 


0.8957 


0.9930 


MEXICANA (MEXICO) 


0.8194 


0.8046 


0.8395 


KOREAN AIR ( KOREAN REPUBLIC) 


0.8843 


0.8034 


0.9765 


CATHAY PACIFIC (CHINA) 


0.8998 


0.7183 


0.9903 


SIA (SINGAPORE) 


0.9147 


0.8776 


0.9719 



Table 7.6. Technical efficiency per Company 


Companies 


Average 


Minimum 


Maximum 


Total sample 


0.8838 


0.5831 


0.9952 


European Companies 


0.8641 


0.5831 


0.9952 


North American Companies 


0.8916 


0.7239 


0.9930 


Asian Companies 


0.8975 


0.7183 


0.9908 



7.6 

Conclusions 

It is customary in economic theory to link processes of deregulation and 
liberalisation with incentives directed at improving the productivity and efficiency 
of Companies. In this context, the objective of our research was to explore the 
relationship between the liberalisation of the air traffic market that took place in 
the period under study (1992-2000) and the productive change that was brought 
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about in Companies in the sector. To this end, a study was made of the behaviour 
of 19 international airline Companies (European, Asian and American), utilising 
the stochastic frontier model designed by Battese y Coelli (1995), which makes it 
possible to estimate the development of technical progress and technical 
inefficiency, in addition to introducing variables related to liberalisation as 
determinants of technical inefficiency. 

The results obtained in this investigation bear out the hypothesis we started out 
from, that a direct relationship exists between market liberalisation of airlines and 
the productive improvement that has been seen in Companies in this sector. 
Moreover, the time period in which this deregulation process was carried out had a 
significant influence on improvements in the levels of technical efficiency, 
American and Asian Companies being most efficient of all. Nevertheless, in the 
case of the latter there were factors that had a positive effect on efficiency, such as 
the good reputation and quality approach of some Asian Companies, coupled with 
advantages in labour costs and the fact that their SCR is much more competitive 
than that of the Europeans. Consequently, within the international market, the 
Asian air market is amongst those with the greatest potential for future growth, 
while American and European Companies are naturally investigating ways to 
increase their share in these markets. 
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8.1 

Introduction 

One of the main structural changes in all OECD countries has been a shift in 
employment and total output towards Services. This trend indicates that Services 
are a dynamic part of the economy and make an increasing contribution to 
employment and economic growth. Thus it could be argued that the “service or 
tertiary society” is the stage towards which all industrialised countries are moving. 

Although the share of Services in the economy has increased, there is a lack of 
understanding about their industrial performance. Traditional views considered 
that employment in Services was low-skilled and unproductive, while tertiary 
activities were considered neither dynamic nor innovative. In fact, several of the 
Service activities require highly skilled jobs, show increasing productivity and are 
highly innovative and dynamic. Although Services have been classified as non- 
tradable, they are increasingly exposed to competition and are becoming more 
tradable. The increasing importance of Services in the so-called “industrialised 
countries” and the poor understanding of the process of Services expansion require 
a review of the interpretations of the expansion of Services. 

The primary interpretation of the growth in Services was the “theory of stages” 
(Kindleberger, 1958, Rostow, 1960) whose explanation was based mainly on the 
pattems of final consumption (Petty’s Law and, specifically, Engel’s Law: as 
income per capita increases, final demand shifts towards superior goods including 
Services. This reflects a shift in consumer demand due to a high income elasticity 
of Services). Furthermore, the growth of the service sector generated diverse 
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interpretations about the u post-industrial societies” (Bell, 1974), according to 
which the Service sector is gradually taking the place of industry as the new engine 
of growth. However, a different set of questions were raised about the growth in 
Services with Baumol’s (1967) “theory of unbalanced growth”. Amongst the most 
important questions posed were those relating to the definition of the role of the 
pattem of consumption and the role of productivity differentials in relation to 
Service output and employment growth. Fuchs (1968) showed that the pattem of 
consumption had a less important role than that of productivity differentials 
(relatively slow productivity growth in some Services). As a consequence of the 
lower productivity growth of the Service sector as a whole compared with the 
manufacturing sector, plus the low skilled labour intensive characteristics of many 
Services (“cost disease of personal Service”), a secondary set of questions appeared 
when the intersectoral comparisons are made in constant or current prices 
(Gershuny & Miles, 1983 and Kravis et al ., 1983, Baumol et al 1989). In general, 
the increase in the contribution of Services to GDP during the last five decades has 
stemmed more from changes in relative prices than from an increase in output. In 
Order to understand this trend, however, it is necessary to distinguish between 
different Services. 

One of the problems encountered when examining the “Services” is how they 
can be most accurately defined. Since the Service sector is highly heterogeneous, it 
must be broken down into different categories according to the functions 
performed, the transformation processes and to the market served: producer or 
consumer Services; distributive, social, personal or business Services; market or 
non-market Services; and physical, person-centred or information Services. Most 
of these branches have been incorporated gradually into the Systems of National 
Accounts to improve the Classification of Services for analytical purposes (CEC- 
EUROSTAT, IMF, OECD, UN & World Bank 1994 and OECD 1995a & B). In 
order to understand the shift towards the Service economy, it is necessary to 
analyse what is being produced in the economic System and how it is being 
produced. 

Classification of Services is not the only difficulty, traditional analysis is too 
biased on manufacturing performance. The process that drives Services activities 
in many cases is different to that of manufacturing. In the same way, the key 
factors in Services have a different relevance to those in manufacturing (R&D, 
innovation, organisational change or human Capital). Notwithstanding, several 
Services activities are becoming similar to manufacturing and vice versa (Standard 
process and mass production), while the differences among Services are as varied 
as those among manufacturing. In fact, the focus on the different categories 
between manufacturing and Services is becoming less interesting and it seems 
more relevant to analyse the interaction between sectors and activities (Tomlinson 
1997). 

Standard indicators of technology intensity show that Services make a 
contribution to total R&D expenditures which is relatively limited compared with 
the size of the sector in total employment. Certainly there is a problem of 
measurement, since most countries have only recently covered Services in R&D 
and innovation surveys because it was assumed that manufacturing was the source 
of technological change and traditional measures do not usually capture key 
factors of innovation in Services (such as patent registration). 
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One advantage of macroeconomic analysis based on Input Output techniques is 
that it enables overall production to be disaggregated by sector and by sub-system. 
It then becomes possible to examine the growth of different Services and industries 
in relation to the process of structural change of the economic System. 
Furthermore, since the evolution of a sector or sub-system is not independent of 
the rest of the economy, it is necessary to evaluate the links between sectors 
(intersectoral relations) in terms of changes of final demand and technical change 
(Diaz Fuentesl993). 

With these general premises in mind, and with the specific methodological 
approach, the following sections of this paper aim to explore a number of key 
questions: 

- What has been the extent of structural change in employment in Services in the 
main “industrialised countries” over the past four decades? 

- Which technological trends explain the directions of structural change in 
employment in manufacturing and Services? Was there a relationship between 
direct technology intensity and employment growth by sector? 

- What is the embodied contribution of R&D and innovation expenditures to total 
technology intensity in manufacturing and Services? What is the sectoral 
importance of the acquisition of technology embodied in inputs generated by 
other industries? 

- Is the growth of employment in Services explained solely by the relative 
increase in the final demand of these Services, or is it also necessary to consider 
the growth of intermediate demand of Services? Has the growth of Services 
been caused by a greater demand for their use in the production of 
manufactured goods? 

While accepting that there is a positive correlation between economic growth and 
employment in Services, this research departs frorn an optimistic point of view 
about Services. This seeks to analyse the increasing integration between 
manufacturing and Services in technological innovation. This integration includes 
the outsourcing of manufacturing to specialised Services activities and Services 
that were not performed previously by other firms. 

The following sections discuss these questions and demonstrate the importance 
of detailed analysis of inter-industry performance. Part 8.2 considers the trends in 
employment in the main “industrialised countries” taking a global measure of the 
extent of structural change of the major sectors share frorn 1960 to 1980, and frorn 
1980 to 1997 (OECD, 1992a & annual a). Part 8.3 discusses the relationship 
between conventional measurements of technology intensity and the trend’s 
directions in employment in manufacturing and Services since 1980 (OECD, 
annual b and 1998). Part 8.4 goes beyond examining the extent and direction of 
structural change in Services and provides an inter-industry analysis of innovation 
which, on the basis of the most recently published input-output table, estimates 
technology flows and the acquisition of innovation by sector in 1994 (INE 1997a 
and 2000, EUROSTAT 1997). Part 8.5 presents the results of the changing 
composition of employment in Services due to final and intermediate demand 
(Garcia et al. 1994, EUROSTAT 1987, 1992 and OCDE 1995c). These results 
have been obtained through the use of inter-industry analysis. This analysis is 
founded on the notion of vertical integration of the sub-systems and provides an 
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enhanced view of the intersectoral relationships. Finally, some implications of the 
results are drawn in part six. 



8.2 

Extent of Structural Change in Services Employment 



Düring the last four decades the world economy has exhibited extensive economic 
structural changes. At the same time, economic performance has varied 
significantly; the real increase in OECD GDP and productivity during the period 
from 1980 to 1997 was half that experienced between 1960 and 1980. The slow- 
down of economic growth in OECD countries since the end of the 1970s has been 
accompanied by three significant recessions, two oil shocks, growth of 
international trade, globalisation of financial markets and the diffusion of a set of 
new technologies (Freeman & Soete, 1994). As a consequence of these changes, 
the economic structures of these countries have been transformed markedly, 
reflecting structural, as opjposed to cyclical, shifts in the composition of 
employment and production. Notable economists consider that the reasons for 
this slow-down and higher unemployment are structural and are either caused by 
restrictions in markets, or by the lack of technological innovative capabilities. 
Certainly, despite the considerable economic growth in the postwar period (GDP 
growth rates for OECD countries averaged around 4 per cent between 1960 and 
1980 and 2 per cent between 1980 and 1990), labour absorption has been limited 
throughout the three decades, since employment average growth rates were around 
1 per cent in both periods considered. This trend has been even more significant 
for EU countries, whose GDP growth averaged 4 per cent (1960-1980) and 2 per 
cent (1980-1997) while total employment growth rates were only 0.2 and 0.5 per 
cent for the respective periods. Considering the problems of labour absorption 
during the “golden age”, it is interesting to explore employment trends since the 
tuming point of the middle of the 1970s (OECD, annual a). 

Taking a global measurement of the extent of structural change of the major 
sectors shares for the main industrialised countries including Spain (Table 8.1), it 
is clear that there are basic similarities in the pattems of structural change among 



"Whilst economic theory has pointed to compensation mechanism generating new employment 
to replace jobs which are lost through technical change, no one has claimed that this process is 
instantaneous or painless. Economists differ however on the extent to which they would rely 
on self-adjusting market-clearing mechanism or on active public investment and labour market 
2 policies". See Freeman & Soete (1994, pp. 17-38). 

Over the last decade the European unemployment rate has averaged 10 per cent, which is a 
much more serious matter than the fluctuations around the average. Conventional business 
cycles account for relatively little of the history of unemployment. Most of the annual 
variations in unemployment come from the long-frequency fluctuations between half decades 
rather than from the short-ffequency fluctuations within half decades. This is because there are 
long term changes in social institutions, and the shocks (wars, oil or financial crisis) have long- 
^ lasting effects, see Layard et al. (1994, pp. 91-109). 

These Information and Communication Technologies, although they have a vast ränge of 
present and future applications, do not yet easily match the inherited previous skill profile, 
management Organisation, industrial structure or institutional framework. See Freeman & 
Soete (1994, pp. 47-66). 
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countries during the whole period: agriculture is declining while Services are 
increasing as a share of overall employment and in all the three areas the sectoral 
changes were shifting towards Services. However, in the period 1960- 80 the rising 
share in employment Services was mainly correlated to the declining share of 
agriculture rather than in manufacturing or industry as in the period 1980-97. 
Additionally, countries differ widely in the sectoral composition of employment, 
in the proportion of structural adjustment, and in the degree of flexibility which 
work Organisation displays in response to changes. On the one hand, fforn 1960 to 
1980, some countries showed an increase in structural changes and industrial 
employment (Japan, Spain and Italy) and GDP shares (Japan and Spain), implying 
a significant catching-up in comparison with the leader country. On the other 
hand, from 1980 to 1997, the contribution of Services to GDP has declined in the 
two countries with higher productivity growth (Japan and Germany). Furthermore, 
structural change can be considered as a source of growth. This applies, in 
particular, to countries in which employment is high in agriculture and 
productivity 4 is low, since labour can be reallocated to other sectors of higher 
productivity. 



Table 8.1. Structural trends in employment 1960-97 

Agriculture Industry Services Wholesale Transport, Finance, Community, 

& retail storage & insurance, Social & 

trade, commun - real Personal 

restaurants ications estate, Services 

& Hotels business 



Services 





1960-80 1980-97 1960-80 1980-97 1960-80 1980-97 


1980-97 


1980-97 


1980-97 


1980-97 


USA 


-4.9 


-0.9 


-4.0 


-7.4 


8.9 


8.3 


1.2 


-1.3 


7.5 


1.1 


Canada 


-7.8 


-1.5 


-4.2 


-5.3 


12.0 


6.8 


1.1 


-0.9 


3.4 


3.5 


UK 


-2.1 


-0.7 


-10.1 


-10.7 


12.2 


11.4 


0.3 


-0.7 


4.6 


1.9 


Germany 


-8.4 


-2.4 


-2.9 


-7.6 


11.3 


10.0 


2.9 


-0.6 


1.4 


3.0 


France 


-13.8 


-4.2 


-1.7 


-10.3 


15.5 


14.5 


1.7 


-0.7 


4.2 


10.1 


Italy 


-18.3 


-7.5 


4.0 


-5.9 


14.3 


13.4 


3.3 


0.6 


0.8 


11.5 


Spain 


-19.5 


-10.8 


5.9 


-6.2 


13.7 


16.9 


5.7 


-0.1 


1.9 


7.6 


Japan 


-19.8 


-5.1 


6.8 


-2.3 


13.0 


7.4 


-1.1 


-0.6 


1.7 


6.9 


Australia 


-4.5 


-1.3 


-8.0 


-8.8 


12.5 


10.1 


1.8 


-0.6 


5.5 


1.7 



Source: OECD (annual a & b). 



Within Services, the share of finance, insurance, real estate and business Services 
(FIREBS) and community, social and personal Services (CSPS) increases 
proportionally in all cases while that of transport, storage and Communications 
(TSC) and, in certain cases, Wholesale and retail trade, restaurants and hotels 



4 

There is some evidence, based on 30 countries for the period 1960-90, that the GDP growth 
rates correspond negatively to increasing Services shares. See Chenery (1986). 
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(WRTRH) lose their share in total Service employment. Breaking down these 
figures shows that, from 63 to 85 per cent of growth was due to Service sectors 
and that among these activities FIREBS and CSPS have the greatest impact. The 
increasing share of business and other intermediary Services makes it gradually 
more difficult to measure compositional change, since it becomes necessary to 
look at structural links among sectors instead of showing only change in 
employment or GDP shares by sector. 

8.3 

Innovation and Employment Trends in Structural Change 

The Classification of the main sectors into three parts is a conventional but limited 
method of measuring the extent of structural change. A more accurate approach 
would be to disaggregate the sectoral evolution of GDP and employment in 
relation to the different trends in each individual sector, so it becomes possible to 
distinguish growing, medium and declining growth activities. A complementary 
measurement of structural change in terms of direction could be obtained by 
classifying the branches according to their technology intensity. 

Contemporary theories of economic growth and international trade have stressed 
the role of innovation as a fundamental source of growth, employment and 
productivity, the capacity to innovate depends on multiple factors (Grossman, G. 
& Helpman, E. 1994). Technology Investments are developed in a few 
manufacturing industries, however, the overall performance of the economic 
System depends on putting technology to work by using ideas and products 
developed in other activities. However, most firms and industries, in particular 
Services such as FIRBS, acquire technology by purchasing and assimilating Capital 
embodied technology machinery. This fact has changed the EU attitude to Science, 
technology and innovation: ” Support to innovation should be broadened from 
“mission-orientated” projects with specific research outcomes, such as a new 



In OECD (annual b) the industries were classified according to their annual growth rate (1974- 
90) in the main industrialised countries as: High-Growth: 1. Computers and Office machinery, 
2. Aerospace, 3. Communications, 4. Finance and insurance, 5. Business Service, 6. 
Government, 7. Rubber plastic, 8. Pharmaceutical, 9. Social and personal Service, 10. 
Instruments. Medium-Gr owth: 11. Chemical, 12. Trade, 13. Transport, 14. Agriculture, 15. 
Electrical machinery, 16. Paper and printing, 17. Electricity, gas and water, 18. Non-ferrous 
metals, 19. Food, drink and tobacco, 20. Motor vehicles, 21. Hotels and restaurants. Low- 
Growth: 22. Mining, 23. Non-electrical machinery, 24. Construction, 25. Fabricated metals, 
26. Stone, clay and glass, 27. Textiles, 28. Petroleum refining, 29. Wood and fumitures, 30. 
ß Ferrous metals, 3 1 . Shipbuilding. 

In OECD (annual b) the 21 manufacturing branches across 11 industrialised countries are 
ranked according to the R&D expenditures to gross output as a rough estimate of technological 
sophistication, with the following scheme: High-Tech: 1. Aerospace, 2. Computer and office 
machinery, 3. Communication equipment, 4. Pharmaceutical, 5. Instruments, 6. Electrical 
machinery. Medium-Tech: 7. Motor vehicles, 8. Chemical, 9. Non-electrical machinery, 10. 
Rubber and plastic, 11. Non-ferrous metals, 12. Other transports. Low-Tech: 13. Stone, clay 
and glass, 14. Food, drink and tobacco, 15. Shipbuilding, 16. Petroleum refining, 17. Ferrous 
metals, 18. Fabricated metals, 19. Paper and printing, 20. Wood and fumitures, 21. Textiles, 
footwear and leather. 
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combat aircraft, to “diffusion-orientated” programmes, such as educating small 
firms about new products and process" (Working Group on Innovation and 
Technology Policy 1999). 




Fig. 8.1. R&D intensity and employment growth 1980-95 

Source: OECD (annual b & c, 1998) 



Technology intensity concems the degree to which technology is produced and 
used within different industries. Activities with relatively high R&D or innovative 
expenditures per unit of output or value added are classified as high technology 
industries or technology intensive industries. However, technology generation and 
use are hard to measure, the roughest estimation is R&D intensity that measures 
the expenditure directly incurred by an industry, while technology use is estimated 
by R&D expenditures incurred by the acquisition of intermediate and Capital 
goods. While the shares of manufacturing in employment and output have 
declined over the last two decades, the share of high technology industries in 
manufacturing output has increased steadily. However, employment shares of high 
technology industries have risen less than output shares, which means that labour 
productivity has risen in these industries. Fig. 8.1 provides some evidence of R&D 
effort and relative employment growth in the manufacturing industries for the 
period 1980-1995 for 15 OECD countries. The graph shows a correlation between 
high technology industries and employment. This indicates that technological 
change has accompanied structural change, favouring the emergence of 





110 



D. Di'az-Fuentes 



employment in high technology sub-sectors such as Aircraft, Office and 
Computing machinery, Drugs and medicines, and Professional goods. 

In the Services, the evidence shows that the most rapidly growing sectors in both 
output and employment terms are FIREBS and CSPS; in certain cases WRTRH 
has noted a significant expansion (in particular in countries specialised in 
tourism). Although the Service sectors are increasingly recognised as being 
important innovators, the existing indicators of technological intensity do not 
reveal their significance (Andersen & Howells 1998). Innovation in Services is 
less based on R&D and more linked to acquired technology and changes in 
process, markets and Organisation. In the past many countries focused R&D 
surveys on manufacturing because it was assumed that this was the source of 
technological change and innovation (Young 1996). 

The latter approach would be usefiil in helping to evaluate the direction of 
structural change, but would fail to identify the transformations taking place 
between industries below the aggregate levels. The economies are also undergoing 
a different structural change, as the flrms change the Organisation of their 
production and source inputs. These changes affect the linkage between and 
within flrms, industries and sectors. These sorts of structural changes can be 
analysed by looking at the structure of production in each industry and sector. 
These can be done applying input-output techniques, which provide a picture of 
inter-industry relations and linkages. Additionally, it could not connect the change 
in the structure to other factors such as shifts in domestic demand, foreign trade, 
technical change or input productivity. 

8.4 

Embodied Technology and Technology Diffusion in 
Services 

Input-output techniques make it possible to analyse the economic System based on 
intersectoral relations of innovation expenditures that represent the technology 
flow of embodied R&D or innovation. This methodology applied to technology 
diffusion allows the measurement of total innovation intensity by sectors (Sakurai 
et al 1996). 

In order to understand the shift towards the Service economy it is important to 
examine the principal factors affecting both what is produced in the whole 
economy and how it is produced. This means that, when examining the Service 
sector, it is not enough to look at this sector alone, since structural changes in the 
pattems of final demand, the intermediate demand and technical change must also 
be considered. 

Services form a network through which economic activity takes place. The 
process of structural change can influence the degree to which Services are 
required throughout the whole economic System. The supply of Services, in tum, 
makes it possible to attain greater specialisation and division of labour. These 
factors are essential in reinforcing the observed shift towards Services 
employment. 
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In order to present the methodology, a sector will be defined as a cluster of 
branches (firms) producing Commodities in agreement with a Standard 
Classification (NACE), and a sub-system as a group of different activities which 
are required in the economic System to produce a specific product or Service. In 
the context of inter-industry relations, a sub-system consists of activities of 
different branches, all of which directly and indirectly contribute to the production 
of a specific final output. 

The balance equation of output in an input-output table can be defined as: 
x = Ax + f (8.1) 

X = LF (8.2) 

where x is the vector of output by branch, A is the matrix of technical coefficients 
whose typical elements ay = Xy / E Xy represent the value of the ith input needed to 
produce one unit of industry j's output; L =(I-A) _1 , known as the Leontief inverse, 
represents the total requirements per unit of final output in terms of gross output; 
X is the diagonalised vector x; and F is the diagonalised vector f of the final 
demand by industry. The intersectoral relationships presented in equation (8.1) for 
an open static economic System can be written as: 

X = A d x + f d +e 

X = L d F (8.3) 

where A d is the matrix of domestic coefficients, f 1 is the domestic final demand 
vector for domestic Outputs, e is the foreign demand or exports vector and the 
inverse L d =(I-A d ) _1 represents the total requirements per unit of final domestic 
output whose typical element is l d j. 

Technology intensity for each industry can be defined as the R&D or 
innovation directly incurred expenditure per output for industry i : = rj x l9 

whose diagonalised vector is R. Thus, the vector of technology or innovation 
embodiment ti can be expressed in diagonalised form as: 

TI d = R L d F (8.4) 

This equation connects innovations to final domestic and foreign demand and the 
multiplier of total innovation embodiment per unit of final demand is mi = E rj l d j. 
Finally, in an open economy, imports are a source of technology diffusion, while 
production for exports induces demand for imported inputs. 

m = A m x + f" 1 = A m L d F + F m = M(4) 

TI m = R* [A m L d F + F m ] (8.5) 

As was noted in the third section, conventional indicators of R&D or innovation 
intensity such as direct R&D per unit of output or value added, distort the 
measurement of total R&D and innovation of industries. Indicators that estimate 
total innovation embodiment per unit of output or (such as mi = E r, l d j and R m 
A m ) are more appropriate measures. 
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Table 8.2. Re-estimation of embodied technology intensity indicators by sector 1994 




R&D 


Total R&D 


Innovation 


Total Innovation 




direct 


intensity 


direct 


intensity 


High technology sectors 


28 Communications 


15,9 


16,5 


80,6 


83,1 


13 Other transport n.e.c. 


36,2 


40,8 


55,3 


64,8 


11 Electrical Machinery and electronic 


24,4 


27,1 


37,3 


43,5 


equipment 










12 Motor vehicles 


8,4 


12,5 


31,2 


43,5 


Medium technology sectors 


21 Rubber & Plastic Products 


5,1 


9,1 


15,5 


25,7 


20 Publishing, Printing & Reproduction 


0,4 


3,6 


12,7 


24,5 


or recording media 










30 Finance and insurance 


0,0 


3,1 


0,0 


13,6 


intermediation 










10 Office and Computing Machinery 


5,3 


8,5 


7,0 


12,5 


24 Recovery and repair 


0,0 


4,7 


0,0 


12,4 


1 7 Leather products and footwear 


0,6 


2,9 


3,4 


11,5 


29 Electricity , gas & water supply . 


2,8 


4,2 


5,2 


10,5 


Lower technology sectors 


1 Agriculture, fishing 


0,0 


2,3 


0,0 


8,0 


23 Construction 


0,0 


2,2 


0,0 


7,3 


26 Hotels & restaurants 


0,0 


1,4 


0,0 


7,0 


15 Tobacco products 


2,0 


2,6 


4,5 


6,9 


27 Transport 


0,0 


2,3 


0,0 


6,4 


33 Public Services 


0,0 


1,9 


0,0 


5,0 


32 Private social and personal Services 


0,0 


1,3 


0,0 


3,4 


31 Business Services & real state 


0,0 


0,6 


0,0 


2,8 


25 Wholesale & retail trade 


0,0 


0,7 


0,0 


2,7 



Sources: Elaborated by the author based on R&D and innovation data (INE 1997 & 1999) and Input- 
Output tables for 1994 (INE, 2000). 



Considering that R&D and innovation direct expenditures are reduced in most 
Services with the exceptions of Communications and electricity, gas and water 
distribution, the Classification by sectors is blurred when considering the purchase 
of R&D and innovation inputs. The change is of particular significance for R&D 
users such as financial and insurance intermediaries and recovery and repair, and 
to a lesser extent for the rest of the Services. In terms of innovation, table 8.2 
shows similar results but the contribution of acquired innovation is relatively 
larger in sectors such as construction, hotels and restaurants and transport. 

R&D expenditures are widely considered as indicators of technology intensity 
and the main determinant of economic growth and productivity. In Spain as well 
as in other “industrialised countries” these expenditures mainly originate in a few 
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manufacturing industries such as: electronic machinery and electronic equipment, 
Chemical products (including pharmaceutical) and other transport (including 
aerospace) (Papaconstantinou et al. 1996). A large share of the Outputs of these 
“high technology” manufacturing sectors are demanded as intermediate inputs into 
the production process by different sectors and are also sold to final demand. In 
this way the direct R&D expenditures of the provider industries become embodied 
in products, process and Services across the economic System. 

The estimation of performed and acquired technology is estimated using input- 
output techniques and R&D and innovation expenditures. The shares of R&D and 
innovation are expressed in relative terms to the total indirect expenditures 
embodied in output using the methodology explained above (equations 8.4 and 8.5 
for TI). 



Table 8.3. Technology clusters of R&D and innovation. The 7 largest sectors in Spain in 1994 



R&D providers 


share 


R&D users 


share 


11 Electrical machinery and electronic eq 


18,05 


24 Recovery and repair 


6,8 


7 Chemical products 


15,92 


21 Rubber & Plastic Products 


6,4 


13 Other transport n.e.c. 


12,09 


12 Motor vehicles 


6,1 


12 Motor vehicles 


11,85 


7 Chemical products 


4,5 


28 Communications 


11,01 


32 Private social and personal Services 


4,1 


9 Machinery n.e.c. 


6,90 


4 Basic Metal Ferrous 


4,1 


29 Electricity , gas & water supply. 


3,45 


30 Finance and insurance int 


3,4 




79,28 




35,3 






Total Services 


28,3 


Innovation providers 


share 


Innovation users 


share 


28 Communications 


17,31 


1 Agriculture, fishing 


6,2 


12 Motor vehicles 


13,76 


24 Recovery and repair 


6,0 


14 Food & beverages 


13,23 


12 Motor vehicles 


5,7 


7 Chemical products 


10,43 


26 Hotel & restaurants 


5,6 


1 1 Electrical Machinery and electronic eq. 


8,59 


21 Rubber & Plastic Products 


4,8 


13 Other transport n.e.c. 


5,74 


30 Finance and insurance int 


4,6 


9 Machinery n.e.c. 


4,14 


32 Private social and personal Services 


4,5 




73,21 




37,4 






Total Services 


35,1 



The picture that emerges in terms of R&D is, on the one hand, a concentrated 
cluster of “high technology” industries providing most of the R&D in the 
manufacturing sector and in certain Services such as communication and 
electricity, gas and water distribution. On the other hand, the cluster of users is 
more disperse and includes different Services and manufacturing activities, such as 
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recovery and repair Services, private and personal Services and finance and 
insurance institutions. 

In terms of innovation expenditures the picture that emerges is slightly 
different. The cluster of providers is still concentrated and includes some of the 
same industries than in R&D but also others such as food and beverages. The bulk 
of innovation is in communication and motor vehicles rather than in “high 
technology” manufacturing providers. Last but not least, the cluster of innovative 
users includes Services that are not usually considered “high or medium 
technology intensive” such as: agriculture, forestry and fishery; recovery and 
repair Services; hotels and restaurants; finance and insurance; and private social 
and personal Services. Thus, many Services industries act as the main users of 
technology and innovation and constitute significant R&D and innovative clusters. 

8.5 

Compositional Structural Change in Empioyment (Final 
and Intermediate Demand for Services) 

In the second section the extent of the structural change has been examined and, in 
the third section, the direction of these changes in the face of technology 
intensities. A more precise definition of compositional structural change considers 
changes in the sectoral integration of an economy: output and empioyment shares 
reported for different sectors, and the changes in the inputs used by them. The 
advantage of this method is that it provides a detailed image of how the structure 
of an economic System and its linkages are at one moment, and how they have 
unfolded over time. The evaluation of both the extent and direction is connected to 
the broad sources of change for each sector, namely: final demand, import 
Substitution and pattem of inter-industry linkages in the economy (referred to as 
technical change), and this represents the path of change followed to reach a 
specific sectoral structure. 

The Input-Output (IO) technique enables changes in output and empioyment to 
be estimated, and it is also usefixl in helping to evaluate the relationship between 
empioyment and technology. The temporal variations of IO intermediate 
coefficients themselves reveal significant information about the technical change 
that operates in an economic System. In so doing, they represent an extension of 
previous measurements of structural change. A complementary IO analysis of 
empioyment and innovation is based on the concept of the sub-system and the 
notion of a vertically integrated sector, which was introduced by Sraffa (1960) and 
Pasinetti 1981) for theoretical purposes, but can also be used in applied terms 
(Sakurai 1993, Diaz Fuentes 1999). 

The intersectoral relationships between branches and sub-systems from an IOT 
will be represented in a single table as matrix operator B. The last two terms on 
the right hand side of (8.2) L d F correspond to the actual amount of all domestic 
input that is directly and indirectly required for the production of a final 
commodity (column j). When X is multiplied by the matrix L d F the operator B is 
obtained: 

B = XL d F 



( 8 . 6 ) 




8 Technological Innovation and Employment 



115 



This is the matrix of shares of production or B operators, the elements of which 
(by) shows the share of total output Xi which is required in the sub-system j. The 
sum of these elements is one. The results of matrix operator B can be utilised to 
re-analyse variables associated with the production by branch such as 
employment, R&D and innovation. This can be disaggregated to the highest level 
(Terleckyj 1974, Barker 1990 and Sakurai 1997). In this case, the diagonalised 
vector of employment ü by branch has been used to calculate the matrix of 
employment U: 

U = üB 

U shows, by rows, the amount of employment that each branch contributes to each 
sub-system (and of which the sums are the same total of u), and the columns of U 
show the employment of each sub-system. These matrices B and U disclose the 
direct and indirect shares of output and employment by sector and sub-system, and 
the indirect shares of output can be separated by replacing the Leontief inverse in 
(8.3) and (8.6) by L - I. With this methodology, employment in Services can be 
separated by final and intermediate requirements. 



Table 8.4 



Employment in intermediate and final demand market Services 
(thousands of employees and percentage of total employment in Services) 





1975 


1980 


1985 


1989 


1994 


Intermediate 


1.089 


1.213 


1.408 


1.506 


1.643 




26,0 


30,0 


33,4 


31,1 


31,5 


Final 


3.103 


2.836 


2.810 


3.329 


3.566 




74,0 


70,0 


66,6 


68,9 


68,5 




4.192 


4.049 


4.218 


4.835 


5.209 


Employment in intermediate market Services in manufacturing as a 
employees and percentage of manufacturing sub-system) 


sub-system (thousands of 


Intermediate market Service in Manufacturing 


402 


390 


378 


407 


418 


% of manufacturing Subsystem 


16,1 


18,0 


19,1 


22,1 


28,1 



The general results that have been extracted can be summarised by the following 
two main points. First, table 8.4 (top part) shows that the growth of market Service 
employment in Spain was due mainly to an increase in intermediate demand; and 
this cannot be explained exclusively by the “stages of growth theory”. This 
corresponds to the trends exhibited by the main advanced European countries that 
experienced constant growth in the share of employment due to intermediate 
market Services since the 1960s. Therefore the growth of market Service 
employment in the European economies is not directly, nor mainly, due to an 
increase in final demand, but rather to an increase in the intermediate demand of 
Services. Second, table 8.4 (bottom part) shows that the expansion in intermediate 
demand for Services is accompanied by their increasing use in the production of 
manufactured goods in Spain, and this corresponds with the trends exhibited by 
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the principal European economies throughout the considered period. This implies 
that the production of manufacturing goods goes beyond the industrial sector and 
requires increasing Services. 

8.6 

Conclusions 

This chapter has examined the expansion of employment in Services that 
constitutes one of the most significant features of long term structural change in 
“industrialised economies”. Four issues regarding structural change in 
employment in terms of innovation have been presented: structural change 
extension; direct technological intensity effects; embodied innovation; and the 
inter-industrial dimension in terms of the changes in the final and intermediate 
demand. 

The most rapidly growing sectors in terms of employment in most 
“industrialised economies” have been Services and in particular: FIREBS and 
CSPS. However, the increase in employment in WRTSRH has been significant in 
Spain, France and Italy, which have been at the same time the countries with the 
largest increases in Service employment between 1980 and 1997. The observed 
pattems of the extension of structural change in employment do not indicate 
whether the growing role of Services reflects a change in final demand, business 
demand or outsourcing to specialised Services sectors. 

Assuming the importance of technological innovation in international economic 
growth and productivity, and the fact that a few manufacturing industries 
concentrate most of the R&D surveyed expenditures, the correlation between 
direct technology intensity and employment was evaluated. It was observed that 
technological change determined structural change in employment in sub-sectors 
such as aerospace, office and computing machinery, telecommunications 
equipment, drugs and medicines and Professional goods. Moreover, some Services 
sectors such as FIREBS and CSPS acquire technology by purchasing and 
assimilating inputs embodied in innovation related to the previously mentioned 
high technology intensity manufacturing sub-sectors. However, it is notable that, 
although Services are increasingly innovators, the existing indicators of R&D and 
technological intensity do not reflect the whole scope of innovation in Services 
because traditional surveys have been designed for manufacturing. 

One important element of innovation is the acquisition and not only the direct 
expenditures on R&D. To identify the inter-industrial flows of innovation taking 
place between sector below the aggregate levels input-output techniques were 
applied, in first instance, to analyse the embodied R&D and innovation 
technological flows with the purpose of measuring “total technology intensities” 
by sectors and sub-systems. Folio wing a defined methodology it was noticed that 
the acquisition of technology was an important component of the innovation 
expenditures in Services. Given the low level of direct technology intensity in 
Services, the total innovation and R&D embodied intensity increased significantly 
in all Services categories, but in particular in: finance and insurance, recovery and 
repair, electricity, gas and water distribution, hotels and restaurants and transport, 
which correspond with the aggregated trends envisaged in the second and third 
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section. Moreover, given that a limited number of manufacturing sectors are the 
main providers of technology, the estimation of performed and acquired 
technology based on input-output techniques shows a concentrated cluster of high 
technology industries including communication Services and a cluster of R&D and 
innovation users that includes recovery and repair, private and personal Services, 
finance and insurance, and hotels and restaurants. Thus, many Services act as the 
main users of technology and constitute a technology cluster which importance 
should be reconsidered in national survey on innovation (OECD 1995c). 

This result points to the importance of policies of innovation diffusion in 
Services. Although new technologies are concentrated in a small number of 
manufacturing industries that spend directly in R&D, the new process, products 
and Services created in that industries generate benefits that become widespread 
through diffusion and use. The performance of an economic System depends on 
applying technology by using and adapting products, process and Services 
generated elsewhere. This ability of the firms and industries is critical for the 
economic system's productivity and growth. 

Finally, a more precise analysis of structural change was presented in the fifth 
section. This considered not only the extension and the technological direction 
component but also the broad sources of structural change from the intermediate 
and final demand side. Growth in Service employment must be explained by 
considering the increasing Integration between sectors (services and industry). 
This relationship can be explained by the following four factors. First, the 
increasing specialisation among sectors, which requires a complex network of 
Services such as Communications, transport, banks and insurance, recovery and 
repair, and after-sales Services that link the different sub-systems. Second, the 
expansion in the foreign trade of goods and Services, which is another perspective 
of specialisation. Third, the augmenting regulations (Standards of quality or 
environment) which require specialised Services (such as legal, tax, engineering, 
Publicity, training, accounting, or finance and insurance). Fourth, the emergence 
of new economies of scale in the production of Services, which induce a process of 
extemalisation or outsourcing to specialised Service sectors. 

The application of intersectoral analysis in this research on Spain and other 
European countries has generated a different conclusion from that derived by the 
aggregated analysis, since the methodology allows the interpretation of the links 
between manufacturing and Services, and the sets of innovation expenditures and 
employment that are directly and indirectly utilised in the productive Systems. 
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The assessment of the relevance of the new theories of trade as opposed to the 
traditional theoretical approach, when it comes to explain real trade pattems, is 
essentially an empirical question. Nevertheless, there are two problems (which 
remain unsolved) which influence the results obtained when measuring intra- 
industry trade and specialisation. The first problem refers to the very existence of 
the phenomenon and is related to the definition of “industry” and the selection of 
the level of data dissagregation more appropriate to study such a phenomenon. 
The second problem, closely connected to the former, comes from the objective 
difflculty of finding a convenient quantitative measure. 

This chapter is devoted to the analysis of the obstacles faced by empirical 
treatment of intra-industry trade. With this purpose, we present first in 9.1 a 
general critical overview of the main measurement indices proposed by the 
literature. We deal next in 9.2 with two important problems of measurement: the 
adjustment of global trade imbalance and the question of categorical aggregation. 
Finally, in 9.3 we offer a summary of the main conclusions reached regarding 
these aspects of intra-industry trade. 

9.1 

Measures of Intra-industry Trade and Specialisation 

Intra-industry trade may be defined as simultaneously importing and exporting 
products belonging to a particular country and the same industry. This 
phenomenon comes usually accompanied by intra-industry specialisation, which is 
the concentration of production factors in particular groups of products inside an 
industry, at the expense of other production lines. Although both expressions are 
offen indiscriminately used, they refer to two different aspects. 
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Therefore, at a conceptual level, it is possible to establish a distinction between 
intra-industry trade and intra-industry specialisation. According to that distinction 
and following the Classification provided by Greenaway and Milner (1986) and 
Kol and Mennes (1983, 1986), we can differentiate two broad categories of 
indices regarding empirical measurement of the aforementioned concepts: 

i) Indices determining the extension of intra-industry trade by measuring the level 
of overlap existing in trade fiows. 

ii) Indices of intra-industry specialisation, which are measuring the degree of 
similarity of trade pattems (relative structure of imports and exports). 

We next present separately the main indices used for the measurement of both 
aspects. Later, a comparison is made among the diverse measures and a Suggestion 
is proposed as to which may be more appropriate according to the purpose of the 
analysis to be performed. Finally, we summarise in subsection 9.1.4 the later 
developments on the measurement of intra-industry trade and mention the current 
unresolved issues on the subject. 

9.1.1 Intra-industry Trade Indices 

The question of intra-industry trade measurement and the inherent problems was 
first explicitly analysed in early works of Grubel and Lloyd (1971, 1975). 
Nevertheless, some previous papers (Verdoon, 1960; Michaely, 1962; Kojima, 
1964; Balassa, 1966) already provided several procedures which could be applied 
to the measurement of intra-industry trade; although, that was not the main 
purpose of these works. 

Verdoom (1960) used the following indicator to examine changes in trade 
pattems in the Benelux: 

X i 

s i =— 

J M j 

where Xj and Mj are exports and imports of commodity j and Sj takes values inside 
the interval [ 0 , +oo ). When the value of the index comes closer to unity along 
time, it would be an indicator of increasing intra-industry specialisation. 

The disadvantage of this measure, as Gmbel and Lloyd (1975) already pointed 
out, is that any fraction and its inverse measure the same level of intra-industry 
trade, what makes it complicated to compare among industries. It can be seen that, 
except for the case where Sj = 1, this indicator does not measure directly the 
degree in which imports and exports overlap in a specific group of products. 

Kojima (1964) solved this problem by defining the degree of “horizontal 
trade” in a group of products j between two countries A and B as follows: 

X i 

D.- =— if Mi > X; 

J Mj 



and 
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iVl; 

D, = -*- if M, < X: 
1 *i 



Dj takes values inside the interval [0, 1], making it easier comparison among 
industries. The value of the index approaches unity as the degree of “horizontal 
trade” increases . Although the index D, represents an improvement compared to 
the Verdoom (1960) measure, both indices share a problem: the use of ratios 
between trade flows is not giving a direct measure of intra-industry trade as a 
Proportion of total trade. 

Kojima (1964) also provided an aggregated indicator of “horizontal trade” 
between two countries. The aggregated index was a weighted sum of indices Dj 
using as weight the share of trade of commodity group j in total trade. This 
aggregated indicator was defined as follows: 



n 

D=I 

Hl 



Xj Xj +Mj 



M • n 

J Z(Xj + Mj) 

V j =1 



= ZDj • Wj , if Mj > Xj 

j=l 



n 

D=X 

Hl 



Mj Xj+Mj 



X • n 

J Z(Xj+Mj) 
v J =1 



: X D j • w j , if Mj < Xj 
j=l 



where n is the number of groups considered and w j = 



Xj+Mj 



I(Xj+Mj) 

j=l 



2 



In a work aimed at analysing the effects of the completion of the Common Market 
on international specialisation in EC countries, Balassa (1966) proposed the 
following indices to measure the degree of overlapping between exports and 
imports: 



Note that Dj is equal to Sj when Mj > Xj and it is equal to Sj _1 when Mj < Xj . 
Note that it is in fact a weighted mean since the sum of weights is unity. 

( \ 



n n 

IWj = X 

j=l Hl 



Xj + Mj 



= 1 . 



KXj+Mj) 

VJ =1 

For this reason and from now on, we will call mean the type of weighted sums in which the 
sum of weights is 1 . 
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Xj-Mj 

A j = for commodity group y, 

X i +M j 



and 

— 1 n 

A ; = — X A i , as an aggregate index for n groups of products. 
n j=l 



where Aj is an arithmetic mean of indices Aj for a specific level of 

disaggregation. Therefore, this procedure of aggregation has the following 
shortcoming: it gives the same weight (1/n) to all industries independent of their 
share in total trade. 

The Aj index is inversely related to the degree of intra-industry trade. It takes 
values from 1 (when all trade is inter-industrial) to 0 (when all trade is intra- 
industrial). It has the advantage compared to previous measures of giving 
information on the proportion of total trade that is of the intra-industry type. 

Another important point to mention is that, as a result of expressing net trade as 
a proportion of total trade in a specific group of products, absolute values of 
imports and exports that may be quite different may produce the same value of the 
index . 

Moreover, suppose that commodity group j is formed by a number of 
subgroups m (which will be denoted by the subindex i). In this case, Aj may be 
expressed as follows: 

m 

KXij-Mjj) 

Aj =— 

J m 

HXij+Mij) 

i=l 



Assuming that either (Xy-My) <0,Vi=l,2,..,m,or (Xy-My) >0,Vi=l,2,..,m, 
then: 



m 



Z(Xjj -Mjj) 



2|Xij - M ij 



and we can write: 



A j = 

J m 



m . 

ZXj. 

i=l ' 



-M: 



KXij+Mjj) 



i=l 



m Xjj - Mjj 
ä(Xj+Mj) 
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Multiplying and dividing by (Xy + My): 



m Xjj-Mjj (Xjj+Mjj) 
JlCXjj+Mij)’ (Xj+Mj) 



m 



lAjj-Wy 



Xij-Mij 

where = is the correspondent index for subgroup i and where 

(Xjj+My) 



W; 



(Xjj+Mjj) 

(Xj+Mj) 



(Xij+Mjj) 



m 

KXij+Mij) 

i=l 



is the weight for that subgroup i. 



Therefore, Aj is a weighted average of the indices of the subgroups that form a 
group j only when the sign of trade imbalance is the same for all the subgroups. 
This weighting effect is lost, nevertheless, when there are trade imbalances of 
opposite sign. 

In general, a weighted average may be more suitable when we want to obtain a 
compendium measure, which reflects the relative importance of intra-industry 
trade at a particular level of aggregation . This is so specifically when such a 
compendium index refers to an economy in a wide sense, including all trade. On 
the other hand, if the level of aggregation that has been selected to calculate 
individual indices closely corresponds to the researcher’s view of homogeneity 
inside a group or industry, weighting individual indices may be inappropriate. 

Greenaway and Milner (1986) propose the following procedure to guarantee 
the weighting effect of the index: 



Ai 



m | i 

IXij-MJ 

i=l 

Xj + Mj 



Note that multiplying and dividing by (Xy + My) we have: 



m Xjj-Mij (Xjj + Mjj ) 
JiCXjj+Mij)' (Xj+Mj) 



m 

ZAjj -wjj 

i=l 



where w^ 



(Xjj+Mjj) 
(Xj+Mj) ' 



Note, for instance, that Aj is always 0 when Xj equals Mj , independent of the actual value of 
imports and exports. 

For instance, if we want to show the relevance of intra-industry trade in subgroups i, we 
would want Aj to be a weighted average of indices Ay. 
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Therefore, when all the subgroups trade imbalances have the same sign, the 
relationship Aj = Aj holds. 

Using the same procedure it is possible to obtain global weighted averages of 
all indices Aj, that is: 

— . i=i' 1 n 

A J=-i =SA j -w j 

Z(Xj-Mj) J=1 

H 



where wj = 



(Xj+Mj) 

n n 

ZXj + ZMj 
j=l j=l 



That is to say, the weight employed is the ratio volume of trade in industry j to 
total trade. 

Grubel and Lloyd (1975) proposed an index according to which the share of 
intra-industry trade in total trade flows of an industry j can be computed by the 
following expression: 

B (X j +M j )-| x j - M J -| t |Xj ~Mj| 

J (Xj+Mj) (Xj+Mj) 



where 0<Bj < 1. 

The average of indices Bj for a set of industries (j=T, 2, 3, ..., n) may be 
computed as follows : 



B; 



IWj -Bj =s — 

;_i J J :_1 n 



n (Xj+Mj) (Xj+Mj)- Xj-Mj 



j=l 



j=l 



I(Xj+Mj) 

j=i 



(Xj+Mj) 



S(Xj+Mj)-s|Xj-Mj| i|Xj-Mj| 
j=l j=l , j=l 



Z(Xj+Mj) 

H 



Z(Xj+Mj) 

j=l 



where 

j=l 






n 



Z(Xj-Mj) 

j=l 



when the sign of trade imbalances is not the 



same for all industries. 



Note that the same weighting procedure as in index Aj is used. 
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The advantage of index Bj when compared to Balassa’s A j is that the latter gives 

the same weight to all industries, whereas the former takes into account the share 
of industry j trade in total trade. 

A potential shortcoming of the aggregated G-L index B j is that it introduces a 

downward bias in the measurement of intra-industry trade when total trade is not 
n n 

balanced, that is when £Xj * £Mj . In this case imports and exports cannot 
j=i J =1 

match exactly in all industries and therefore the index can never reach its 
maximum value 1 . Grubel and Lloyd propose a correction procedure to deal with 
this problem, which we will discuss in section 9.2, totally devoted to trade 
imbalance adjustments. 

On the other hand, it can be noted that index B; is a modification of Balassa’s 
Aj , sharing the same properties conceming the weighting effect. However Bj is 
directly related to the level of intra-industry trade. It reaches its maximum value 
when total trade is of the intra-industry type and its value is zero when there is no 
matching at all between exports and imports in a specific industry. 

The same comments we have pointed out for the Balassa (1966) index 
conceming the trade imbalances opposite sign effects hold for the Bj index. If we 
consider the industry j composed of m subgroups /, we can rewrite the Bj index as 
follows: 



m / \ 


m / \ 


1 

+ 

w: 


X (Xjj - Mjj t 


1=1 i 


i=l 



ifXij+Mjj) 

i=l 



There is, however, an alternative aggregating procedure, which will provide the 
following index: 



B = Z Wy -Bjj 
J i=l 



ifXij+MijJ-lIXy-Mj 

i=l i=l 



ifXij+My) 

i=l 



(Xij+Mjj) (Xjj +My) Xy -My 

where w« = = and B,. =1 

'J m >J (Xjj+Mjj) 



As XlXij -Myl >[ 



i=l' 



IX +M i m 

( ^ KXy+Mjj) 

i=l 

y (x ,j - \1 , j | , we have that Bj > Bj* 
|i=l 



6 
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As we commented before, the question to consider is which procedure is more 
suitable. If the trade flows to aggregate are not homogeneous from the intra- 
industry trade point of view (that is, if they are not belonging to the same industry) 
it doesn’t seem desirable trade imbalance to cancel among themselves. In this 
case, Bj would be a more convenient measure to use. On the other hand, if trade 
flows to aggregate are homogeneous, the use of Bj would be preferable. 
Nevertheless, the problem in practice consists of finding out if a group of products 
included under a concrete level of aggregation by a Classification System can be 
considered as homogeneous from the perspective of intra-industry trade. That is to 
say, if this commodity group constitutes an industry . International trade flows 
Classification Systems may include under the same group products with very 
different factor requirements. This gives rises to the denominated categorical 
aggregation problem, which will be treated in section 9.3. 

Moreover, as Bj measures the degree of trade overlapping in a commodity 
group or industry j related to total trade in the industry, the behaviour of the index 
is affected by changes in total trade among industries and/or over time. For 
example, a high degree of trade overlapping may be registered in a specific 
industry in which the volume of trade is relatively small compared to total trade in 
manufactures. Besides, the volume of trade in an industry may or may not vary 
directly with the level of production or the level of sales. 

Other^ features of index Bj are its simetry regarding Xj and Mj and its non- 
linearity . 

A different way of addressing the question of intra-industry trade measurement 
is provided by Vona (1991). This author takes as elementary units the 5 digit SITC 
items and points out that, for this high level of disaggregation, all trade is to be 
considered as intra-industry trade or inter-industry trade depending on the 
existence of two way flows . Vona proposes the following intra-industry trade 
index for bilateral exports between two countries A and B in a commodity group i: 

Ia, b, i = X A , b, i + X B , a, i » if X A , B , i 0 and X B , a, i 0 

Ia, B, i ~ 0 ? if X a, B, i ~ 0 Or X ß, A, i ~ 0 
where 

X a, b, i = exports from country A to country B in commodity group i (5 digit 

SITC item) 

X b, a, i* = exports from country B to country A in commodity group i (5 digit 

SITC item) 



For an illustration of these features of Bj see Carrera (1996). 

Standard International Trade Classification. 

According to this approach industries may be classified under two categories. A first group of 
industries is characterised by scale economies, product differentiation and imperfect 
competition, giving rise to intra-industry trade. A second group is formed by industries under 
perfect competition producing homogeneous goods (according to the H-O-S model), which 
leads to inter-industry trade. 
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Vona calculates an aggregated index for 3 digit SITC groups j from 5 digit SITC 
items i as follows: 



Ia, b, j 



n 



i=l 



X A,B,j +X B,A,j 



•100 



where 

X a, b, j = exports from country A to country B in commodity group j (3 digit 
SITC) 

X b, a, j.= exports from country B to country A in commodity group j (3 digit 
SITC) 

Vona’s indicator takes values in the interval [0, 100]. It is 0 when all Ia, b, i = 0 and 
there is no intra-industry trade and it is 100 when all I A , b, i 0 and all trade is 
intra-industry trade. 

As a feature of this indicator Vona mentions that it is not affected by the 
existence of trade imbalances, as all trade (imports and exports) is considered 
intra-industry trade or inter-industry^trade independent of the exact coincidence of 
the volumes of imports and exports . Nevertheless, the problem of sensitivity of 
the index to the number of items included in every industry still remains unsolved. 
That is, the intra-industry index tends to decrease as the number of subgroups 
included in an industry increases. 

9.1.2 Intra-industry Specialisation Indices 

Michaely (1962) proposed the following indicator to measure the degree of 
overlapping between the share of exports and the share of imports in a commodity 
group y: 



n n 

ZXj IMj 

j=l j=i 

The indices ranges from 0 to l 11 , higher values suggesting a higher degree of 
similarity in the composition of exports and imports. A value of 0 would imply 



— 1 n 

H = 1 X 

2 j=l 



According to this point of view it is not the degree of overlapping between imports and 
exports what matters, but the mere existence of exchange of products in both directions. 

n n 

Note that it is feasible for H to reach unity without global trade balance ( X X j = X M j ) 

j=l j=l 

or without trade balance in industry j (Xj = Mj). Actually H =1 when the following condition 
is satisfied: 
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perfect inter-industry trade, with non-existence of simultaneous exports and 
imports in^a specific commodity group for a country during the time period 
considered . 

Balassa (1966) introduces the concept of “revealed” comparative advantage and 
proposes the following indicator to measure it: 



x i 

RCA = — I 

Xoi 



X 

x 0 



where — — indicates the share of a country in world exports of commodity j and 
x gi 

X 

is the share of a country in world exports of manufactured products. 

x g 

When the value of the index tends to 0, it would be indicating a “revealed” 
comparative disadvantage for industry y, while a value above 1 would identify a 
comparative advantage in the correspondent industry . From a general point of 
view, this index is measuring similarity in trade^ pattems and thus it can be 
employed to measure intra-industry specialisation . As we will show later, the 
index is, in fact, equivalent to the measure proposed by Glejser, Goossens and 
Vanden Eede (1982) aimed at measuring intra-industry specialisation in foreign 
trade. 

Finger and Kreinin (1979) develop an indicator that uses the share of every 
industry’ s exports in a country’ s total exports to measure the degree of similarity 
of export pattems from two countries (a and b) to a third market (c). The indicator 
was computed as follows: 



X; 






12 



13 



14 



n n 

IXj 
j=l 



ZMj 

j=l 



The index originally proposed by Michaely was actually the following: 



n 

D=I 

j=ii 



Xi 



M ; 



n n 

IXj 
j=l 



IMj 

j=l 



This indicator ranges from 0 to 2. Index H (which ranges from 0 to 1) is just an adaptation 
of D in order to make it easier comparison to other measures. 

A value of RCA = 1.1, for example, would indicate that the share of a country in exports of 
commodity j is 10% higher than its share in total exports of manufactures. 

Aquino (1978) is using the Standard deviation of Balassas’s index as a measure of the 
intensity of inter-industry specialisation in a country. A deviation of 0 would indicate the 
absence of inter-industry specialisation while higher dispersion of the values of the index 
would point at higher inter-industry specialisation. 




9 The Measurement of Intra-industry Trade and Specialisation: a Review 



133 



S (ab, c) = £ min [Xj (ac), Xj (bc)], 
j=i 



where Xj (ac) is the share of industry j in the exports from country a to country c 
and Xj (bc) is the share of industry j in the exports from country b to country c. 
The indicator takes the same values than the index of Michaely (1962) and, in fact, 
both measures are equivalent, as it is shown in Kol (1988). For this reason, 
Michaely’ s index can be used to compare trade pattems from many points of view. 
It can be employed, for instance, to compare imports and exports pattems in a 
country, exports (or imports) pattems for two countries or a group of countries, 
among them or related to a third market. 

Glejser, Goossens and Vanden Eede (1979) mark a break with previous 
methodology and propose a conceptually different approach to quantify the 
magnitude and variations of intra-industry trade. Unlike other measures that 
consider simultaneously imports and exports, these authors make a distinction 
between supply specialisation (exports) and demand specialisation (imports). The 
indicators are built on the assumption that a country is specialised in a specific 
industry when it exports (or imports) relatively more than a group of conveniently 
selected countries. 

The comparative supply specialisation for a given time period can be computed 
as follows: 



1 n 

S=- I log 

n j= i 






X 



ß ) 



1 n 

= 1 1 Sj 

n j=l 



where: 



n = number of industries considered 

Xj = exports from industry j made from a country to the group of countries 
considered 

X = total exports from a country to the group of countries considered 

X a = total exports from industry j made by the group of countries considered 

(except for the country that is being analysed) 

X g = total exports from the group of countries considered (except for the 
country that is being analysed) 

The index of comparative demand specialisation is similarly computed as follows: 



1 n 

M = - S log 
n j=1 



Mi M 



/ 



gj 



M M 



ß ) 



1 £ 

= " 2 Hj 

n j=i 



where: 

n = number of industries considered 

Mj = imports from industry j made from a country to the group of countries 
considered 

M = total imports from a country to the group of countries considered 
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Mgj = total imports from industry j made by the group of countries considered 
(except for the country that is being analysed) 

M g = total imports from the group of countries considered (except for the 
country that is being analysed) 

A higher divergence between X/X and Xgj/X g and between M/M and Mgj/M g is to 
be expected with a higher degree of inter-industry specialisation. On the other 
hand, if intra-industry specialisation is dominating, the quotiens (Xj/X)/(Xgj/X g ) 
and (Mj/M)/(Mgj/M g ) approach unity in every industry j and, consequently, the 
1 n 1 n 

unweighted means £= — £ and ju = — £ Hi approach zero. 
n j=i J n j= i 

The variability among industries, that is, the variances of £, and p can be 
computed as follows: 



n j=i 



n j=i 

Kol and Mennes (1986) and Kol (1988) establish an equivalence between this 
indicator and the one introduced by Balassa (1966) to measure “revealed” 
comparative advantage based on the share of exports from every industry in total 
world exports . Although the purpose of those authors was not the same, both 
indicators are able to measure “revealed” comparative advantage and intra- 
industry specialisation in the form the authors define such concepts. 

9.1.3 Comparison of Measures 

In previous subsections we have examined two groups or families of measures and 
evaluated the most commonly used indices of intra-industry trade and 
specialisation. The first group of indicators is mainly^ based on the degree of 
overlapping of imports and exports in every industry . The second group of 
measures takes into account the degree of similarity in relative structure of imports 
and exports. 

Conceming the first family of measures, the following considerations can be 
mentioned: 

- The index Sj introduced by Verdoom (1960) has the following shortcomings: 
a) Every ffaction and its inverse measure the same degree of intra-industry 



Operating in RCA we obtain the following relationship with the index of Glejser, Goossens 
and Vanden Eede: 

= log (RCA). 

16 Except for the index of Vona (1991). 
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trade, making it complicated to compare among industries. 

b) The index is not giving a direct measure of intra-industry trade as a 
Proportion of total trade. 

- The index D, introduced by Kojima (1964) is solving the former problem in a) 
but the problem in b) remains unsolved. 

- The measure Aj introduced by Balassa (1966) has, compared to the previous 
ones, the advantage that it gives information on the proportion of total trade 
that is of the intra-industry type. However, it has the unattractive feature of 
being inversely related to the degree of intra-industry trade. Moreover, the 
aggregated measure Aj is giving the same weight to all industries, independent 
on their share in total trade. 

- Grubel and Lloyd (1975) introduce the index Bj, which is the most widely used 
in the empirical measurement of intra-industry trade. The Bj index shares some 
features of the Aj index of Balassa (1966) but has the advantage of being 
directly related to the level of intra-industry trade. Moreover, the aggregated 
measure Bj is preferable to Aj because it takes into account the share of every 

industry in total trade unlike the later, which gives the same weight to every 
industry. 

- The index I Ai B , i of Vona (1991) gives a different approach to the problem of 
intra-industry trade measurement. It has the advantage of not being influenced 
by trade imbalance. However, its use at the empirical level is complicated 
because the information on the high level of data disaggregation required is 
not easily available and it is difficult to handle. 

From the previous analysis it can be concluded that the indicators proposed by 
Grubel and Lloyd (1975) are the more convenient among those indices that form 
the first group of measures. 

Conceming the second family of measures a distinction can be made among those 
that are using differences in relative shares of imports and exports in total 
respective flows (Michaely, 1962; Finger-Kreinin, 1979) and those, which are 
using quotients among such shares (Balassa, 1966; Glejser, Goossens y Vanden 
Eede, 1979). Following Kol (1988) and with the purpose of making it easier to 
compare both groups of measures we will first express the index of Balassa (1966) 
in terms of the indicator of Michaely (1962): 



RCA = 



2Ei 

X M 



This indicator has some shortcomings. First, it has to be corrected in case M, = 0, a 
case not uncommon at high levels of data disaggregation. Second, while a total 
similarity of trade pattems gives a value zero for the index, in case of increasing 
dissimilarity the value of the indicator tends to zero or oo. This unlimited ränge of 
values complicates the interpretation of results conceming the degree of 
dissimilarity of trade pattems. A third problem is the following: when changing 
numerator and denominator we obtain completely different values of the index 



17 Vona proposes the 5-digit SITC level of aggregation. 
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(and the variance of the sample). The degree of similarity is however the same . 
Finally, the use of quotients to compare trade pattems gives rise to a fourth 
problem. Imports and exports with a comparatively small value but very different 
shares in total respective flows may influence the value of the index in an 
undesirable way, by increasing it in case of high similarity of exports and imports 
pattems . 

On the other hand, the index provided by Michaely (1962) has none of these 
disadvantages. Values of Xj = 0 or M, = 0 are not a problem. The interpretation of 
results is easier because the index ränge from 0 (complete dissimilarity) to 1 
(complete similarity). The third and fourth problems do not apply here because 
differences are used rather than quotients. Because of all these reasons it can be 
concluded that, when trying to measure the degree of similarity of trade pattems, it 
is preferable to use Michaely’ s index. As stated before, this indicator has a wide 
ränge of applications in numerous situations and it is suitable to compare trade 
pattems from diverse points of view. 

Two criteria have been thus basically used when it comes to measuring intra- 
industry trade ’s intensity: the magnitude of trade flows and the similarity in 
relative stmcture of imports and exports. These two criteria give rise to two 
families of indices, which are not strictly comparable. However, Silber and Broll 
(1990) show that both types of measures can be altematively expressed as 
measures of distance or measures of similarity, deriving the corresponding family 
of indices. The empirical section of their work shows that, although the two 
groups of indices measure different types of distances or equalities, both are 
highly correlated. 

9.1.4 Later Developments 

In later years a great deal of work on IIT measurement has concentrated in two 
refinements of the traditional Grubel-Lloyd index. One of them is aimed at 
measuring what has been called marginal IIT, an issue related to the analysis of 
adjustment costs associated to trade liberalisation. The other one is targeted at 
adjusting IIT measures to disentangle vertical and horizontal IIT. 

9. 1.4. 1 Measuring Marginal Intra-industry Trade 

It is generally assumed that the costs of the adjustment process following trade 
liberalisation may differ depending on whether the new trade generated is 
classified as inter-industry or intra-industry trade. If emerging trade has an intra- 
industry nature the adjustment costs will be probably lower, because in this case 
reallocation is to take place within every industry rather than among different 

18 

Using the logarithm of the indicator may solve this third problem. This is, in fact, what 
Glejser, Goossens and Vanden Eede are doing. In this case, when changing numerator and 
denominator, the indicator changes its sign but not its value. However, the aggregated index 
still ranges from 0 to . The second problem remains thus unsolved. 

19 

Kol (1988) is giving an illustrative example of this feature. 

20 

That is, those that are based on the degree of overlapping of trade flows and those based on 
the similarity of trade pattems. 
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industries. When analysing the adjustment consequences of trade expansion, what 
is relevant is not the current level of IIT but the composition of changes in exports 
and imports. In other words, how IIT changes at the margin. 

To capture this important issue, Hamilton and Kniest (1991) developed a 
measure, which calculated the share of IIT in new trade flows as follows: 



MUT = 



Xt Xt ~ n for M t - M t _ n > X t - X t _ n > 0 



M t -M t _„ 



x t -x t _ n 



for X t - X t _ n > M t - M t _ n > 0 



where Xj and Mj refer to exports and imports of commodity group j and t and t-n 
refer to the two points in time taken into consideration. The index is calculating 
the proportion of the increase in imports or exports that is matched. When all new 
trade is matched (the increase in imports totally matches the increase in exports 
for a particular industry) the index will be unity. When all new trade is inter- 
industry trade (there is no matching at all between new imports and new exports 
for a concrete industry) the index will be zero. 

The Hamilton and Kniest measure, although being a considerable improvement 
for the analysis of adjustment issues, has some important shortcomings, as several 
authors have pointed out [see, for example, Brülhart (1994) and Greenaway, Hine, 
Milner and Elliot (1994)]. First of all, the index cannot be calculated when the 
change in exports or imports is negative. Second, the measure is unscaled, that is 
to say, it is not related to total amount of new trade or to the initial level of trade or 
to the value of production in the industry considered. Third, the index measures 
the changes in nominal terms rather than in real ones. To address these 
shortcomings the aforementioned authors have proposed several alternatives 
providing a menu of indices (which deal with the scaling problem and are always 
defined) to choice depending on the purpose of the job at hand. 

The Brülhart (1994) index has been widely used in later works on adjustment 
issues. This dynamic measure can be expressed as follows: 

A t |AX i -AM i | _ i l(X t -X t _ n )-(M t -M t _ n )| 

’ lAXil + lAMfl |X t -X t _ n | + |M t -M t _ n ’ 



where Xj and Mj refer to exports and imports of commodity i and t and t-n refer to 
the two points in time taken into consideration. The index has values ranging ffom 
0 to 1, a value of 0 indicating that new trade is entirely of the inter-industry type 
and a value of 1 showing complete matching between new exports and imports in 
industry j and so lower transitional adjustment costs. The measure is dynamic in 
the sense that it provides information on the proportion of changes in total trade 
flows that are of an intra-industry type. 

Recent work by Thom and McDowell (1999) shows that the proposed marginal 
intra-industry trade indicators [such as Brülhart (1994)] may be underestimating 
the extent of intra-industry trade as they cannot distinguish between inter-industry 
trade and vertical intra-industry trade. These authors propose an alternative 
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procedure of classifying marginal trade flows based on the joint application of the 
Brülhart’s index and the aggregate measure defined as follows: 

Iax^amJ 

Aj = 1 i ii— 

llAXil+SlAMil 

i=l i=l 



where Aj measures the extent of total intra-industry trade in industry y, which has n 
subindustries, aggregating vertical and horizontal intra-industry trade. The 
difference between Aj and Brülhart’s index would give vertical trade and inter- 
industry trade would be given by the residual. 

9 . 1.4.2 Disentangling Vertical and Horizontal Intra-industry Trade 



The distinction between horizontal and vertical differentiation of product is an 
important one when dealing with adjustment issues. As several authors point out 
[see Greenaway, Hine and Milner (1994, 1995)] horizontal IIT is likely to lead to 
lower adjustment pressures than vertical IIT. This is so since different industry 
and country characteristics can be associated with the exchange of products 
depending on whether such products involve similar or different qualities 
(horizontal or vertical differentiation). 

Greenaway, Hine and Milner (1994, 1995) propose a methodology to identify 
vertical and horizontal IIT built upon previous work of Abd-el-Rahman (1991). 
This approach assumes that quality is reflected in price and price can be proxied 
by unit values. The share of vertical and horizontal IIT in total trade is computed 
as follows: 



mj 1 = 1 





(9.1) 



HTj =1 



I 

i=l 



Xy-My 



S( x ü +M ü 



i=l 



I XX+Mj 



i=l 



x J+Mj 



(9.2) 



where i refers to the 5 th digit SITC products in a commodity group or industry j. 
Total IIT is then decomposed in vertical IIT and horizontal IIT as follows: 



IIT* = IIT- V + IIT- h 

j J J 
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IIT* = 1 - 



i ( x e + M„ ) 



1* 

Horizontal intra-industry trade (IITj ) is given by (9.1) for the items i in 
commodity group j where the following condition holds: 



1 -a 



<u\^ 

uv™ 



< 1 + OC , 



where UV X and UV m refers to unit values of exports and imports and a is a given 
dispersion factor. 

Vertical intra-industry trade ( IIT^ ) is given by (9.2) for the items i in 
commodity group j where the following condition holds: 



u yj 

uv™ 



< l-a 



or 



uv* 



uv™ 



> 1 + a . 



The dispersion factor a has been given several values in the diverse works where 

this procedure is employed giving different price wedges within which IIT is 

considered as horizontal and outside of which it is defined as vertical. However, 

even when large price wedges have been used, vertical IIT has been found to be 

very significant in empirical studies. The results obtained suggest that the 

distinction between horizontal and vertical IIT should be taken into account in 

22 

econometric models, as the determinants of both types of IIT may differ . 



9.2 

Trade Imbalance Adjustment 



9.2.1 Grubel and Lloyd (1975) 

One of the shortcomings of the aggregated index ß , is that, in case of global trade 

n n 

imbalance (that is, ZXj £Mj ), it can introduce a downward bias in intra- 
j=l j=l 



Abd-el-Rahman (1991) uses a = 0.15 and Greenaway, Hine and Milner (1994, 1995) employ 
both a = 0.15 and a = 0.25. 

22 

Some work on this subject can be found in Menon, Greenaway and Milner (1999). 
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industry trade measurement. This is due to the fact that, under this assumption, 

n j * 

IXj-M j > 0 and, therefore, B j can never reach its maximum value unity. 
j=i 

To deal with this problem, Grubel and Lloyd (1975) proposed an alternative 
adjusted measure, which expresses intra-industry trade as a proportion of total 
trade minus the value of trade imbalance (in absolute terms), as follows: 



i(X j +M j )-l|x j -M j | _ 
p = m = __i 





n 


1-K 


KXj+Mj)- 


KXj-Mj) 




j=l 


j=l 





where K = 



Z(Xj-Mj) 
.H 

Z(Xj+Mj) 

j=l 



— 23 

and where 0 < c j<i . 

As it can be observed, Cj increases as the value of K (value of trade imbalance as 

a proportion of total trade) gets higher. On the other hand, the index refers to 
aggregated trade flows and it has no counterpart for individual industries. 

Moreover, when trade imbalance has the same sign in all industries, C j reaches 

24 

the value 1 , independent of the size of such trade imbalances . 

In short, the adjusted measure Cj intends to indicate what the level of intra- 

industry trade would be when no trade imbalance would exist. Nevertheless, 
Greenaway and Milner (1986), among others, argue that the need for the index to 
be adjusted when that is not the case should be based on wider theoretical grounds 
than the mentioned fimctional restriction of Bj . Specifically, it has to be 



23 



24 



Note that X(Xj + Mj) > £ Xj — Mj 



H 



Note that, in this case, 






SCXj-Mj) 
j=l 



Z(Xj-Mj) 

j=l 



n | 

= X X j 

j=l 



■M; 



The index reaches thus its maximum value unity even when not all trade is intra-industry 
trade in the analysed country. Vona (1991) points out that this problem is not, nevertheless, 
very important at the empirical level. The probability for this problem to appear is inversely 

related to the number of elementary items used to compute Bj or C j , what makes it 

difficult to happen when using a high level of data disaggregation. 
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considered that a Situation of global macroeconomic equilibrium may be 
compatible with the existence of trade imbalance for a group of industries. 
Moreover, in case we actually have a Situation of disequilibrium, the direction of 
the adjustment may be different from the one proposed by Grubel and Lloyd. The 
restoring forces could imply an increase in trade balance and, in this case, Bj 
would be upwards biased. 

Finally, it is worthy to mention that Grubel and Lloyd’s argumentation on the 
need to correct in the presence of global trade imbalance implies the Classification 
of trade flows in three categories: inter-industry trade flows, intra-industry trade 
flows and trade imbalances. The third of these categories is considered as a 
“disturbing” factor that has to be excluded of the analysis. Nevertheless, and if we 
take into account that most empirical works on intra-industry trade refer to the 
manufacturing sector, it is difficult to find a theoretical reason why trade in 
manufactured products should be balanced. A country may have, by instance, a 
deficit in this sector while having a superavit in other sectors. This matter becomes 
more evident if we consider bilateral trade flows, because it is not clear which 
bilateral position would be consistent with a Situation of multilateral 
equilibrium. 25 . 

The critics to the Grubel and Lloyd’s adjustment refer not only to the suitability of 
the proposed procedure [Aquino (1978)], but also to the very need to adjust global 
trade imbalance [Greenaway and Milner (1981, 1983), Kol and Mennes, 1983). 
Other contributions to the debate may be found in Greenaway (1984), Pomfret 
(1985), Greenaway and Milner (1986, 1987), Kol (1988), Kol and Mennes (1989) 
and Vona (1990, 1991). 

9.2.2 Aquino (1978) 

Aquino (1978) argues that, if Bj then there is a downwards biased when there 

exists a trade imbalance. It is precisely because individual indices Bj are 
themselves downward biased measures of intra-industry trade at the industry level. 
Consequently, this author proposes to adjust individual indices Bj under the 
assumption of equiproportional distribution of trade imbalance effects among all 
industries. The adjustment procedure begins by estimating the value of exports 
and imports in case no global trade imbalance would exist. For this purpose, the 
following expressions are computed: 

i(x j+Mj ) 

X je =X J ---^ (9.3) 

^ -X 7 



The equilibrium condition could be thus imposed only when measuring intra-industry trade 
including total exchange of goods and Services between a country and the rest of the world 
(excluding, to simplify, Capital flows). 
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i(Xj+Mj) 

M i e = M r~ tL - n ( 9 - 4 ) 

IM- 

j=i 



Operating in (9.3) and (9.4) we have: 

£x je =IM je =2. X(Xj +Mj) 

j=l j=l 1 j=l 



Using these “artificial” values the following intra-industry trade index for 
commodity group j can be computed: 



X ; 



Mi 



Qj 



X je M je 
(X je +M je ) 



n n 

IZXj ZMj| 

lj=l j=l 



J e ' 



X; 



Mi 



n n 

SXj SMj 

>1 j=l 



Finally, based on individual indices Qj the following aggregated index is 
computed: 



S(X j+ Mj)-i|x je -M je | 

Qj=— (9-5) 

Z(Xj+Mj) 

j=l 



Examining (9.5) it can be deduced that this indicator uses as a criterion for 
measurement the degree of similarity of relative shares ojfimports and exports. 
Taking into account this feature, an equivalence between Qj and the Michaely’s 
indicator H can be obtained as follows: 
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I(Xj +Mj)-T I(Xj + Mj) - 1 
j=l 1 j =I j =1 


Xj Mj 


n n 

IXj IMj 




j=l j=l 



ICXj+Mj) 

j=i 




i 

j=i 



X 



n 

ZXj 

j=i 



Mi 



n 

IMj 

j=l 



= H 



The use of this adjustment procedure does not imply a Variation in total trade 
flows as the following holds: 

i(Xj +Mj)= £(X je +M je ) (9.6) 

j=l j=i 



Exports and imports pattems are also preserved after adjustment, that is: 



X je Xj 



n 


n 


SX je 


IXj 


j=i 


j=i 


and 




M je 


_ M J 


n 


n 


IM je 


IMj 


j=l 


j=i 



The above equality (9.6), however, does not hold necessarily at the industry j 
level. In fact, when total trade is not balanced, the fulfilment of (X je +Mj e = Xj+Mj) 
implies that: 



n n 

IMj IXj 

j=l j=l 

On the other hand, as Greenaway and Milner (1986) point out, the direction of the 
adjustment of Bj for every industry depends on the relationship between the sign 
of trade imbalance in the industry and total trade imbalance. 
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Xj > Mj and Z X j > Z M : imply that Qj > Bj; 
j=i j=i 



Xj > Mj and Z X j < Z M j imply that Qj < Bj; 
j=l j=l 



Xj < Mj and Z X j > Z M j imply that Qj > Bj; 
j=l j=l 



Xj < Mj and Z X j < Z M j imply that Qj < Bj. 
j=i j=i 

This variability in the direction of the correction affects the ranking of industries 
according to their intra-industry trade index. The election of the group of 
transactions, which will be the basis for the adjustment of trade imbalance, 
becomes thus critical, especially when the adjusted index is to be used in 
econometric analysis. Other critics to this procedure are the following: first, the 
fact that a trade imbalance in a group of industries is considered as a 
disequilibrium Situation and, second, the assumption that the effects of restoring 
forces are equiproportionally distributed among industries . For these reasons, it 
would be advisable, when analysing the determinants of intra-industry trade, a 
sensible selection of the time periods included in the sample rather than the use of 
adjusted indices. Another alternative proposed by the aforementioned authors is 
the use of an average measure of the indices corresponding to a carefully selected 
period of time. 

9.2.3 Balassa (1979) 

Balassa (1979, 1986) utilises a correcting procedure similar to the one used by 
Aquino (1978). However, while the latter is using as a basis for adjustment the 
group of transactions analysed (manufactured products, total trade^or other sets of 
transactions), Balassa only takes into account total trade imbalance . 

Technically, the Balassa adjustment is identical to Aquino’s correction: 



26 

Aquino (1981) argues that such a feature can be acceptable in absence of other elements that 

help to identify an alternative criterion. 

27 

Balassa (1986) justifies this selection “in order to allow for inter-industry specialisation 
between primary and manufactured products”. 
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X b =X j — 
J 2 



S(Xj+Mj) 



J 



n 

IXj 

j=l 



and 



SCXj+Mj) 

=M; • — 

J J 1 n 

IMj 

H 

Like in the previous case, the correction is equiproportionally distributed among 
industries . However, the election of a specific trade imbalance as the basis for 
the adjustment represents an improvement regarding the problem of ambiguity 
commented for Aquino (1978). 

9.2.4 Loertscher and Wolter (1980) 

Loertscher and Wolter (1980) use the same correcting procedure as Aquino 
(1978). However, they apply it to bilateral trade flows of manufactured products 
instead of total trade of a country with the rest of the world. We have already 
mentioned, nevertheless, the lack of theoretical support for the existence of 
bilateral equilibrium with specific countries in a concrete set of transactions. 
Greenaway and Milner (1981) cornment that, particularly, bilateral disequilibrium 
may be the result of factors leading to inter-industrial and intra-industrial 
specialisation. It may be thus inconvenient to adjust on the basis of bilateral 
disequilibria because it can hide the influences that one pretends to measure. 

9.2.5 Bergstrand (1982) 

Finally, Bergstrand (1982) considers total trade imbalance as the basis for 
adjustment, as in Balassa (1979), but^opposite to the later, he uses bilateral trade 
flows to measure intra-industry trade . 

The indicator proposed for the measurement ^of bilateral intra-industry trade 
between countries i and j in products of industry k is the following: 



28 

In case of a deficit (superavit), all exports are increased (diminished) in the same proportion. 

The same applies for imports. 

29 

Bergstrand justifies this decision arguing that in a model with multiple countries, multiple 
goods, two factors of production and non equality of factor prices, the Heckscher-Ohlin 
theorem is always fulfilled for bilateral trade flows between pairs of countries and thus the 
focus of interest should be the presence of bilateral intra-industry trade, which is “non 
expected” according to the mentioned theorem. 

Adjusted on the basis of each country’ s multilateral trade imbalance. 
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IITy * = 1 - 



where 

X lc«_l (X; +M,) (Xj+Mj) ! xk 

ij 2 2Xj 2Mj ’J 

x k*_l f ( x j +M j) (Xj+Mj) ! xk 
ji 2 [ 2Xj 2Mj J' 1 

Xy = value of exports from country i to country j of products in industry k 

n m , 

Xj = X Z X = total exports of country i 
k=l j=l J 




n m , 

Mj = 2 XX” = total imports of country / 
k=l j=l 

Finally, X » , Xj and Mj are similarly defined. 



Taking, for example, the case of exports of industry k from country i to country j it 



1 (Xj+Mj) (Xj+Mj) 



can be noted that the correcting factor 



arithmetic mean of two elements. The first of them, (Xi +M t ) j ncrease 

2Xj 



exports from i when such country presents a global trade deficit (Xj+Mj > 2Xj) and 

(Xj+Mj) 

to reduced them in the opposite case. Regarding the second element, — ’ 

it tends to increase exports from i to j when country j has global trade superavit 
(Xj+Mj > 2Mj) and to reduce them in the opposite case. A similar analysis applies 

to X J . 

In short, both elements of the correcting factor act in the same direction when 
one of the countries has a deficit and the other one has a superavit. However, both 
elements tend to offset when the sign of trade imbalance is the same in the two 
countries. 
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On the other hand, it can be noted that the correcting factors for Xy and X ^ 

31 

are independent on the industry considered k . This means that the correction is 
equiproportionaly applied to all industries, the same as in Aquino’s and Balassa’s 
procedures. 

The correcting procedure stops when countries i and j present a multilateral 
equilibrium (independent on the existence of bilateral disequilibria). It is thus an 

iterative procedure to compute X~* and X^* in which bilateral trade flows are 

changed until all countries reach a multilateral equilibrium. The amount of 
information required to apply this method is much higher than the one needed by 
the other indicators commented. 

9.2.6 Comments and Conclusions 

We have revised in previous subsections different methods, which can be applied 
to correct for trade imbalance, either in global trade or in specific sectors, for 
bilateral trade flows or multilateral exchanges. 

The correction of the index Bj was introduced to allow the indicator to reach its 
maximum value unity, even in the presence of trade imbalance. However, this 
argument has no theoretical foundations. Moreover, the adjusted index Cj has no 

counterpart at the industry level, as opposite to indices Bj and Bj . Finally, we 
have to mention that, in case all trade imbalances have the same sign in all 
industries, the index Cj reaches its maximum value 1 independent on the size of 

such trade imbalances 32 . 

Conceming the Aquino (1978) index, although it is an improving compared to 
the previous one, it suffers from several shortcomings. In concrete, we have 
already commented that the ranking of industries according to their level of intra- 
industry trade experiments significant variations, depending on the trade balance 
of the set of transactions considered as a basis for the adjustment. On the other 
hand, the correction is equiproportionaly distributed among industries, what may 
imply the lost of important information on specific industry characteristics related 
to intra-industry trade. Moreover, the Aquino (1978) index always takes the same 
value when the share of imports and exports in every industry in total imports and 
exports, respectively, remains constant. 

In Aquino (1978) and Loertscher and Wolter (1980) the adjustment is made on 
the basis of trade imbalance in manufactured products but, again, there is a priori 



The exports of every industry k are multiplied by the same factor and the same applies for 

imports. 

32 

This is due to the following: 
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no theoretical foundations for trade in a specific group of transactions to be 
balanced. 

Loertscher and Wolter (1980) correct in terms of bilateral trade flows but the 
situations of bilateral disequilibrium may be consistent with multilateral 
equilibrium and, in fact, they may be reflecting a trend towards inter or intra- 
industry specialisation. For this reason, it is even more hazardous to use and 
adjustment procedure of the Aquino (1978) type on a bilateral basis because it 
may hide the very influences that one pretends to measure. 

In our opinion, a convenient method to correct for trade imbalance should 
satisfy, at least, the following conditions: first, it should be applicable to several 
levels of data aggregation; second, it has to use total trade as a basis for the 
correction (and not only a particular set of transactions) and third, it has to 
simulate a Situation of multilateral equilibrium. Taking all this into account, the 
methods proposed by Balassa (1979) and Bergstrand (1982) seem to be more 
suitable than the other procedures commented. Balassa’ s adjustment has, 
compared to Bergstrand, the advantage of its simplicity. 

Nevertheless, none of the exposed procedures considers the causes underlyin^ 
the disequilibrium situations considered, conveniently adjusting the procedure . 
In fact, these methods could create more distortions than the ones they pretend to 
eliminate and have no clear relationship with theoretical foundations. However, in 
some cases, particularly when it comes to comparing among countries with 
different trade imbalances or for a country with very different trade imbalances 
along the time period considered, the use of adjusted indices may be convenient. 
Some alternative methods may be used [Greenaway and Milner (1981)]. A 
sensible selection of the years considered could be made, which avoids periods of 
evident global trade imbalance. On the other hand, an average of intra-industry 
trade indices may be used for a carefully selected time period. 

In any case, the distortions that trade imbalance may introduce in the 
measurement of intra-industry trade are probably insignificant when compared 
with those brought about by the definition of “industry” and the consequent effects 
of categorical aggregation. We deal with this question in next section 9.3. 

9.3 

Categorical Aggregation 

9.3.1 Definition of the Problem 

One of the most important problems of intra-industry trade ’s measurement is the 
unknown influence of categorical aggregation. The identification and 
measurement of intra-industry trade clearly depends on the degree of homogeneity 
of the products classifled under the same group in trade statistics and the nature of 
such homogeneity. The bias introduced by categorical aggregation on intra- 
industry trade measurement emerges when essentially heterogeneous products 

33 



The equiproportional distribution of the correction among Industries reflects the lack of 
flexibility of the proposed procedures. 
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(that is, products than cannot be considered belonging to the same industry) are 
classified under the same group in trade statistics. In general, the existence of bi- 
directional trade flows for a commodity group will identify “genuine” intra- 
industry trade and Statistical aggregation (simultaneous exchange of products 
classified under the same category but with different factor requirements). 

The most convenient procedure to eliminate this problem seems to be, at first 
sight, the reclassification of data in a way that all resulting categories correspond 
as close as possible with the theoretical construction of an “industry”. The 
problem is that there are multiple criteria of regrouping the data and, even is a 
specific criterion is systematically used, reclassification is an extremely arduous 
job . 

International trade data may be classified according^to different criteria. The 
most commonly used System is the SITC nomenclature , in which products are 
grouped at various levels of aggregation identified by digits. The alternative to 
reclassification would be the election of a specific level of data disaggregation in 
official classifications as the best approximation to the concept of industry. In this 
case, a double question emerges. On the one hand, the most convenient level of 
data disaggregation has to be identified (that is, the level that more narrowly 
corresponds to the concept of “industry”). On the other hand, the influence of 
categorical aggregation for the selected level of data disaggregation has to be 
determined. 

9.3.2 Categorical Aggregation Tests 

Due to the difficulties of product reclassification under a systematic and coherent 
with the economic criterion procedure, several alternatives have been proposed. 
Particularly, three alternative methods to assess the influence of categorical 
aggregation on intra-industry trade indices have been suggested : 

i) To perform measurement at higher levels of Statistical disaggregation. 

ii) To execute diverse measurements based on different Classification Systems. 

iii) To compute an adjusted intra-industry trade index. 

The first of these procedures has been the most widely employed. It is based on 
the consideration that a substantial fall in average intra-industry indices when 
descending from a specific level of aggregation to a higher one may be indicating 
the presence of categorical aggregation problem. 

However, an important point has to be taken into account. There is no specific 
Standard pattem to evaluate if the fall in the value of the index is significant and so 
the results of this test are not at all conclusive but have an indicative character . 



An additional problem would be the Classification of pieces and components. 

^ Standard International Trade Classification. 

Greenaway and Milner (1983). 

37 

Usually, a fall in the value of the index of about 20% when passing from the 3-digit SITC to 
the 5-digit SITC level of aggregation, while maintaining a reasonably high volume of intra- 
industry trade for the 5-digit level, is considered to indicate that the main component of 
recorded intra-industry trade is not categorical aggregation. 
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Conceming the second type of categorical aggregation test mentioned, it is also 
instructive to observe the sensibility of indices Bj to Classification Systems based 
on different criteria (like processing characteristics or product features). The 
stability in the ranging of indices computed using different classifications may be 
illustrative but, like in the previous case, no conclusions can be obtained regarding 
the absolute significance of the categorical aggregation error. Moreover, the need 
to select a specific level of data aggregation for the indices’ computation still 
remains. 

The third type of categorical aggregation test involves a more systematic way 
of evaluating the categorical aggregation problem associated with trade 
imbalances of opposite sign in subgroups of products with different factor ratios 
and low degree of Substitution . If we considered commodity group j formed by 
m subgroups denoted by i, this procedure consists of computing an adjusted intra- 
industry trade index according to the following expression: 

m V ( X ij + M y ) - X |x ij M j j | 

B * = Zwy -Bjj & 

1=1 KXjj+Mij) 

i=l 



where Wy = 



_X, + M S Ft-“« 



m X • +M • 

KXij+Mjj) A J +M J 

i=l 



Xij+My 



The aggregation procedure to compute the intra-industry trade index for 
commodity group j consists of computing By indices for the subgroups i and 
averaging them according to the share of every subgroup in total trade of j. 



As X Xy-Mij £ 
i=l 



KXii-Mii) 



i=l 



holds, we have 1 > Bj > B j >0. 



This means that, the more trade imbalances of opposite sign cancel between 



subgroups, the higher Bj will be related to B j and, in case all trade imbalances in 
the subgroups have the same sign, we have Bj = B j . 

The adjustment is thus based on the assumption that categorical aggregation is 
associated with opposite sign trade imbalances in the subgroups i of which a 
commodity group j is composed. When this Situation exists, a measurement of 
intra-industry trade at a higher level of data disaggregation generates a lower value 
of intra-industry trade and helps to correct for the influence of Statistical 
aggregation. 
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9.3.3 Critics and Conclusions 

The available evidence collected in the literature devoted to examine the influence 
of Statistical aggregation on intra-industry trade measurement and, specifically, the 
existence of high volumes of intra-industry trade even at high levels of data 
disaggregation, shows that this is a phenomenon that cannot be explained solely 
by Statistical aggregation. 

There is no doubt that the interpretation of intra-industry trade measurement 
becomes more complicated with the existence of categorical aggregation. The 
more convenient way of facing this problem could be the reclassification of data 
building product categories as homogeneous as possible, (that is, according to 
what can be considered an “industry”). However, this regrouping is problematic 
due to the lack of a unique reclassification criterion, among other factors. 

The election of a concrete level of data disaggregation to measure intra-industry 
trade is not either exempt of problems. If we select a too low level of data 
disaggregation most intra-industry trade registered may be due to categorical 
aggregation. On the other hand, when the selected level of disaggregation is too 
high, we run the risk of losing the economic meaning of the subgroups employed. 
Some products with similar factor requirements may appear in different groups 
and part of the recorded inter-industrial trade would be, in fact, intra-industry 
trade. In most empirical works on intra-industry trade it is considered that the 3- 
digit SITC level is a reasonable approximation to the economic concept of 
“industry”. 

With the purpose of assessing the influence of categorical aggregation on intra- 
industry trade measurement at a particular level of data disaggregation, several 
procedures have been proposed. The first of them consists of performing 
measurement at higher levels of Statistical disaggregation. Among other 
shortcomings, it has the problem of dealing with the elevated volume of 
information needed to compute indices at high levels of data disaggregation. The 
second procedure, consisting of computing indices based on different 
Classification Systems, is similarly arduous. Moreover, the problem of selecting a 
concrete level of disaggregation to perform the measurements still persists. 

Finally, the third procedure involves the computation of an adjusted intra-industry 

* 

trade index: Bj . This adjusted index is giving a more convenient measure of 

intra-industry trade than index Bj in the presence of categorical aggregation. Yet, 

if categorical aggregation is not associated with opposite sign trade imbalances in 
$ 

the subgroups, thenBj would be a downward biased intra-industry trade measure. 

In the absence of information about the existence of categorical aggregation, it is 
not possible to know which of both measures is more suitable. 

Finally, it is worthy to mention that none of the mentioned procedures gives 
conclusive results regarding the determination of the aggregation error ’s absolute 
relevance. However, it is convenient to employ at least one (or more) of these 
methods for indicative purposes, in order to reduce, in part, the arbitrary 
indiscriminate use of a specific level of data disaggregation to compute intra- 
industry trade indices. 
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9.4 

Summary and Conclusions 

In this chapter we have tried to carry out a critical review of the main indicators 

used for the purpose of intra-industry trade and specialisation measurement. 

Moreover, we have discussed the two main problems related to such 

measurement: trade imbalance adjustment and categorical aggregation. 

The analysis performed may be summarised in the following points: 

i) Two different criteria have been identified conceming intra-industry exchange, 
giving rise to two families of indices (not extrictly comparable): the degree of 
overlapping in trade flows and the degree of similarity in trade pattems (or 
relative structure of imports and exports). Among the indicators ^elonging to 
the first group, the indices proposed by Grubel and Lloyd (1975) seem to be 
the more suitable. Conceming the second^roup of indicators, it can be 
concluded that the index of Michaely (1962) is the most convenient measure 
when it comes to measuring the degree of similarity in trade pattems. 

ii) Regarding the procedures aimed at correcting the distortions introduced by 
trade imbalance in intra-industry trade measurement, a common feature to all of 
them is the lack of theoretical support to justify the need for adjustment. It is 
thus preferable to use unadjusted indices for carefixlly selected periods of time. 

iii) The problem of categorical aggregation remains unsolved. Although there is 
no definitive solution, several methods have been developed to test for the 
influence of categorical aggregation. With the possible exception of the product 
regrouping method, these tests are easy to perform and, even if not conclusive, 
they have an indicative character, which allows us to eliminate in part the 
arbitrariness of the use of a specific level of data disaggregation when 
measuring intra-industry trade. 

iv) Two issues have dominated the literature on intra-industry trade measurement 
in the last years: marginal intra-industry trade and disentanglement of vertical 
and horizontal intra-industry trade. 

The measurement of the so-called marginal intra-industry trade is related to the 
analysis of adjustment costs associated to trade liberalisation. The indices 
developed are based on the assumption that when assessing the relevance of 
intra-industry trade during an adjustment process, it is the proportion of intra- 
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industry trade in new generated trade what matters, rather than the level of 
intra-industry trade. 

On the other hand, the distinction between vertical and horizontal intra-industry 
trade is an important one as the determinants of both may differ. Diverse 
indices, which disentangle horizontal and vertical intra-industry trade, have 
been developed and used in empirical analysis. 

Finally, it is important to mention that recent work on this subjects addresses 
the issue that the indicators of marginal intra-industry trade so far developed 
may be underestimating the extent of intra-industry trade, as they identify intra- 
industry trade with horizontal trade only. In other words, they cannot 
distinguish between inter-industry trade and vertical intra-industry trade. A 
method to deal with this question can be found in Thom and McDowell (1999). 
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10.1 

Introduction 

One of the features which characterises recent international trade flows is the 
prevalence of intra-industry trade. This phenomenon may be defined as 
simultaneous imports and exports of products belonging to the same industry and 
constitutes an important part of total world trade. The lack of concordance 
between the predictions of the traditional trade theory and the empirical evidence 
found has given rise to the appearance of new theories, which try to explain the 
structure and pattems of international trade by introducing imperfect competitive 
markets where the existence of scale economies and diversity of consumer 
preferences constitute essential ingredients. 

The assessment of the relevance of the new theories of trade as opposed to the 
traditional theoretical approach, when it comes to explain real trade pattems, is 
essentially an empirical question. Nevertheless, there are two problems (which 
remain unsolved) which influence the results obtained when measuring intra- 
industry trade and specialisation. The first problem refers to the very existence of 
the phenomenon and is related to the definition of “industry” and the selection of 
the level of data disaggregation more appropriate to study such a phenomenon. 
The second problem, closely connected to the former, comes ffom the objective 
difficulty of finding a convenient quantitative measure. 

This chapter is devoted to the analysis of one of the main obstacles faced by 
empirical treatment of intra-industry trade: the problem of categorical aggregation. 
The aim of this study is to effectuate a measurement of Spanish intra-industry 
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trade for several groups of products and country areas and to evaluate the 
influence of categorical aggregation applying several tests. With this purpose, we 
introduce first in section 10.2 the methodology employed. We offer next in section 
10.3 the main results obtained. Finally, in section 10.4 we present a summary and 
the main conclusions reached regarding this aspect of intra-industry trade 
measurement. 

10.2 

Methodology 

10.2.1 Measurement of Intra-industry Trade 

Intra-industry trade may be defined as simultaneous imports and exports of 
products belonging to the same industry. As stated in previous chapter, a variety 
of intra-industry trade indicators has been introduced in the literature. Among 
them, we have employed the most widely used, which is the Grubel and Lloyd 
(1975) index. This index measures the share of intra-industry trade in total trade 
flows of an industry j\ which is computed by the following expression: 

(Xj+Mj) — |Xj-Mj| i |Xj ~Mj| 

J (Xj+Mj) (Xj+Mj) 

where 0 <B, < 1. 

The average of indices Bj for a set of industries (j=l, 2, 3, n) is computed as 
follows: 



_ n n (X i + M i ) (X | + M j ) - X j - M j 

Bi = Zw r B i =z - - - 

* ^i(x j+Mj) ( X i +M ^ 

j=i 



S(Xj +Mj)— X Xj — Mj V X j - M j 



j=l 



j=l' 
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I(Xj+Mj) 
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The problem in practice lies in finding out if a group of products included under a 
concrete level of aggregation by a trade Classification System can be considered as 
homogeneous ffom the perspective of intra-industry trade. That is to say, if this 
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Commodity group constitutes an industry 1 . This gives rises to the denominated 
categorical aggregation problem, which will be treated in next section. 

10.2.2 Categorical Aggregation: Definition of the Problem and 
Assessment Methods 

One of the most important problems of intra-industry trade’ s measurement is the 
unknown influence of categorical aggregation. The identification and the 
measured levels of intra-industry trade depend on the degree of homogeneity of 
the products classified under the same group in trade statistics and the nature of 
such homogeneity. The problem brought about by categorical aggregation on 
intra-industry trade measurement emerges when heterogeneous products 2 are 
classified under the same group in trade statistics. The existence of bi-directional 
trade flows for a commodity group will identify “genuine” intra-industry trade and 
“spurious” intra-industry trade or Statistical aggregation 3 . 

While the most convenient method to eliminate this problem seems to be the 
regrouping of data in a way that all resulting groups correspond as close as 
possible with the theoretical construction of an “industry”, the problem is that 
there are multiple criteria of accomplishing that regrouping. And even if a specific 
criterion is systematically used, reclassification is an extremely arduous task . 

International trade data may be classified according to different criteria. One of 
the most extended Systems is the Standard International Trade Classification 
(SITC), in which products are grouped at various levels of aggregation identified 
by digits. The alternative to reclassification has been the election of a specific 
level of data disaggregation in official nomenclatures as the best approximation to 
the concept of industry. In this case, two problems have to be solved. First, the 
most convenient level of data disaggregation has to be identified 5 . Second, the 
influence of categorical aggregation for the selected level of data disaggregation 
has to be determined. 

Taking into account the difficulties of product reclassification using a 
systematic and coherent with the economic criterion procedure, several alternative 
methods have been proposed. Particularly, three alternative methods to assess the 
influence of categorical aggregation on intra-industry trade measures have been 
suggested : first, performing measurement at higher levels of Statistical 
disaggregation; second, executing diverse measurements based on different 
Classification Systems; and third, computing an adjusted intra-industry trade index. 
In this work, we have begun by performing intra-industry trade measurement 



Most studies made on this subject consider that the Standard International Trade Classification 
3-digit groups are a reasonable approximation to the economic concept of industry, existing a 
certain degree of consensus in this respect. 

2 

That is, products which do not belong to the same industry. 

3 

Simultaneous exchange of products classified under the same group but with different factor 

requirements. 

4 

An additional problem would be the Classification of pieces and components. 

^ That is, the level that more narrowly corresponds to the concept of “industry”. 

See Greenaway and Milner (1983). 
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using the 3-digit groups of the SITC. Next, we have used the first and the third 
methods commented above to assess the influence of categorical aggregation on 
the obtained measures. 

The first of these procedures 7 is based on the consideration that a substantial 
fall in average intra-industry indices when descending from a specific level of 
aggregation to a higher one may be indicating the presence of categorical 
aggregation problem. There is no specific Standard pattem to evaluate if the fall in 
the value of the index is significant and so the results of this test are not conclusive 
but have an indicative character. Usually, a fall in the value of the index of about 
20% when passing from the 3-digit SITC to the 5-digit SITC level of aggregation, 
while maintaining a reasonably high volume of intra-industry trade for the 5-digit 
level, is considered to indicate that the main component of recorded intra-industry 
trade is not categorical aggregation. Consequently, we have performed 
measurement at the 4-digits and 5-digits levels of data disaggregation and 
evaluated the fall in the 3-digits indexes of intra-industry trade. 

The other type of categorical aggregation test performed involves a more 
systematic way of evaluating the categorical aggregation problem associated with 
trade imbalances of opposite sign in subgroups of products with different factor 
ratlos and low degree of Substitution. If we considered commodity group j formed 
by m subgroups denoted by z, this procedure consists of computing an adjusted 
intra-industry trade index according to the following expression: 

m — ( X i j + M jj )- Y | Xjj - M jj| 

B * = IWij • B, = ^ * 

1=1 ZCXy+Mij) 

i=l 



where Wy = 
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So, the aggregation procedure to compute the intra-industry trade index for 
commodity group j consists of computing By indices for the subgroups i and 
averaging them according to the share of every subgroup in total trade of /. 

As commented in previous chapter, the adjustment is thus based on the 
assumption that categorical aggregation is associated with opposite sign trade 
imbalances in the subgroups i of which a commodity group j is composed. When 
this Situation exists, a measurement of intra-industry trade at a higher level of data 



Which has been the most widely employed. 
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holds, we have 1 > Bj > B j > 0 . 



This means that, the more trade imbalances of opposite sign cancel between subgroups, the 
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higher Bj will be related to Bj and, in case all trade imbalances in the subgroups have the 
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same sign, we have Bj = B j . 
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disaggregation generates a lower value of intra-industry trade and helps to correct 
for the influence of Statistical aggregation. 

Applying this method, adjusted indexes have been calculated for every SITC 3- 
digits group j consisting of m subgroups i (4-digits). 

10.3 

Main Results Obtained 

Table 10.1 shows the Grubel and Lloyd intra-industry trade indexes calculated for 
several levels of data aggregation for the 10 sections of SITC and for the 
aggregates of primary products, manufactured products and total products. The 
indexes have been calculated for several country aggregates: European Union, rest 
of world and total world. 

Looking at the 3-digits indexes we can underline the importance of intra- 
industry trade in Spanish trade flows, especially if we consider the exchange of 
manufactured products with the EU countries. According to these indexes, more 
than 60 per cent of total trade is of the intra-industry type. The share of intra- 
industry trade is particularly high for some sections, including beverages and 
tobacco, manufactured goods classified chiefly by material, machinery and 
transport equipment and Chemicals and related products. Additionally, we can 
observe that intra-industry trade seems to be more significant in the exchange of 
manufactured products than in primary products for all the country aggregates 
considered and that the proportion of this trade is higher for trade with the EU 
countries than for the rest of the world. 

As can be seen from table 10.1, the influence of categorical aggregation varies 
depending on the group of products and the country aggregates considered. 
Considering the narrowest groups of products (5-digits), we can indicate that the 
share of intra-industry trade is still relatively high in the exchange of 
manufactured products (58.0 per cent for total trade and 59.8 per cent for trade 
with the EU). These figures are high enough to sustain that observed intra-industry 
trade in this case is not only the result of categorical aggregation 9 . 

The proportional fall experimented by intra-industry trade index when passing 
from 3-digits to 5-digits reaches a value of 9.9 per cent for the case of exchanges 
with EU countries and 10.3 per cent for total trade, being lower, in general, than 
the one observed for other countries. This fact, together with the maintenance of 
high proportions of intra-industry trade even for the highest disaggregation levels, 
confirms the existence of “genuine” intra-industry trade and not only categorical 
aggregation for the 3-digits groups of SITC: 

On the other hand, comparing the effects of aggregation on the groups of 
primary products and manufactured products, we observe that the former shows 
the influence of this problem with higher intensity when we consider the 
exchanges with EU countries, while the opposite holds for trade with the rest of 
the world. The less intense falls in the indexes are to be found for the exchange of 
manufactured products with the EU. 



9 



That is, the heterogeneity of products classified under the same industry. 
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Table 10.1. Bj intra-industry trade indexes for several data aggregation levels. Spanish SITC 
trade data 1992 







EU-12 






Rest of World 






Total World 




SITC Rev. 3 Sections 


3 dig. 


4 dig. 


5 dig. 


3 dig. 


4 dig. 


5 dig. 


3 dig. 


4 dig. 


5 dig. 


0 Food and live animals 


0.357 


0.283 

(7.4)* 


0.243 

(4.0)** 

(11.4)*** 


0.377 


0.242 

(13.5) 


0.169 

(7.3) 

(20.8) 


0.415 


0.317 

(9.8) 


0.265 

(5.2) 

(15.0) 


1 Beverages and tobacco 


0.848 


0.164 

(68.4) 


0.115 

(4.9) 

(73.3) 


0.272 


0.247 

(2.5) 


0.116 

(13.1) 

(15.6) 


0.781 


0.229 

(55.2) 


0.152 

(7.7) 

(62.9) 


2 Crude materials, inedible, 
except fuels 


0.523 


0.447 

(7.6) 


0.350 

(9.7) 

(17.3) 


0.207 


0.165 

(4.2) 


0.127 

(3.8) 

(8.0) 


0.439 


0.374 

(6.5) 


0.281 

(9.3) 

(15.8) 


3 Mineral fuels, lubricants and 
related materials 


0.667 


0.457 

(21.0) 


0.406 

(5.1) 

(26.1) 


0.132 


0.117 

(1.5) 


0.082 

(3.5) 

(5.0) 


0.240 


0.193 

(4.7) 


0.181 

(1.2) 

(5.9) 


4 Animal and vegetable oils, 
fats and waxes 


0.539 


0.497 

(4.2) 


0.388 

(10.9) 

(15.1) 


0.251 


0.209 

(4.2) 


0.207 

(0.2) 

(4.4) 


0.401 


0.363 

(3.8) 


0.325 

(3.8) 

(7.6) 


5 Chemicals and related 
products, n.e.s. 


0.560 


0.537 

(2.3) 


0.473 

(6.4) 

(8.7) 


0.718 


0.557 

(16.1) 


0.429 

(12.8) 

(28.9) 


0.640 


0.588 

(5.2) 


0.507 

(8.1) 

(13.3) 


6 Manufactured goods 
classified chiefly by material 


0.750 


0.651 

(9.9) 


0.566 

(8.5) 

(18.4) 


0.580 


0.450 

(13.0) 


0.370 

(8.0) 

(21.0) 


0.732 


0.639 

(9.3) 


0.561 

(7.8) 

(17.1) 


7 Machinery and transport 
equipment 


0.707 


0.684 

(2.3) 


0.665 

(1.9) 

(4.2) 


0.572 


0.503 

(6.9) 


0.456 

(4.7) 

(11.6) 


0.696 


0.671 

(2.5) 


0.648 

(2.3) 

(4.8) 


8 Miscellaneous manufactured 
articles 


0.550 


0.519 

(3.1) 


0.478 

(4.1) 

(7.2) 


0.426 


0.347 

(7.9) 


0.304 

(4.3) 

(12.2) 


0.525 


0.483 

(4.2) 


0.446 

(3.7) 

(7.9) 


9 Commodities and 
transactions not classified 
elsewhere in the SITC 


0.032 


0.032 

(0) 


0.032 

(0) 

(0) 


0.030 


0.030 

(0) 


0.030 

(0) 

(0) 


0.034 


0.034 

(0) 


0.034 

(0) 

(0) 


Total primary products (0 to 4) 


0.465 


0.326 

(13.9) 


0.272 

(5.4) 

(19.3) 


0.228 


0.171 

(5.7) 


0.121 

(5.0) 

(10.7) 


0.390 


0.286 

(10.4) 


0.238 

(4.8) 

(15.2) 


Total manufactured products (5 
to 9) 


0.667 


0.638 

(2.9) 


0.598 

(4.0) 

(6.9) 


0.555 


0.460 

(9.5) 


0.398 

(6.2) 

(15.7) 


0.668 


0.623 

(4.5) 


0.580 

(4.3) 

(8.8) 


Total products (0 to 9) 


0.639 


0.583 

(5.6) 


0.540 

(4.3) 

(9.9) 


0.441 


0.359 

(8.2) 


0.301 

(5.8) 

(14.0) 


0.602 


0.544 

(5.8) 


0.499 

(4.5) 

(10-3) 



* Fall in intra-industry trade index 3-digits to 4-digits 



** Fall in intra-industry trade index 4-digits to 5-digits 

*** Fall in intra-industry trade index 3-digits to 5-digits 

Source: Own elaboration with trade data provided by Hamburg Institute for Economic Research 



Another fact to be mentioned is that the measures obtained for trade in 
manufactured products with non-communitarian countries are more affected (that 
is, they suffer higher falls) when passing to higher disaggregation levels than those 
corresponding to trade with EU countries. On the other hand, the influence of 
categorical aggregation is more significant in the measurement of trade in primary 
products with communitarian countries than with the rest of the world except for 





10 Measurement of Intra-industry Trade: a Categorical Aggregation Exercice 161 



section 0 (food and live animals). This fact suggests that trade with these 
geographical areas involves different types of products. 

Focusing on the effects that categorical aggregation has on the different 
sections of SITC, we observe that this problem does not affect equally to all of 
them. 

Within the category of manufactured products, section 7 (machinery and 
transport equipment) registers the lowest falls in intra-industry trade index when 
passing to higher disaggregation level. 

A case to be equally mentioned is section 5 (Chemicals and related products). 
When measuring intra-industry trade at the 3-digits level we obtain for this 
section, contrary to the general trend, an index of intra-industry trade for non- 
communitarian countries higher than the one obtained for trade with EU countries. 
However, when we carry out the measurement at a higher disaggregation level (5- 
digits), this characteristic stops being and exception. While in the case of EU trade 
this section is less affected by categorical aggregation, the fall experimented by 
the share of intra-industry trade with non-communitarian countries reaches 28.9 
per cent. This suggests, on the other hand, that intra-industry exchanges within 
this section affect different types of products depending on the group of countries 
considered. 

Conceming primary products two cases deserve attention. On the one hand, 
section 0 (food and live animals), for which a higher intra-industry trade was 
obtained for trade with non-EU countries, loses this character when carrying out 
measurement at higher levels of disaggregation. Contrary to what can be observed 
for the rest of primary products, this section is more affected by categorical 
aggregation in the exchanges with non-communitarian countries than with EU. On 
the other hand, section 1 (beverages and tobacco), which showed an unusually 
high level of intra-industry trade when measured at the 3-digits level, seems to be 
dramatically affected by categorical aggregation, experimenting falls in the 
indexes of 73.3 per cent (for EU countries) and 62.9 per cent (for total trade) when 
measured at the 5-digits level. 

With regard to adjusted intra-industry trade indexes test, table 10.2 shows the 
percentage of industries which present falls in intra-industry trade index lower 

$ 

than those mentioned when using adjusted index B j instead of index Bj 10 . 

The results of the calculations accomplished allow to state that categorical 
aggregation does not have a big influence in the measurement carried out, 
according to this test. 

From the information contained in table 10.2, we can indicate that more than 93 
per cent of industries present decreases in intra-industry trade lower than 25 per 
cent when using adjusted indexes for the case of total trade. This figure is of 92.66 
per cent for EU trade and more than 85 per cent for exchanges with rest of the 
world. Approximately 70 per cent of SITC groups (63 per cent for non- 
communitarian countries) present a decrease lower than 5 per cent in the share of 
intra-industry trade when using adjusted indexes. Even if we observe the number 
of industries in which intra-industry trade falls less than 1 per cent when using 



io 



The full series of indexes are available at Carrera-Gömez (1996). 
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adjusted indexes, it accounts for more than 60 per cent (52 per cent for non-EU 
countries). 



Table 10.2. Percentage of industries with intra-industry trade decrease after using adjusted index 



Percentage of industries 


Reduction of intra- 
industry trade 
measured by B, 


EU-12 


Rest of World 


Total World 


less than 25% 


92.66 


85.33 


93.44 


less than 20% 


89.19 


81.08 


89.58 


less than 15% 


84.94 


77.61 


83.78 


less than 10% 


80.31 


70.27 


78.76 


less than 5% 


70.66 


62.55 


69.50 


less than 1% 


60.23 


51.74 


62.16 



Source: Own elaboration 



Finally, we can mention that correlation coefficients between the two series of 
indexes (unadjusted and adjusted) have been calculated for all the country groups 
considered. The high value of the coefficients (higher than 0.9 in all cases) 
indicates that the order of industries under the criterion of intra-industry trade 
intensity does not vary substantially when using adjusted indexes. This suggests 
that they are likely to produce similar results to those corresponding to unadjusted 
indexes when being used as dependent variable in cross-section analyses. 

10.4 

Summary and Conclusions 

This work is devoted to one of the main controversies conceming intra-industry 
trade: the problem of categorical aggregation. This problem is caused by the fact 
that essentially heterogeneous products are grouped together in trade statistics 
introducing a bias in intra-industry measurement. The more convenient way of 
facing this problem could be the reclassification of data, building product 
categories as homogeneous as possible, (that is, according to what can be 
considered an “industry”). However, this regrouping is problematic due to the lack 
of a unique reclassification criterion, among other factors. 

In this study we have first performed a measurement of intra-industry trade at the 
3-digits level of the SITC, considering that this level is a reasonable 
approximation to the economic concept of “industry”. 

Next, with the purpose of assessing the influence of categorical aggregation on 
the measurement carried out, several procedures have been employed. The first of 
them consisted of performing measurement at higher levels of Statistical 
disaggregation (4-digits and 5-digits levels of the SITC). The results obtained 
confirm the maintenance of an elevated proportion of intra-industry trade even 
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when measured at the highest disaggregation levels for the case of manufactured 
products. The share of intra-industry trade in primary products is considerably 
lower, being worthy to mention that, in this case, measurement is not highly 
affected by categorical aggregation either except for particular cases. A featured to 
be highlighted is the fact that manufactured products show a higher influence of 
categorical aggregation in trade with non-EU countries than in exchanges with the 
rest of the world, while the opposite holds for trade in primary products. On the 
other hand, intra-industry trade in manufactured products seems to be less affected 
by categorical aggregation than intra-industry trade in primary products for trade 
with communitarian countries, as opposed to the Situation presented by trade with 
the rest of the world. This suggests that, as can be expected, trade with these 
geographical areas involves different kinds of products. 

The second procedure employed involved the computation of adjusted intra- 
$ 

industry trade indexesBj , which may give a more convenient measure of intra- 

industry trade than index B, in the presence of categorical aggregation associated 
with opposite sign trade imbalances in the subgroups. The results of the 
calculations accomplished allow to state that, according to this test, categorical 
aggregation does not have a big influence in the measurement performed. 

In conclusion, the evidence shown by the data and the results of the tests 
performed and, specifically, the existence of high volumes of intra-industry trade 
even at high levels of data disaggregation, demonstrates that this is a phenomenon 
that cannot be explained solely by Statistical aggregation. 

Finally, we have to mention that the problem of categorical aggregation 
remains unsolved. Although there is no definitive solution, several methods have 
been developed to test for the influence of this problem on intra-industry trade 
measurement. While none of them gives conclusive results regarding the 
determination of the aggregation error’ s absolute relevance, it is convenient to 
employ at least one (or more) of these methods for indicative purposes, in order to 
reduce, in part, the arbitrary indiscriminate use of a specific level of data 
disaggregation to compute intra-industry trade indices. With the possible 
exception of the product regrouping procedure, these tests are not difficult to 
perform 11 and, even if not conclusive, they have an indicative character, which 
allows us to eliminate in part the arbitrariness of using a specific level of data 
disaggregation when measuring intra-industry trade. 
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section Analysis 
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A wide ränge of studies devoted to recent evolution of international trade 
highlights the fact that a significant and increasing proportion of flows consists of 
bi-directional change of products belonging to the same industry. This 
phenomenon is known as intra-industry trade. 

Among the several interesting aspects arising from the study of this type of 
trade we can mention the following. On the one hand, it is difficult to find a 
satisfactory explanation for this phenomenon inside the traditional theoretical 
ffame. On the other hand, it is important to analyse the factors determining this 
kind of trade as they differ from the conventional determinants of international 
trade and, consequently, the effects of trade policies and structural adjustments 
may be different to those obtained under the traditional assumptions. The former 
question has given rise to the appearance of a wide ränge of explanatory models, 
what gives an idea of the diversity of conditions under which this phenomenon 
may arise . On the empirical side we find, on the one hand, a set of studies aimed 
at determining the relative importance of intra-industry trade through the use of 
several methods to measure it. On the other hand, we have a series of papers 
whose purpose is to identify the industry and country characteristics leading to 
intra-industry trade . 

This paper is among the latter and it is aimed at establishing the factors 
explaining intra-industry trade in the Spanish economy, from an industry 
perspective, given the increasing relevance that this type of flows has acquired in 



A comprehensive review of this literature can be found in Greenaway and Milner (1986), 
2 Tharakan and Kol (1989) and Carrera (1996). 
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P. Coto-Millan (ed.), Essays on Micro economics and Industrial Organisation 
© Physica- Verlag Heidelberg 2004 - 




166 



G. Carrera-Gömez 



3 

Spanish foreign trade . With this purpose, in section 11.1 we deal whith the 
measurement of Spanish intra-industry trade and calculate intra-industry trade 
indexes for several industrial sectors. Section 11.2 contains an analysis of the 
determinant factors of Spanish intro-industry trade including empirical contrast of 
a series of hypotheses related to industry characteristics which may influence the 
extent of intra-industry trade in Spanish manufacturing industrial sectors. Finally, 
section 1 1.3 synthesises the main conclusions of the paper. 

ii.i 

Spanish Intra-industry Trade Measurement 

A variety of intra-industry trade indicators have been introduced in the literature, 
the index proposed by Grubel and Lloyd (1975) being the most widely used for 
this purpose. The G-L indicator expresses the proportion of intra-industry trade in 
a given industry, which will be referred as i, as follows: 

|Xj - Mj| 

B i = 1-7 r 

(x ; + M ; ) 



where 0 < Bj < 1 . 

For a set of m industries an aggregated index can be constructed according to 
the following expression: 

m 

I 

B = 1 - — 

m / 

B 

where 0 < B < 1 . 

From a theoretical point of view, we can consider an industry as a group of 
products with similar factor requirements and a high degree of Substitution. 
Nevertheless, this is not easy to take into practice and the available international 
trade Classification Systems contain aggregations of products that can be 
heterogeneous from an intra-industry trade perspective . As a consequence, the 
main problem we face when trying to measure intra-industry trade is the difficulty 
of defming the industry from an operative point of view. The usual way of treating 
this problem consists of choosing a particular level of data aggregation, measuring 
and evaluating the influence of Statistical aggregation, often calculating intra- 
industry trade indexes for higher levels of data disaggregation. Most works made 
on this subject use the SITC 3digit groups as a reasonable approximation to the 
economic concept of industry, being possible to talk of a certain degree of 
consensus in this respect. 



pCi-Mil (Xi + Mj) 

\ _ m / X : 

Xi + MJ i=1 Z(Xi + MJ 



( 11 . 1 ) 



^ See Martin (1992) and Carrera (1996). 

This fact causes the denominated categorical aggregation problem. 
Standard International Trade Classification. 
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Table 11.1 shows the aggregated intra-industry trade indexes for Spain for a 
number of sectors which have been constructed from SITC 3-digit groups. 
Additionally, intra-industry trade indexes have been calculated for SITC 5-digit 
items. The ränge correlation coefficient between both series of indexes is 0,96, 
indicating that the ranking of sectors according to their intra-industry trade 
intensity has practically no variations between both cases. 

Table 11.1. Intra-industry trade indexes for Spain (1992) 


Sector 


Denomination 


IIT index 


1 


Solid fuel, coke, hydrocarbons, radioactive minerals and petroleum refineries 




0.226 


2 


Electric power, water supply and gas 




0.430 


3 


Metallic mineral products 




0.630 


4 


Non-metallic mineral products 




0.675 


5 


Chemicals 




0.616 


6 


Fabricated metal products (except machinery and transport material) 




0.806 


7 


Agricultural and industrial machinery 




0.632 


8 


Office and computing machinery, optical accuracy instruments and similar tools 




0.438 


9 


Electrical and electronic machinery and material (except Computers) 




0.585 


10 


Automobiles, pieces and accessories 




0.800 


11 


Other transport material 




0.648 


12 


Food, drinks and tobacco 




0.519 


13 


Textiles, leather and footwear 




0.570 


14 


Wood and cork 




0.541 


15 


Paper, graphic arts and Publishing 




0.716 


16 


Rubber and plastic products 




0.734 


17 


Other manufacturing industries 




0.563 



Source: Own elaboration with foreign trade data provided by Hamburg Institute for Economic Research 



11.2 

Analysis of the Determinants of Intra-industry Trade in 
Spain 

11.2.1 Hypotheses 

The bulk of intra-industry trade explanatory theories conforms an analytical 
outline from which we can extract some basic hypotheses conceming industry 
characteristics influencing the degree and extent of intra-industry trade. 
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One of the elements the literature is emphasising as responsible of intra-industry 
trade share of total trade in a sector is the differentiation of product. What we 
could call “genuine” intra-industry trade refers to the exchange of horizontally 
differentiated products, where the diversity of tastes and preferences of the 
consumers plays a significant role . On the other hand, technological 
differentiation of the product may lead to a specialisation of an inter-industry 
nature, where industries better endowed with technological knowledge produce a 
higher number of new products and enjoy a comparative advantage in this sense. 
Nevertheless, if the new products are included in international trade classifications 
together with the old ones, the exchange of both will appear as intra-industry 
trade. 

Consequently, the product differentiation hypothesis may be formulated on the 
following terms: 

• Hypothesis 1: The higher the degree of product differentiation in a particular 
production sector the greater the proportion of intra-industry trade in the sector. 

On the other hand, the presence of scale economies is considered as an important 
condition for the emergence of intra-industry trade in differentiated products. It is 
precisely this element that causes a firm finding more profitable to produce only 
one or a small set of varieties of the product to export, importing meanwhile the 
other varieties. This way the consumers are allowed to enjoy the benefits inherent 
to the increase in their ränge of election and a certain degree of intra-industry trade 
is generated. These issues allow us to state the following hypothesis: 

• Hypothesis 2: The share of intra-industry trade is higher in those industries 
with larger economies of scale. 

In addition to the aforementioned, several theoretical foundations may be found 
that induce us to expect a positive association between the participation of 
multinationals in a productive sector and the degree of intra-industry trade in that 
sector. The diverse cost minimisation strategies employed by this type of firms 
under certain assumptions may generate intra-industry exchange. This makes it 
possible to formulate the following hypothesis: 

• Hypothesis 3: The proportion of intra-industry trade in a particular industrial 
sector is larger the more intense the activity of multinational firms operating in 
that sector is. 

Finally, it has been frequently noted in the literature on intra-industry trade that 
this type of trade is more sensitive to trade barriers than traditional inter-industry 
trade, what allows us to suggest the following hypothesis: 

• Hypothesis 4: The share of intra-industry trade is higher in those industries 
facing a lower degree of trade barriers. 



On the other hand, it is often impossible to separate vertical and horizontal differentiation in 
empirieal work (mainly because both types of differentiation can be simultaneously found in 
products). In any case, a positive association between both forms of product differentiation and 
intra-industry trade can be considered. 
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1 1 .2.2 Data, Variables and Econometric Model 

The set of trade and induslry data used to construct the variables employed in the 
analysis were recorded in various Classification Systems or nomenclatures . 
Therefore, and due to the lack of official tables of correspondence among such 
classifications, it was necessary to carry out a previous task of conversion among 
the different nomenclatures with the aim of having the data homogeneity required 
for carrying out the work. The data were finally grouped together into the 
industrial sectors listed in table 11.1 (see section 11.1) and the test carried out is 
referred to a cross-section of these industrial sectors for the year 1992. 

The construction of the variables employed was made as follows: 

The dependent variable used in our analysis is the level of intra-industry trade. 
For every sector an aggregated index has been calculated according to expression 
(11.1) as it is shown in table 11.2. 



Table 11.2. Dependent variable 

Dependent variable 



Intra-industry trade share in total 
sector j trade 



Notes: 

j = 1, 2, 3,..., 17. 

i = 3-digt-SITC groups for every sector j 



Proxy 



Grubel-Lloyd (1975) aggregated index: 



. = 1 - 



Xy-M y| 



JjXij + M. 






As regards the explanatory variables, the following have been used: horizontal 
product differentiation, technological product differentiation, scale economies, 
involvement of multinational firms in the sector and level of tariff protection. The 
construction of the independent variables is summarised in table 1 1.3. 



The following Classification Systems have been used: 

-Stantandar International Trade Classification, Rev. 3. 

-Industrial Survey 1989-1992 sectors (National Institute of Statistics). 
-National Classification of Economic Activities (CNAE-74). 

-Tariff sector used in Melo and Mones (1982). 

-Industrial Situation Survey (Ministery of Industry and Energy). 
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Table 11.3. Explanatory variables 



Explanatory 

variable 


Proxy 


Horizontal product 
differentiation 


Advertising expenditures (AE) to sales 
(S) ratio 


4ES S,! 

J is, s i w l s ' S J 
1=1 


Relative AE 


AE; 

ZAEj 

j=l 


Technological product 
differentiation 


Research and development expenditure 
(RD) to value added (VA) ratio 


S RDr ^RDr VA r ) 

IVA r VA i K VA r VA jJ 

r=l 


Technological development expenditure 
(TD) to VA ratio 


ZTD r j D s f TD v A ) 

ZVA r VA i K VA r VA jJ 

r=l 


Scale economies 


Minimum efficient size (MES)* to cost 
disadvantage (CDR)** ratio 


SE = MES, , = £ s P|_. 

' CDR, J iti 1 Pj 


Involvement of 
multinational firms 


Foreign direct investment (FDI) 


FDIj = XFDIj. 

r=l 


FDI to VA ratio 


i FD, E FDI, 

FDIVA, =— = i- 

J s VA • 

Z VA r VA J 

r=l 


FDI to total direct investment (TDI) ratio 


FDI, 
FDIT, = J 

J TDIj 


Level of tariff 
protection 


Tariff protection (T) 


1 m ^ 1 m _ 

-ZT k -Mf -IT k -Mf 
T _ m k=l 3 _ m k=l J 

3 Mf + Mf U Mj 



Notes: 

j = 1,2, 3, 17 

i = 3-digt-SITC groups for every sector j 
1 = INE Industrial Survey sectors for every sector j 

r = CNAE-74 sectors (2 digits) for every sector j. INE Research and Development Survey 
k = Melo and Mones (1982) tariff sectors for every sector j 
T k = Average Common Extemal Tariff (Melo y Mones, 1982) 

Mj EU = Imports ffom the European Union 
Mj r = Imports from the rest of the world 
Mj = Total imports for sector j 

* MES is computed as the size (value of production) of the median establishment divided by the value of the 
sector’ s total production. By median establishment we mean the establishment which corresponds to the median 
of the accumulated distribution of the sector’ s production. 

** RDC is computed as the quotient whose numerator is the value added per employee in establishments smaller than 
the MES and whose denominator is the value added per employee in the remaining establishments. 



Once the variables have been constructed, we have to associate the intra-industry 
trade indexes with the set of industry characteristics in table 1 1.3 using regression 
analysis. But one of the difficulties we face when doing so is the limited ränge of 





